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16755 ‘HE GENERAL HYDROLOGIC SYSTEM ae 
KEY WORDS: ‘Data uniformity; Hydrologic models; Instantaneous unit 

hydrograph; Rainfall- runoff relationships; Runoff; Watersheds = 
ABSTRACT: The general Hydrologic System (GHS) model is considered. It is found 
that there are inconsistencies in equations for the instantaneous unit hydrograph (IUH) 
which were derived using the GHS Model. The existence of these inconsistencies is 
supported by an analysis of the data reported in the literature dealing with watershed © ti) 


simpler, and the above inconsistencies no longer remain. The reported data supports 


REFERENCE: McCann, Roger C. (Assoc. Prof, Dept. of Mathematics, Mississippi 
"State Univ., Mississippi State, Miss. 39762), ~ and Singh, Vijay P., “The General 

_ Hydrologic System Model,” Journal of the Hydraulics Division, ASCE, ‘Vol. 107, No. 
Basi Proc. Paper 16755, December, 1981, pp. 1581-1592 mB 


s; Channels 


“realistic, initial conditions are assumed then the GHS model becomes considerably : 


ABSTRACT: A theory of flow resistance for large- -scale roughness is developed, "thin be 
main points of which are supported by a flume study from which the outline of a 
process-based resistance equation is determined. Rivers with large-scale roughness have © 
steep slopes and depths of the same order of magnitude as the bed material size. i 
resistance depends on the form drag of the roughness elements and their disposition in — 
the channel. Theoretical considerations show resistance to be a function of Reynolds © 
number, Froude number, roughness geometry, and channel geometry. These processes 
are examined using the results of flume experiments based on five different roughness | 
beds and a wide range of flows. Semi-empirical analysis supports the theory and allows t 
the development of a theoretical power law resistance equation which accounts for 


most of the processes. A summary of the data is presented for use in other research 
eee. Bathurst, James C. (Sr. Scientific Officer, Inst. of Hydrology, — 
Wallingford, Oxon, United Kingdom), Li, _Ruh- Ming, and Simons, Daryl B., 
“Resistance Equation for Large-Scale Roughness,” Journal of the Hydraulics Division, — 
ASCE, Vol. 107, No. HY12, Proc. Paper 16743, December, 1981, pp. 1593- 1613 


" 6742 MEASUREMENT OF AIR IN FLOWING WATER 
KEY WORDS: Aeration; entrainment; Gamma _ rays; ‘Hot “film 
Hydraulic jump; Hydraulics; Two phase flo flow; ratio; 


_ ABSTRACT: Procedures are developed using gamma ray and hot-film anemometer 
| methods for determining void fraction and sizes of air bubbles entrained in flowing © 
water. The hot-film data are analyzed by computer. The gamma ray data are used to. 
_ determine the duration factor, a term accounting for the additional time that must be | 
; _ added to the time sensed by the hot-film anemometer to arrive at the true bubble 7 
_ duration. The duration factor is found to be a function of velocity. The procedures are — 
"developed using a hydraulic jump with Froude number 6. The continuity equation is — 
used in conjunction with void fraction and mean velocity data to compute the location , 


of air entering an ond the j jump. a 


REFERENCE: Babb, ‘Alan. F. Albrook Hydr, Lab., ‘and Prof. of Civ. 


Environmental Engrg., ‘Washington State Univ., Pullman, Wash. 99164), and Aus, 


Henry C., “Measurement of Air in Flowing Water,” Journal of the Hydraulics 
| Division, ASCE, Vol. 107, No. HY12, Proc. Paper 16742, December, 1981, PP. 1615- 
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16766 RADIAL SOURCE-SINK FLOW OCEAN 
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ys 


sources; Fluid Mixing, environments; Oceanic Proces 


wat 
| Solar energy; Stratified 


large sources of momentum discharging into the stratifies tropical ocean. The mixing © 
characteristics induced by such sources in combination with simultaneous sink flows _ 
into the plant intake are considered under the assumption of a discretely stratified 
ocean with radial symmetry about the OTEC plant. The problem is formulated | with» +: 
radial jet integral equations which account for the loss of jet momentum and for the 
pressure gradient, both arising from the sink flow action, and for buoyant damping of _ 
vertical entrainment. For extreme shallow upper layer conditions or for large 
momentum sources, the source-sink flow is predicted to break down into a 
recirculating eddy giving rise to undesirable temperature reductions of the OTEC plant — ae 
intake. The model predictions are verified in laboratory experiments which have a { 
somewhat different geometry but a similar dynamic interaction. Design considerations 
and estimates of maximum OTEC plant sizes are 
REFERENCE: -Jirka, Gerhard H. (Asst. Prof., School of Civ. 
Engrg., Cornell Univ., Ithaca, N.Y. 14853), Johnson, R. Peter, and 


ABSTR: ACT: Proposed ocean thermal energy conversion (OTEC) plants represent 


' - ‘Radial Source-Sink Flow in Stratified Ocean,” Journal of the Hydraulics Division, 


ASCE, Vol. 107, No. HY 12, Proc. Paper 16766, December, 1981, pp. 1631- 16500 


16744 Vv VERT ICAL ROUND BUOYANT JET IN SHALLOW WATER TER 


Ww ORDS: Buoyancy; Entrainment; Environmental neering; 
Hydraulic jump; Hydraulics; Stratified flow; Temperature distribution; 


ABSTRACT: The “stability and bulk “mixing. ‘an axisymmetric 
turbulent buoyant jet discharging vertically into a stagnant waterbody of large 
horizontal extent is studied theoretically and experimentally. A stable discharge + 
configuration is defined as one in which a buoyant surface layer is formed which 
“Tepon. radially from the source and does not communicate with the initial buoyant jet — 


region. On the other hand, the discharge configuration is unstable when recirculation 

cells exist around the jet efflux. A semi-empirical theory shows the discharge stability | 

is only dependent on the dynamic interaction of three near-field regions—the buoyant 

jet region, surface impingement region, and radial internal hydraulic jump region. A 

series of laboratory experiments were performed with a half-jet, using heat as source of | 
“buoyancy, and inserting a plane of symmetry in a large model basin. The experimental | 
_ results are in good agreement with the theoretical predictions, = as regards the 7 
Stability criterion and the near characteristics. 


REFERENCE: ‘Lee, Joseph H. (Asst. Prof., Dept. of Civ. Univ. of 


Delaware, Newark, Del.), and Jirka, Gerhard H., “Vertical Round 1 Buoyant Jet in 
Shallow Water,” Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY 12, 
Proc. Paper 16744, December, 1981, pp. 1651-1675 


16750 TRANSVERSE BED SLOPE IN ALLUVIAL CHANNEL BENDS 
WORDS: Alluvial streams; Channel bends; Channels; Fluid 
ics; Riv Ww ane 
Hydraulics; Rivers; Sediments; Shear stress; Water flow; aterways wok 


| . steady-state transverse bed profile in channel bends, formulated an alternative 


ABSTRACT: The study has reviewed existing analytical models for the prediction of 
model, and compared their predictive capability by laboratory and field data. The \ 


1 | Gatainme model, which is based on a theory for the critical conditions for sediment 


entrainment, gives good agreement between measured and predicted bed slopes. The | 


' - study shows that the key parameter in all available models is grain size. The prediction — 


of the complete bed topography of a given bend for a given flow deta requires 
knowledge of the spatial t f th ss 

REFERENCE: Odgaard, A. Jacob (Asst. Prof., Inst. of Hydr. Research, Univ. _ of lowa, 
lowa City, lowa 52242), “Transverse Bed Slope in Alluvial Channel Bends,” Journal of © 


the eta Division, ASCE, Vol. 107, No. HY 12, Proc. we its 16750, December, — 
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16760 SHEAR DISTRIBUTION AT CHANNEL CONSTRICTIONS 


KEY WORDS: Boundary shear; Channels; Constric 

distribution; Velocity; Velocity distribution = 
ABSTRACT: A Comprehensive laboratory study, employing an appropriate range of 
flow Froude Number and contraction ratio, was undertaken for curvilinear and 

_ accelerating flow occuring against the advancing face and through the gap at end- 7 
dumped channel constrictions. Pertinent results on flow pattern, velocity distribution ra 
and boundary shear distribution, are reported and discussed. The locus of the shear 

: _ maxima, which attains its highest point value on the advancing face of the constriction _ 
is observed to essentially overlap the path of the high velocity filament; a pattern that 

4 is similar to that found in a channel bend. Results show that which an increase in the 
contraction ratio, the corresponding increase in boundary shear over that in the ‘i 
uncontracted channel is exponential, although the increase in velocity is only nominal. 

A followup laboratory study on shear distribution at over-bank channel constrictions — 

REFERENCE: Das, Bishnu P. (Superintending Engr. , (Civ.) Dept. of Irrigation, 

Orissa, India. ), and Townsend, Ronald D., “Shear Stress Distribution at Channel 

_ Constrictions,” Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY12, Proc. 


16763 CONFORMAL MAPPING FOR CHANNEL JUNCTION FLOW a 
WwW ORDS: _Channel bends; Channel flow; Channels (waterways); 
Conformal "mapping; Free Hydraulics; 3 Hydrodynamics; 


¥ ABSTRACT: An analysis of the flow | atid at rectangular channel junctions with __ 
combining flow using the conformal mapping technique has been presented. Equations 
for the location of the Stagnation point and the free-streamline for different ee 
conditions at rectangular channel junctions have been developed. These equations have — 
been solved numerically by using Simpsons rule. It is observed that for each channel 
junction there is a unique value of discharge ratio corresponding to which as 
Stagnation point coincides with the upstream corner of the junction. The value of the — 
discharge ratio depends on the junction angle as well as on the bed-widths of the 
channels. It is also indicated that the zone of stagnant fluid between the free- 
streamline and the channel on the junction angle and the location of 


_ REFERENCE: Modi, Pashupati N. (Reader, Dept. of Civ. Ener. Malaviya Regions 
ok Coll.,Jaipur, India.), Ariel, Pitamber D., and Dandekar, Mohan M., ST 


ware 


_ Mapping for Channel Junction Flow,” Journal of the Hydraulics Division, ASCE, Vol. 
No. HY 12, Proc. Paper 16763, December, 1981, pp. 1713- 
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_ In accordance with the October, 1970 action of the ASCE Board of Direction, which stated stated — 
% that all publications of the Society should list all measurements in both U.S. Customary and 
_ (International System) units, the following list contains conversion factors to enable readers — 7 
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Roger C. - McCann and Vijay P. ‘Singh,? A. M. ASCE 


‘The General Hyd: Hydrologic System (GHS) model i by Chow and 
Kulandaiswamy (4,5,7) has been reported and utilized in several studies (2,3,6,8, — 
%, 10,11,12) dealing with watershed runoff ‘modeling, and flood routing i in open 
channels. In its general form, the GHS model, as expressed by Kulandaiswamy 
ar’* 


surface renoff; I = rainfall excess; - time; Q/dt' = ith 
derivative of Q with respect to ¢; dla’ | = jth derivative of J with respect 
to t; a, = ith coefficient associated with (i + 1)st time derivative of Q; b, 


= jth coefficient associated with (j + 1)st derivative of I; M = nonnegative 
integer; and N = | Ronnegative integer with M coefficients a, and 


watershed runoff data, Chow and Kulandaiswamy (6) conceded that M and 

N could be taken as 2 and 1, respectively. They reported that higher derivatives — 

of J and Q were insignificant for the cases which they studied and could, therefore, 
be ignored without sacrificing accuracy of the model. Further, the coefficients, 
and b,, were assumed to be either constant, oF functions of some 


the GHS was reduced to ering. 


which was further written, using linear notation, 


_ "Assoc. Prof., Dept. of Mathematics, Mississippi State ‘Univ., Mississippi Sta State, ‘Miss. 
af ? Assoc. Prof., Dept. of Civ. Engrg., Louisiana State Univ., Baton Rouge, La. 70803. 
Note. —Discussion open until May 1, 1982. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, — 
ASCE. Manuscript was submitted <or review for possible publication on July 28, 1980. 
er paper is part of the Journal of the Hydraulics Division, Proceedings of the American 
of Civil Engineers, ©ASCE, Vol. 107, No. HY12, December, (1981. 
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Bae between = numerator and denominator. The instantaneous unit ell 

7 graph (IUH), hereafter denoted by u(t), is a solution | of Eq. m when J is the 

_ Chow and Kulandaiswamy (4,5 7) derived i instantaneous unit hydrographs using 
Eq. 3 for four different cases depending upon the nature of the three roots 
of the denominator in Eq. 3. These cases are: (1) The three roots are real : 
and unequal; (2) all are real with two being equal; (3) all are real and equal; 
and (4) one is real and two are complex conjugates. We reproduce their solutions 


for each of the four cases with m, n, and 4 ina the three roots. 


Case 2: m, n, pi Real and and =p. The IUH, u(t), 
u(t) = *(Ae™ - 


inwhich C= es tow 


—mn 


A= 


>— bom b, m ) 


&§ 
Ez 
(5b) 
(6a) 
(6b) 


AB 2c 


I. The TUH, u(t), 
™+ev(B Cos wt wt (Ja) 


(7) 


“What 2 are the initial conditions need. to solve Eq. 3? No such conditions cp 
given (4,5 57). How was the IUH obtained?» A close examination reveals that — 
there is a mathematical inconsistency in the derivation of Eqs. 47 which does: 
not seem to have been observed previously (2,3,6,8,9, 10,11, 12). This inconsi - 
_ tency, which may lead to negative IUH ordinates, is: related to the initial conditions 
to obtain a unique solutiontoEq.3. 
In the following sections it is shown that if realistic initial conditions are 


assumed, the inconsistency ¢ can n be removed, | a realistic IUH can be cody wil 


7 vanish, This considerably simplifies the GHS in Eq. de 


Comments on So.ution or Ea. 
_ The question is: How was a unique solution for Eq. 3 obtained? In Chow 
a. and Kulandaiswamy (5) there are no conditions, initial or otherwise, which 
allow for the determination of a ‘unique solution. However, before giving the 
a IUH and the unit step response (USR), designated as } U(t), Kulandaiswamy 


dQ/dt, d?Q/dt’, 


‘ 
7 The initial conditions used are: the ‘System | is rest to start with, and 


system is of the P(D)/R(D), ‘the above quantities ; are not 


necessarily zero just after the input is applied. 
_ It may be noted that the form P(D)/ R(D) referred to here is _ > Se 


Thus, the initial conditions associated with Eq. 3 must be Q(0) = 


'(t). Thus, regardless of which function “satisfies 
= the initial conditions, we must have u(0) = U’(0) = 0. These initial values — 
not correspond to the the of the given by Chow and 


| | HYDROLOGIC SYSTEM MODEL 1583 
q 
| 
y 
| 
| | | 


‘DECEMBER 


3 
using 4-7 shows the 


In the first case, when m m p, u(0) = +mnj -b ,/a, 
_ 2. In the second case, when m, n, p are real and n = p, u(0) = = —b,/a,. " 
a In the third case, when m, n, p are real and m = n = p, u(0) = b »/@>- 
4, In the fourth case, when m is real, and n and p are complex conjugate 4 
of the 47 rainfall excess-runoff patterns on seven different watersheds analyzed — : 
_ in (5,7) for the GHS model. Therefore, this result is not consistent with what 
« the preceding quote implies. Moreover, if b, 4 0, the solutions in the four 
Since the roots of m, n, and p are functions of a,, a,, and a,, the — 
value of the IUH is ‘not a constant function of the coefficients in Eq. 3. There 
is no apparent reason why this should be the case. 5 a 
2 For the sake of argument, suppose that u(0) # Oi. e., suppose b, # 0. , 
m,n, and p are distinct real aumbers which are close to some number q; 
then 1 + a,D + a, the + =a, (D + m\(D + + p) and u(0) = 
=b,/a,. However, if m m=n p= q, then 1 + a,D + a,D? + a,D° = 
a,(D + q)’, and u(0) = / by. ‘Thus, if u(0) 4 0, the IUH not 
continuously on the coefficients a, and a,. We believe that it is unreasonable 
to assume that the problem of computing the IUH is ill posed in this manner. © 
Therefore, we consider u(0) = 0 to be a reasonable initial condition for the 
- For brevity and purposes of conserving space, we will hereafter restrict our 
_ analysis to the first case, m #4 n # p. Not only is this case the most important © 
of all, but the ensuing analysis, and the conclusions reached for this case, 
will also hold for the remainder of the cases, 
‘The question now is: If the preceding initial conditions were not used oo 
determine the IUH or USR, what conditions were used to obtain a a unique “ 
solution of the differential equation in Eq. 3? No method of solution is indicated s 
_ (5). However, it is indicated that the method of partial fraction was used. Using _ a 
in in which | A, , B, Car are one in Eqs. 4e- d. We nc now consider the case it) = 


= —mnp 


| 


It is that the of the Heaviside function 


158 

i 
7 


on (-~, oo) with £(0) = 1, 1, then F(8(0) = (0). Thus tones 
‘Se e™8(t)dt = = + +) Se = H(t): + Ch up 


 « Pre 


1 which c CL = arbitrary constants. follo lows. from that at for 


‘some of the Comparing Eq. 13 with the 10H 


given in Refs. 5 and 7, and reproduced in Eq. 4a, we conclude that c¢, = 


c, = c; = 1 or, equivalently, C; = 0. Without initial conditions, or 


FIG. 1, -—Typical IUH Obtained from GHS Model of Chow and Kulandaiswamy (4,5,7) 


some other conditions, there is no way to justify this choice of these constants. 
_ Furthermore, Kulandaiswamy (7) graphs his IUH for a number of rainfall 


we will consider only ‘the graphs ir in Ref. 7. In each « case > the graph ‘resembles 
Fig. 1. Let 7,, and ¥,, T, < T,, be the times at which the IUH is zero. 
From the graph in Fig. 1 (and those itis clear that 

u@d<0 . 

is physically unrealistic. ‘Moreover, an examination of the IUH graphs 
: given by Seneeenamy (7) indicates that T, — 7, can be as high as approx — 


Sani of T,, IUH peak time. Thus, the period of negative [UH ordinates is 


significant in in several cases which they studied. 


It may be appropriate to remark that, black-box methods of 
“ten identification (6), it is sometimes found that the IUH, determined from 
a given rainfall excess-runoff event, has some of its ordinates as negative (). 
If this same is used in conjunction \ the ‘same rainfall excess, ‘the 


4 
H(t 
: 
i 
| 


on is identical to the given renoff hydrograph. The 1 negative 
= may arise from one or more of the reasons: (1) Errors in data; (2) tie 
‘a errors in ‘model approximating | the Prototype system; and @)¢ errors in model 
cannot be used for other ‘rainfall excess events. in ‘this respect, it 
be worthwhile to quote Boneh Golan (1): te 
The answer to the question posed in the title of this paper is yes, provided — if 
_ that some restrictions are imposed on the rainfall €XCess which is the 
- input to the surface runoff system. Obviously impulsive inputs are. not 
4 admissible when the instantaneous unit hydrograph (IUH) has negative - 


ordinates. Since the IUH is the direct surface runoff response (known 

as the ‘output’) to an instantaneous rainfall burst of unit volume and 
q such runoff, by definition, cannot be negative. Is an instaneous unit — 
hy drograph with negative ordinates (1UH WHO) desirable? Here the answer 

Be. is no, because of the restrictions that have been imposed on the input. Pat 


Conomions ano Coerricients ano b, 


_ A simple, yet realistic initial condition can be assumed as u(0) = 


q used in the solution of Eq. ‘the coefficients and | b, vanish, as will 


be clear from the ensuin 
a Suppose that there is a aguas octets of Eq. 3 which we rewrite as. 


"+a, Bo u ‘+ u = by - b, 38 

subject | to u(0) = «wan (0) = a; -u”(0) = B for some numbers a and B. If 

we multiply each side of Eq. 12 by e~”, and integrate with ‘Tespect to to ¢ from | 


to in which > | 0, then we have 


“dt =1+b s ss 


since | (t) f(t) dt = (- ~1) ‘f'(). Since Eq. 16 for a Owe must 
have (where L is the Laplace transform of u): mel 

(a,s°+a,s* + a,s+1)L[u] -(sa+B +a) = 


in which m, n, and p are a roots on, = 
Using ‘the method of partial fractions, we can write 


4 


/ 
| 
i 
4 
0. Most 
E 
> 
ers 
q 


al 


for some choice of the number Ay B, and C, whenever the memerstor and 


A+ B+C=-b 


a Once A, B, and ‘eC have been determined, the solution u(t) of Eq. 15, with 


given initial conditions, can be obtained by inverting L as in 


u(t) = — (Ae™ + Be“ + Ce") ...... 


B+ 


‘ Thus, Eq. Is has a continuous solution with u(0) = 0 aa b, 
be: An analysis of the results reported in (5,7), for 47 rainfall excess- soaiaalt storm 
seven different watersheds, shows that ef vad 


30.8 =0.6 =<0.2 50. 1 30.08 =0.06 <0.02 an 


Thus, over half of the cases |u(0)| = 0. 2. This supports our argument ay | 
_ u(0) should be zero. Moreover, it must be emphasized that the coefficient 4 
was obtained by using a difference approximation to the second derivative of Ss 
E rainfall excess (5,7). Since rainfall excess is given as a piecewise constant ng 
(histogram), this difference approximation may introduce significant error ~ 
the calculation of b,. This accounts, at bon partially, for the relatively large 
To show that u '(0) = = and u"(0) =0, ‘the. following argument is 
The of the IUH is to write the solution of Eq. 2 
Since / lis given in histograph form, J (t) = some constant c >o on some 
interval O< t<t ,- Using Leibnitz’s rule, we have on the interval 0 < ¢ 


edt B= vale ow Inds ti ests sit) bam 


J 


ons is obtained: 
) 
| 
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Ot) = +e \. u(t . wor 


On the interval 0 < < 
+a au’ —1)+ u(t of one 


If we now use the identities = (- f(t) fort > 0, we 


a: 


=(- ™ — mp Be” — mnCe”' 

= 


BY 
in which a = a) B= 


—~0. Therefore, i = a or 6B = O. If B = 0, then by Eq. 31, - = 0. Thus, 


= @€636——sdIn none of the 47 reported cases (in none of the | Six possible models - 
each ca case) are a, and a, even close to Thus, we must have 
: whe 0, i.e. 


@ he 
(0) = = 0, @=0 


initial conditions are in the direct 
(DSRO) reported in Ref. 7. In all 36 ba the DSRO’s reported, the graph of — 
the DSRO’s resemble Fig. -2(a) near t = O. In at least 22 of the 36 cases, 
the graph of the DSRO has the oly as shown in Fig. 2(b) near t= 0. Let 
R(t) denote the | the DSRO at time Then i in at least 60% of the cases, lim = 0, 


to model the runoff near = 0, it that we should 1UH so 


4 that lim (t)=0 0 whenever is by Eq. 25. From . 26 we neve 
“ey 
4 


and in the remaining cases it appears that that w we hee have | lim R R*@) = = 0. Thus, 


(26) 1 

| 
a 
| 
i _ Inorder that volume be preserved under convolution with the IUH, we must 
= 1, we have from Egs. 22 and 21c: 
| 


HYDROLOGIC SYSTEM MODEL 
im Q” (aca... 
BS: Thus, choosing u(0) = = 0, u u’(0) = 0, and u u (0) = = 0 appear to be the logical — 
_ initial conditions for determing the IUH. Therefore, a = 0 and Bp =0. ity 
Furthermore, if a = 0 in Eq. 18, we may solve the system of Eqs. 2la-c 


A=——— ¢ 


hae — a, t 


(p —n\(p — m) . 


Te 
provides the, 
Study of 


2. 


we a, | 


FIG. 2. —Beginning of DSRO Obtained from GHS Model of Chow and Kulandaiswamy — 


| 
a 
ying 
= 
ok 
| 


= 0 only if b, = 
~ If b, = 0, the data reported (5, 7) al also indicate that b, should also be — iy 
001 <0.00 
—| <0. <0.07 <0. $0.03 = =0. ol s 50.005 0.001 <0.0005 


‘Thus, b, is very small when with. a, 
Susceptible to error when computed empirically. In fact, b is so small that 


is = 0, then from Eq. 36, u lead 


43, 0 146 38, 


*-2.09 


Upon simplification we obtain: 
| FIG. 3.—Comparison of IUH’s for Data of Kulandaiswamy (7, p. 149; Storm No. — 4 
7 1, N. Br. Potomac River Basin) for Two Cases: (1) When the Coefficients 5, and 
Are Not Zero, and(2)b,=6,=0 | 


To illustrate the difference between the two cases: Including (1) 5 b 


that the c case se with =), = 


reducing the accurac of the GHS model. 


e dra 


Y The following conclusions can be drawn from ee 


a va . The derivations « of the TUH ‘reported by Chow and Kulandaiswamy | (4, 5 D 
are inconsistent with regard to the initial conditions. That is perhaps “ om 


Pap 


An analysis of the reported in Ref. 7 that the second 


4. Their results (4,5,7) indicate that the period of negative IUH ordinates — 


be as high as 40% of UH peak time. 
35 Removal of the inconsistency in Ref. 4 (i.e., u(O0) = by = b, = 0) makes 
the GHS model much simpler and more attractive from a ‘practical standpoint. 
The work upon which this study is based was supported in part by fends 


provided by t the National Science Foundation under the project, “‘A 
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The following symbols are used in this paper: ated ¥ 
coefficient associated with derivative of rainfall excess; 
d‘[]/dt' = ith derivative of the quantity [ with: respect 
= number of derivatives of runoff; ad 
= root of the denominator in Eq. 3; — 
N= = number of derivatives of rainfall excess; 
on = root of denominator in Eq. 3; 
PB = root of denominator in Eq. 3; 
= Laplace transform variable; 
y = 


(t) ith derivative of u with respect tot; 


= real part of root; and rt 
= im imaginary part of of complex. root. 
where 


iol =» 


| 
4 
| 
| 
| 
a 
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RESISTANCE EQUATION FOR LARGE- SCALE : 


bess ate 


Simons, F. ASCE 


on 


received relatively little attention compared with channels with finer bed 
However, with the increasing development of mountain regions, especially in 
_ less industrialized countries, such channels are beginning to be featured in 
engineering projects and, consequently, there have recently been a few studies — 
of their general characteristics ¢,83, 16 ,20). This paper is a further contribution — 
to the subject and reports the results of a flume study of large-scale roughness _ 
carried out at Colorado State University (2). _ The aim of the study, which was 
intended to be of general scientific interest as well as of engineering hg 
to identify the processes affecting the flow hydraulics, construct a theory of 
j flow resistance, and test the theory by using the flume data to develop a theoretical 
q flow resistance equation. | Some of the work, particularly the description 
Vs flow pre processes, has been published elsewhere (3). This p paper, therefore, concen- 
4 4 trates on the flume work and the subsequent development of the resistance — 
equation and provides only a brief survey of the theory. A summary of the 
aa is provided which, it is hoped, may be of use to other research efforts. pritve 
The roughness i is large- -ecale if the roughness ¢ elements affect: the free ‘surface, — 
7 a condition requiring that the relative submergence, or ratio of depth, d, to | 
7 = =_= element height, be less than about four (7). The roughness is intermediate-_ 
scale if the relative submergence | lies between about four and 15, the latter 


& marking the boundary with small-scale roughness and the more familiar ? 


4 ‘Sr. Scientific Officer, Inst. of Hydro., Wallingford, Oxon, United Kingdom. 

? Assoc. Prof. of Civ. Engrg., Colorado State Univ., Fort Collins, Colo. sy 
Assoc. Dean, Coll. of of Engre., and Prof. o of | Civ. Engrg., Colorado State Univ., Fort 

Note. —Discussion open until May 1, ‘1, 1982. To ext extend the ‘closing | date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, — 
ASCE. Manuscript was submitted for review for possible publication on November 17, | 
se This paper is part of the Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, ©ASCE, Vol. 107, No. HY12, December, 1981. 


; ; 8 Boulder-bed channels with large-scale roughness, for which the size of the 
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sediments it is assumed that element height is given by the . short axis (4,12,13). 
However, for typical bed materials the equivalent son. in terms of the -” 


“type of flow which may be described by layer henry. For 


roughness | — <4 


n% of axes; and D, ofthe shor defined for the median axis. 
_ Under natural conditions, boulder-bed channels are likely to have slopes greater _ 
than 0.1-0.4% (8), but less than about 5-10%, and Froude numbers can attain 
supercritical values (17). . Entrainment of air can be significant (9,18) but more 
probably at very steep slopes, so it is not considered here. yt 


_ Previous studies have indicated that the flow resistance of large-scale roughness _ 
i is related mainly to the form drag of the elements and their disposition in 
the channel (1,7,10,13,20). Since the associated resistance processes are different 
from those for small-scale roughness, it is not possible to use the semilogarithmic 
_ resistance equations developed for small-scale roughness. Instead, it appears © 
that for large-scale roughness a power law is more suitable (1,6,13), while for | 
intermediate-scale roughness a semilogarithmic law (different from that 
In developing a theoretical process- -based equation for large-scale roughness, __ 
_ it is necessary to account separately for the processes of fluid mechanics related — 
to vastnenge number and Froude number which determine the drag of individual | 


determine the combined effect of the elements on the flow. Previous equations — 


Reynolds Number.—The form drag of an object v: varies according to whether a” 
7 the boundary layer on the object is laminar, turbulent, or transitional between Oe 
‘these states, and this is determined by the Reynolds number, U D,,/v, in which — 
O= mean velocity of flow, and v = kinematic viscosity. If the boundary — 
layer is laminar or turbulent, the drag coefficient does not vary much with 
~ Reynolds number except at extreme values. In the intermediate transitional — 
region, though, the drag coefficient falls as Reynolds number increases. wi tne inl 


| 
In many cases the transitional region extends o over a relatively narrow "range 


of Reynolds numbers (19) so the effect of Reynolds number is limited. With . 
ane -scale roughness, though, a combination of factors causes the transitional ioe 


| 
] 
| 
4 
| 


region to ts relatively wide, lying Seween: Reyaciée numbers of about 3 x 


10* and 2 x 10° (2,7). As this range is typical of that found in channel flows, 
the flow resistance could be significantly affected, generally decreasing as 
= 
‘Froude Number.—An element which protrudes through or nearly | through the 
- free surface causes distortions in the surface which represent an energy loss. 
The component of form drag generated in this way is referred to here as free ' 
surface drag and varies with Froude number and relative submergence. Pubcon 
Flammer, Tullis, and Mason (7) found that for a constant fguie senben, 
and for relative submergences greater than about 0.8, the drag coefficient 
_ decreases as relative submergence increases, until at a relative amped | 
of about 4 it is magnon. On the other hand, with relative submergence constant, 
a increase in Froude number causes the drag coefficient first to rise to a 
peak and then to ) decrease until, at a Froude | number of ; about 1.5, it is : again = 
negligible. The Froude number at which the peak occurs is about | for relative , 
- submergences less than about 0.5, and about 0.5 for relative submergences 
greater than about 0.8, a variation which is probably related to differential 
formation of hydraulic jumps around the element (2,3). The result of this 
_ dependency is that for a constant Froude number of a relatively low value, — 
and for relative submergences less than about 0.5 8, ‘it is the drag 
coefficient to increase as relative submergence i increases. 
‘Roughness Geometry.—The combined effect of the elements on flow resistance _ 
"depends on the proportion of the bed material which has a significant form | 
- drag, and this is determined by the geometry and disposition of the elements. 
Most resistance equations have accounted for roughness geometry with relative 
submergence only and have left unquantified the effects of size distribution, 
shape, and areal distribution. A more suitable parameter which would account - 
for all these factors, including relative submergence, is the effective roughness 
- concentration, A,, , defined as the ratio of the sum of the wetted frontal cross-sec- 


tional al areas, Ay, of oe elements, , to the area of bed on which they lies 


in which n = the number of significantly projecting elements on an area of 


ne 
to allow identification of such ‘elements, but which p> 
the number of significantly Projecting elements at a site can be isolated and 


an effective. roughness concentration, based on mean element properties, con- 
 Structed accordingly. These factors include: Relative (1); the 
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Theoretical analysis (2,3) indicates that the concentration 
should vary inversely | with relative submergence, the rate of change depending» 
on o. In addition, the absolute magnitude should be a function of roughness» 


shape and the number of elements. Therefore, assuming a power law representa- 
x 
in 1 which x x= a parameter r of roughness shape; ; and fn( )= a function (2,3). 
® flow resistance of the bed is directly proportional to this quantity. 


_ Channel Geometry.—In flows where the roughness elements occupy a a signifi- 7 


an proportion of the channel cross-section, the flow is funnelled between — 
the elements (10). The greater the proportion, the greater is the degree of funnelling 
and the higher is the mean - velocity of flow. This: effect can be “monitored — 


using the ratio of the total wetted cross-sectional aren of vail the elements at a 


Datum level Roughness: Water 


channel bed 


Aut we wetted roughness ere cr tional area 


‘the relative roughness area ea (Fig. 1). Research eleowhere has reached 
the same conclusion (20). 


- _ Since it may not be convenient to measure the relative roughness area directly, 


& would be useful if it could be calculated from the relevant parameters of 
tae shape and roughness geometry. Theoretical analysis (2,3) suggests that 


it should vary directly with the channel width /depth ratio, , and that the rate ‘ 


of change : should depend on the degree to which the roughness elements p project | 


7 into the flow at any given depth (quantified by the effective roughness concentra- 
tion). Still | using a power law, a functional relationship i is therefore es 


wd 

7 in which d = mean flow depth (2,3). This term to wary similar to ie 


| of channel Gerived in other studies (1,13). 


| 
| 
q 
section, A, depth 
| 
| 
/ 
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RESISTANCE E EQUATION 
_A variety of mathematical relationships between flow resistance and its — 
_ determining factors could be specified. However, previous ; research has found : *. 
that a power law is the most suitable for large-scale roughness (1,6,13). Conse- i 
= combining the foregoing analyses, and denoting resistance as in small- 


scale roughness with the function (8/f)'’* from the Darch-Weisbach equation ’ 
6, 11), ‘equation fc or -scale with steady, 


f 


in which f= the Darcy-Weisbach resistance coefficient; an = ait due 
to gravity; R = hydraulic radius; S = channel slope; R = + Reynolds number _ 
7 Mathematical descriptions of the resistance processes are not sufficiently 


advanced to permit theoretical quantification of Eq. 7. } Consequently, it 


to resort to more empirical based on data. 4 


ta were obtained for five different ned to poe beds using the steep 
stream flume at Colorado State University. i, 
Z Bed Material.—The five bed materials are classified as 0.5, 0.75, 1.5, 2. o- 


and 2.5 in. (respectively, 12.7, 19.05, 38.1, 50.8 and 63.5 mm); these figures | 


1 Distributions of Bed Materials 


Bed material Bed material Standard devi deviation | 


Median ALS 6.4 
Long (28.0 | 19.75 


Long “4 
Median | 58.0 | 54.25 to 


Short__ (44.0 | 39. 


| 
| 
a 

= 

| 

Short 23 | 80 | s2| 
80 | 40 | 35) 

re. | 
x 


ao cobbles. The 2.0 and 2.5 in. . materials were cobbles which had been | 

7 collected by hand. Size distributions for the long, median, and short axes of — 
the materials were determined by the Wolman method (21), and standard 

o deviations of the distributions were calculated using Eq. 4 for each axis in — 

— Flume. —The flume is of the open, recirculating type and is 9.54 m (31. 33 
ft) long, and 1.168 m_ (3.833 ft) wide. ‘Discharge is measured by a V-notch — 


weir in the sediment in eet tank at the end of the flume. Further details are 7 


charge, 
cubic 
- | meters 
second 
(2) 


0.00241 
0.01274 
0.03046 
0.05746 
0.07197 
0.00143 
0.00522 
0.01737 
0.03249 
0.04896 
0.00196 
0.00610 
0.01355 
0.03576 
0.06061 
0. 


aa 


meters 
4) 


0.0192 
0.0324 


square 
meters 
per 
second 
(x 10°) 
| 


173° 
1.173 


Froude 
num- 
ber 


Func- 
tion 
of 
effec- 
tive 
rough- 
ness 


con- 


tion 


he 
— 


Each roughness bed was constructed by ‘smearing | masonite boards with 
boards 
“were » then screwed to ) the bed « of the flume. Each bed was one ¢ element thick 
and, as there were no obvious patches of smooth board visible between the = . 
elements, had the approximate appearance of a natural river bed. Ree ae 
| The flume sidewalls (one perspex, the other painted wood) were not covered. = 
7” It was expected that, _ because of the extreme bed roughness and high ner 
_of width to depth, resistance would be little affected by the sidewalls and the = 
flow would be approximately two-dimensional. 


It is assumed (as for river sediments (4,12, 13)) that the height t of each indi, 


| 
| 
| 
| 
y 
veloc- Depth | Rela- | in bach 
Mean | of bed | rough- tance 
eT meters | depth | datum, | ness | coef submer- : 
7 slope second meters |A,/wd’ cient | d/Sso 
(0.02 0.146 | 0.0141 0.2659 0.392 | 1,044 | 2.436 | 0.300 
0.02 0.391 | 0.0279 0.1398 | | 206 (4806) 0.52700 
0.02 0.584 0.0446 | 0.0489 | 0.0874 0.883 0.205 | 7.695 | 0.747 
0.02 0.785 | 0.0627 0.0669 0.0635 1.001 10.804 | 0.942 
0.877 | 0.0702 | 0.0743 | 0.0540 1.056 | 0.143 1.038 
0.05 0.161 | 0.0076 | 0.0129 | 0.4104 0.590 | List | 1.312} 0.177 

0.05 6.296 | 0.0151 | 0.0197 | 0.2318 i= 0.679 | 2.606 | 0.336 
0.05 0.619 0.0240 | 0.0291 | 0.1739 0.246 | 4141] 0450 

0.0338 | 0.0385 | 0.1227 (1.429 | 0.196 | 5.827 | 0.592 

0.05 1.017 | 0.04612 | 0.0461 | 0.1068 | 1.600 7.103 | 0.669 

0.08 0201 0.0084 | 0.0124 | 0.3262 (16 0.701 1,303 i= 0.2270 
0.08 0.392 0.0133 | 0.0179 | 0.2532 1.083 | 0.545 | 2.299 | 0.307 

0.08 0.563 | 0.0206 | 0.0248 | 0.1701 | 1.205 | 1.252 | 0.408 | 3.552] 0439 © 
0.08 0.965 | 0.0317 | 0.0358 | 0.1149} 1.222 | 1.731 | | 5.466 0.600 
0.08 0.0424 | 0.0466 | 0.0903] 1.222 | 1.900 | 0.177 | 7.303 | ons 
0.08 0.0465 | 0.0507 | 0.0829} 1.222 | 1.927 | 0.172 | 8013 | 0.772. 
10.77 sq ft: 1m? = 35.3 cu ft 7 
= 


a have been the average of the long re so and aetien axes and is calculated — 


at each of three different slopes (2, 5, and 8%), thereby covering a wide range — a 
of conditions. At each flow, depth was gaged at a single cross section so that 7 
mean flow and channel properties could be calculated. A Single section was 


bn 


matic 
| 
oe! 


dad’ tive 
rough- tance 
meters | ness coef- 
second | second | meters A,,/wd' P cient 


Depth Rela- 


0.00580 
0.01181 
0.02482 
0.04047 
0.05348 
0.00381 
0.00843 
0.02037 
0.03333 
0.04586 | 0.904 
0.05460 | 0.979 
0.00207 | 0.186 
0.00631 


0.01007 | 0.430 
0.02825 | 0.807 
0.04518 § 1.032 
0.04879 | 1.064 


Note: | m = 3.281 ft; 1m? = 10.77 sq ft; im =35.3cuft. 


_ thought to be representative of the entire bed, because of the uniform ssienitiais 
of the bed, the large ratio of flume width to bed material size (which ensured 
that the flow pattern at the section would be determined by the average properties _ 
of the bed and not by individual and, possibly, unrepresentative elements), 
and the uniform, rectangular cross section of the flume. The section lay about = | 
7.1 m (23.3 ft) downstream from the flume entrance where it could be assumed | 
that, in a macroscopic sense, uniform flow prevailed. wT 
_ In flows with large-scale roughness, the cross-sectional area of flow is” 


‘significantly affected a the projections of the elements into the flow. Conse- 


‘ 
FABLE 3.—Flume Data Torv.75 in. Roughness Bed 
rough- 
Rela- con- 
Chan- submer-| tra- 
6 gence | tion 
slope d/Ss | 
0.02 0.0223 | 0.0282 | 0.2081 | 1.000 0.475 0.708 2.790 | 0.397 
0.0290 | 0.0349 | 0.1696 0.989 0.653 | 0.375 | 3.626 | 04800 
0.02 0.0439 | 0.0495 | 0.1146 | 1.000 | 0.738 0208 5.482 | 0.660 
0,02, 0.0591 | 0.0642 | 0.0801 0.989 0.770 0.270 7.383 | 0.846 
0.02 0.0698 | 0.0746 | 0.0641 0.989 0.792 | 0.255 | 8.728 | 0.975 
4 0.0141 | 0.0204 | 0.3052 | 0.989 0.619 0.269 
0.05 0.0199 | 0.0262 | 0.2411 | 0.989 | 0.822 | 0.591 | 2.484 | 0.349 § 
ae 0.05 0.0299 | 0.0360 | 0.1709 | 0.989 | 1.077 | 0.345 3.736 0.482 4 
0.05 0.0365 | 0.0426 | 0.1433 | 0.978 | 1.308 | 0.234 | 4.557 | 0.560 E. 4 4 
0.05 0.0434 | 0.0491 | 0.1156 | 0.989 1.385 }| 0.209 | 5.428 | 0.655 = 
: 0.05 0.0477 | 0.0536 | 0.1090 | 1.014 | 1.432 0.195 5.965 0.693 — 
0.08 0.0095 | 0.0159 | 0.4031 | 1.014 | 0.608 1.731 | 1.190 0.189 
0.08 0.0142 | 0.0211 | 0.3253 | 1.000 | 1.018 0.617 176 02s 86 
0,08 0.0200 | 0.0288 | 0.2222 | 1.000 | 0.970 | 0.680 | 2.505 0370 
—- 0.08 0.0299 | 0.0363 | 0.1742 | 1.000 | 1.489 0.289 | 3.743 | 0.477, 
0.08 0.0375 0.0435 | 0.1382 | 1.000 | 1.703 0.221 | 4682 | 05795 #8 § 
0,08 0.0392 | 0.0450 | 0.1285 | 1.000 | 1.715 0.218 | 4.905 | 0.605 7 ; 
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intervals of 12. 7m (0. 5 in. > det interval being equal to the smallest size of 
ik bed material. Total flow cross-sectional area, A (Fig. 1), was then calculated : 


“from the bed and free surface elevations by computer integration \ using the : 
trapezoidal rule. Cross-sectional area of the roughness, A,,, was calculated ror 
the bed elevations by the same method, the datum level for the channel bed 
being the top of the boards to which the elements were fixed (Fig. 1). Adding 


A and A, gave the total channel cross-sectional area, wd’, and thus the relativ 


Wetted perimeter is based on the een! level of the bed rather ‘than th 


> outline of the elements as they project into the flow. ‘Consequently, assuming — 


TABLE 4.—Fiume Data tert. in. Roughness Bed 
Kine- 
matic 
vis- 
cosity 
| of 
water, 
square - Rela- 
Mean - | meters tive 
meters | meter depth jatum, | | per f- | submer- 
| second ‘ gence 
slope | second meters | meters "| (x 10°) d/Sso 


0.00250 0.0297 . 
0.00868 0.0425 081 
0.01893 0.0548 . 053 
0.04352 0.0751 . 067 
0.06763 . 0. 0.0921 067 


0.08020 0.1009 1.140 
0.00181 0.0230 
0.00636 | 0. 0.0328 

0.01456 | 0. 0.0416 

0.03073 | 0. 0.0542 

0.06061 0.0716 

0.07421 | 0. 0.0792 

0.00389 .0124 | 0.0249 

0.01092 to. 0.0361 

0.02100 | 0. 0.0436 

0.03126 | 0.721 0.0805 

0.05498 | 0.971. 0.0603 0.322 | 2.550 | 0.511 
0.05574 | 0.883 | 0. 0.0657 | 0. 2 0562 


‘that the flow was two-dimensional, wetted perimeter i is equal to v width. Hydraulic: 

. radius is therefore equal to mean depth, d, given by the ratio of flow cross-sectional a 
: Finally, at the end of each gaging the water temperature was noted so that 
the kinematic viscosity could be obtained. The measured data and their derivatives a 


Posy 
are given in Tables 2-6. __ Saw 
Devetopment or Resistance Equation re 
Quantification of Eq. 7 was carried out, term by term, using the precepts — 


) 
: 
| | 
| 
| 
iis 0.274 | 2135 | 0.967 | 0.2330 
0.02 0.432 | 0.858 | 1.638 0366 
0.02 0595 0.483 2.276 0.471 
7 0.02 0.745 | 0.288 | 3.337 | 0.638 a 7 
0.02 0 0242 | 4225 0.766 
0.02 
0.05 
0.05 
0.05 
3 0.05 
0.05 > 
0.05 
0.08 
0.08 
0.08 
0.08 
» 
) 
iq 
— 
a 


Reynolds ‘Number.—Reynolds number, U Doo /v, v, varied from about 10° to 
> 4 x 10°. In this range, the boundary layer on an object is typically laminar 
-. the drag coefficient is independent of Reynolds number. Consequently, 
the flume flows the function of number in Eq. 7 is a constant. 
Effective Roughness Concentration.—Eq. 6 can be rewritten to . define the 
function of effective roughness concentration, in —_— measured oer 
as—bsothat 


Eq. there s be a unique between the effective roughness 


5. 5.—Flume Data for 2.0 in. Bed 


matic 
vis- 
cosity 


Dis- 


ers 4 Rela- 
cubic i Mean - | meters tive 
meters | meters depth submer- 
per d,in in second; i gence 
second | second | meters | meters J (x 10°) 4/Sso 


0.00329 

0.00837 

0.01158 
0.02541 

0.04047 

0.04949 
0.00329 
0.00713 
0.01413 
0.02068 
0.02941 

0.04368 
0.00247 
0.00565 
0.01077 
0.02187 
0.03249 
0.03724 


and relative submergence for each ghness bed, so values of vf 


bb calculated from Eq. 9 were plotted against d/S,. for each roughness bed — a 
(Fig. 2). In each case the relationship i is given by nthe le nolroow' aif oul! el 


in n which a and c = constants with 


‘From 5 be a function of roughness size distribution, so its 


| t 
> 
| charge veloc- Depth Relas | in | bach ness 
| con- 
 Chan- 
nel tion 
| ay 
0.02 0.100 | 0.0282 | 0.0505 | 0.4413 | 0.946 | 0.189 | 4463 | 0947 | 02200 
‘ ‘" 0.02 0.189 | 0.0378 | 0.0611 | 0.3814 | 0.937 0.311 1.654 1.269 | 0.281 
| 0,02 0.227 | 0.0436 | 0.0665 | 0.3443 | 0.937 | 0.348 | 1.323 1.463 0.324 
; 0.02 0.377 | 0.0578 | 0.0795 | 0.2735 | 0.937 0.500 0.639 1.938 | 0431 
; 0.02, 0.519 | 0.0668 | 0.0892 | 0.2511 | 0.937 | 0.641 | 0.390 2.241 | 0.483 
_— 7 0.02, 0.601 0.0705 | 0.0947 | 0.2553 | 0.926 | 0.722 | 0.307 2.367 0.48600 
0.05 0.132 | 0.0213 | 0.0442 | 0.5179 | 0.946 0.289 4.796 0.715 0.164 7 
] 0.05 0.214 0.0285 | 0.0513 | 0.4450 0.946 0.405 | 2.437 0.956 | 0.218 = 
: ; 7 0.05 — 0.337 | 0.0359 | 0.0575 | 0.3750 | 0.937 0.567 | 1.245 | 1.206 | 0.282 
0.05 0.431 0.0411 | 0.0633 | 0.3508 | 0.937 0.678 0.870 13799 | 033 
a. 0.05 0.542 | 0.0465 | 0.0688 | 0.3252 | 0.937 | 0.803 | 0.621 | 1.559 0.348 
0.05 0.643 | 0.0582 | 0.0788 | 0.2617 | 0.946 0.851 0.553 | 1.952 0.447 
0.08 0.162 0.0130 | 0.0411 | 0.6842 | 0.967 | 3.092 0.436 0.084 
+ 9,08, 0.205 0.0236 | 0.0505 | 0.5330 0.957 0.426 3.522 0.792 0.161 
ty 0.08 0.313 | 0.0295 | 0.0551 | 0.4646 | 0.957 | 0.582 | 1.892 0.989 0.28 @8©=6 = 
a 0.515 | 0.0363 | 0.0659 | 0.4483 | 0.946 | 0.863 | 0.860 | 1.219 0231  £&g 
| 0.08 0.637 0.0437 | 0.0747 | 0.4155 0.946 0.973 0.676 | 1.466 0.267 a 
; 7 0.08 0.712 | 0.0448 | 0.0701 | 0.3615 | 0.957 | 1.074 | 0.554] 1.502 22 8 
> 
= Im = 3.281 ft: 1m? = 10.77 
= 


Variation with: o was studied Fig. 3). The results support the theory 


c= 


4 = bea fu action of cham channel width and bed material cross-st 


erial m 

ra) 


a However, the relationship may not be ree since it is ile for a constant ¥ 


— 6. —Flume Data for 2.5 in. Roughness Bed — 


“water, 

meters 

meters | meters 

per | per 

second | second 


0.00993 
0.01671 
0.02799 
0.04110 
0.04967 
0.00369 
0.00855 
0.01282 
0.02176 
0.03403 | 
0.04896 
0.00397 
0.00605 
0.01128 
0.01775 

0.02737 
0.03319 

0.04485 


b varies directly with it is is proportional 


the effective concentration ) For this reason the negative 


4602 
q 
. 
BA: 
Chan- 
0.02 | 0.00409 | 0.138 | 0.0254 | 0.0567 | 0.5513 | 0.937 0.275 2.111 0.644 | 0.156 ‘ 
0.02 0.223 | 0.0381 | 0.0691 | 0.4489 | 0.957 | 0.365 | 1.201 | 0.964 | 0.234 — | 
0.02 0.301 | 0.0475 | 0.0777 | 0.3879 | 0.946 | O41 | 0.826 | 1.203 | 0296 
0.02 0.409 | 0.0586 | 0.0898 | 0.3469 | 0.946 | 0539 | Ossi | 1485 | 0.354 ie 
0.02 0.500 | 0.0703 | 0.1007 | 0.3021 | 0.937 | 0.602 | 0.441 | 1.780 | 042% ~ 
0.02 0.543 | 0.0782 | 0.1081 | 0.2761 | 0.937 | 0.620 | 0.416 | 1.980 | 0.476 — 
0.05 0.173 | 0.0183 | 0.0489 | 0.6266 | 0.937 | 0.408 | 2.399 | 0.463 0. * 2 
. 0.08 0.283 | 0.0259 | 0.0585 | 0.5575 | 0.957 | 0.562 | 1.268 | 0.655 | 0.153 — . 
7 0.05 0.342 | 0.0321 | 0.0635 | 0.4942 | 0.946 | 0.609 | 1.080 | 0.813 | 0.196 — q 
. 0.05 0.478 | 0.0390 | 0.0707 | 0.4490 | 0.946 | 0.773 | 0.669 | 0.986 | 0235 U7 
0.05 0.611 | 0.0477 | 0.0799 | 0.4034 | 0.937 | 0.894 | 0.501 | 1.206 - iP 
: 005 0.725 | 0.0578 | 0.0889 | 0.3497 | 0.946 | 0.962 | 0.432 | 1.464 | 0.350 i 
0.08 0.210 | 0.0162 | 0.0463 | 0.6503 | 0.957 | 0.527 | 2.309 | 0.410 | 0.101 — 
0.08 0.259 | 0.0200 | 0.0517 | 0.6141 | 0.946 | 0.586 | 1.863 | 0.506 0.120 
Z 0.08 0.374 | 0.0258 | 0.0575 | 0.5512 | 0.937 | 0.743 | 1.159 | 0.654] 0.156 
0.474 | 0.0321 | 6.0630 | 0.4911 | 0.957 | 0.845 | 0.896 | 0.812 01198 
0.08 0.592 | 0.0396 | 0.0705 | 0.4383 | 0.946 | 0.950 | 0.709 | 1.002] 0.24 #8 
0.08 0.669 | 0.0425 | 0.0740 | 0.4265 | 0.937 1.037 | 0.596 | 1.075 0.257 
0.775 | 0.0495 | 0.0810 | 0.3887 | 0.946 1.112 0.518 | 1.254] 0.299 
Note: | m = 3.281 ft; 1 m? = 10.77 sq ft; 1 
Eq. 5 are supported by the data and Eqs. 10-12 a 
ive roughness concentration 
7 
4 


“sign must ‘appear with b in , otherwise the variations of Eqs. 6 and 9 

Relative Roughness Area. —Fer a given depth of flow, the bigger the “i a 

roughness area, the greater is the funnelling of the flow between elements, a. 
a the bigger is the mean velocity of flow, U. It is therefore assumed that © 


the resistance function, 7), varies directly with A,/wd’. This 


an valor of &. Por 
FIG. 2.—Variation of Function of Effective Roughness Concentration with a a 
Submergence for Each ‘Roughness B: Bed (1 in. = 2! 25. 4 mm) be 


D MATERIAL 
= 
40.5 inch 


4 


06. 648 


— 


Distribution (1 in. = 25.4 mm) ~ 


- mathematical dependency i is supported by | earlier studies (1,13), ni it is iit 
7 j “applicable only to large-scale roughness. As relative submergence becomes larger, 
the value of A ,,/wd’ becomes so small that a power law such as Eq. 7 breaks 


- Roughness Geometry. —Ino order to quantify the the funct unction of of roughness itt 


i 
: wit 
» 
— 
| 


NOT SHOWN 


FIG. 4. —Variation of Parameter a' with Parameter Y,./w (1 in. = = 25.4 nm) 


Effective Roughness Concentration Function , — 


2 

FIG. 5.—Variation of (0 /(gdS)'/*)/(A,/ wd’) with Function of Effective Roughness 

Bey soya for 1.5 in. Bed; Corresponding Froude N Numbers, Multiplied by 100, a 
Are Given by Each Data Point (1 in. = 25.4 fom) 7 


4 
» 
| 
Zz 
| 
|| 
| | 


“size, in inches is 


Constant 


don 


*Calculated with Egs. 16 


172 
‘in Eq. 7, the ‘mie 1/(A../w wd’) was plotted against the 


a effective roughness concentration on separate graphs for each roughness 
bed. An example is shown in Fig. 5. Generally, the data points lie about one 


line for each | bed, but scatter is obvious at the lower values of b. For these © 


‘is presumably significant (7). This is supported by the variation of their auld 
_ Froude numbers which shows that for é a given value of b, the lower the rane a 
number, the higher isthe resistance. j= | 
4 In order to define the variation related to roughness geometry only, the Variation i 
oe related to Froude number must be eliminated. This is achieved by ‘choosing — 
7 only those points for which free surface drag is negligible. Since this condition 
requires that no element must be affecting the surface, it is necessary to a 
the selection procedure on the size of the biggest element on the bed. Assuming 
that the | maximum of the bed material is Si and that 
_ then appropriate points may be selected using the measurements of Flammer, cl 


Data points with F > 1. 
4 2. Data points with d/S,o = 1.5 and F > 1.3. — 
3. Data points with d/S,,.. = 2.5 and F > 1. 

4. Data points with d/S\o. = 3.0 and F > 1.0. 


This segregation reduces the number of available ite points t to eight for the 
a 4 0.5 in. material, six for the 0.75 in. material, one for the 1.5 in. material, and 


none for the 2.0 and 2.5 in. —, For these — limited data it was — 


“4 
Z =" 
= 
i 
| 


ing 


"on 
g with bed material properties Fig. 4 
Their values are given in Table 7. 

* With values available for only th three bed materials, the functi g 


ons defining 
e and m cannot be derived with great “meaning ot the 


inch 


0.75 inch 
dine 


FIG. 6.—Variation of Parameter e with Parameter Y,,/w (1 in. = 25.4 mm) = 


ers 


15 


| 


‘FIG. 1—Ve —Variation of Parameter with Parameter (1 in. = 254 4r 


on) 


~ constants is not clear and it may be that t they a are a reflection of the ; shortcomings _ 
in the analysis concerning roughness shape. However, it is also possible that, 
since e affects the € magnitude of the function, ‘it depends on the ‘number Of 


Vow 
ee 
q | aor’ 


HY12 


elements at a section, given 
Fig. 6 which shows that 


Ay 


Fig. 7 shows that m be given as a although ‘on 


relationship is less definite because of the limited range of values of the ‘oll 


oints: 
7 


Eqs. 16 and 17 were used to calculate values of e e and n m n for the 2.0 ond 25 
Free Surface Drag.—Since Reynolds number effects are assumed to be absent, 

any variations in the data still unaccounted for should be related to free surface 4 
drag only. The resistance coefficient, 4 is a “direct function of the ‘element 
drag ‘so from eqs. 7 and 15: 


d ‘in which C= the accounting for = 


From the outlined “the deag, coctticien the bed should 
vary inversely with Froude number at a constant effective roughness concentra- 
; - except at low Froude numbers. To quantify the relationship, values of a 
the function fn(1/C,,) were calculated from Eq. 18, and plotted against Froude 
_ number for different ranges of b, these being defined by a. margin of 0.05 z 
above and below values of b equal | to 0.1-0.7 (Fig. 8). The range of the data, 
and the mean value of b for each range are noted in ‘the figure. y Meee & 
There is some scatter since the results of this stage of the analysis depend : 
on the empirically derived results of the stages concerned with roughness geometry > 
and relative roughness area. However, the trends in Fig. 8 indicate clearly 
the e drag coefficient does as the relationship being: ay ( 


nes 


which i and j = constants which vary with b. war 


_ Since free surface effects diminish as relative submergence (and therefore 


should be lees at higher values of b than at lower values. Mathematically, this 
requires an inverse dependency between j and b. Confirmation is supplied by 
the data (Fig. which show 
ath 


4 __ RESISTANCE EQUATION — 
his theory is tested 
| 
— a 
7 
— 
4 
ee surface drag with 


50747 


206 
1 Froude 


= 
FIG. 8.—Variation of Function of Free Surtace Drag Coefficient with Froude Number 
and of Concentration in, = 4 mm) 


jel log (0.755/b) 


| 


+ 


at! to 


whe 


Function 


"” at a constant Froude number. Theoretically, at low relative en a 
(or values of b) drag should increase and the value of i should decrease as 
ez submergence increases, but the effect should disappear at higher r relative .. 
submergences. This is supported by the observed variation in the value of i — 
(Fig. 8) and, for the range of the data (Fig. 10): > ‘oust See 


‘Thus, data ‘generally ‘support the theory for free surface drag show 


Once rises value of 0.755, the constant j becomes zero and the fonction 
_ of the drag coefficient equals unity, corresponding to the disappearance of ie 
surface dong. For bigger values of b condition ‘should be 


1/j 


ons: 


Effective Roughness Concentration Function, 


10. —Variation of Parameter i'” with of Effective Roughness Concentra- 


@ 15 inch 
inch 


percent margin 0a 
j Line of perfect agreement » 


Measu 
FIG. 11.—Comparison of Measured Values of Resistance Function (8//)'/? v with Values 
Calculated Using Eqs. 9, 13, and 23 (1 in. = = 25.4 mm) woe Isaitercndt lgranog A 

maintained s since, physically, free surface drag continues to be absent. A value 
for b of 0.755 therefore defines the boundary between large-scale and interme- 


diate-scale roughnesses. 


- 
sha 

SS 


Rp 


in which b is defined by ‘Eq. | 13; and A wd! can be calculated for channels 
‘Using all the flume data, values of the resistance function (8/f)'” 
calculated from Eqs. 9, 13 and 23 and compared with the measured values — 
7 (Fig. 11). Since the data are those from which the equations are derived, teeement 
is naturally reasonable. It is encouraging, though, that the calculated values = 
are generally within 10% of the measured v values and do not ‘exhibit trends 
away from the line of perfect agreement. 


APPLICATION 4 


ie transport effects, ‘under the following conditions: ee 


0.41 < d/Sy< 12.10. 
3. 13<w/d< 153, 
4.0.19 < F< 1.93. es 
The values of o and R are 5 aa lower than those encountered in natural 
: rivers, while the effects of varying channel width and element shape have _ 
not been properly tested. Consequently, the equation is not fully representative _ 
of natural conditions. Despite these flaws, though, the formation of the equation | 
has allowed much of the theory for resistance of large-s scale — to — 


_ of Thompson and Campbell (20), and Kellerhals (14) and Hey’ s small- scale 
~ roughness equation (1 1), is shown in Table 8, based on the river data of Bathurst 
(1). In this limited test, Eq. 23 compares favorably, although its predictions 
_ are a little low, possibly because of the higher Reynolds numbers in the fiver — 


"compared with the steep stream flume. A more survey of large-scale 


_ Ageneral theoretical resistance equation such as Eq. 3 must inevitably | appear 
_ complex because the processes which it describes are complex. Simpler equations 


— 
1610 = HY12— 
Resi 
[8 
me 
| 
; 7 or steady, uniform flow over closely packed roughness 3 
| 
f= can be devised (1,5,14,20), but, while these may be very suitable under certain _ 7 
“a circumstances, they are less general and may need to be calibrated for individual > 
i 7 sites. They do not fully elucidate the processes of flow, nor do they show ad 7 


which ial are the most important in na given | flow, and which 1 may ) safely i 
be neglected when making simplifications. In addition, the limits of applicability * 
_ of a process-based equation can be more easily defined in advance of its use, — 
=. since the range of the processes which it describes is known. Nevertheless, 
it is -Fecognized that a general equation is not always the most appropriate 


and that Eq. 23 ‘Tequires further refinement before it can be considered to 
be of practical use, rather than as a means of testing the theory of flow resistance. - 


TABLE 8.—Comparison of Observed and Predicted Values of (8//)'/? Using D — : 


charge, | Predicted Value of (8/f) 
Using Equation of 


R/ Deo 


Cronkley B 


waves, sediment movement, and other effects alter the alter om) cane of resistance. te 


Flow | resistance of large-s scale roughness depends on the | form drag of the 
_ elements and their disposition in the channel. The relevant processes of flow 
- be described as functions of Reynolds number (affecting the element boundary 
layer), Froude number (affecting free surface drag), roughness geometry (affecting 
_ the effective roughness concentration), and channel geometry (affecting the 
relative roughness area). Full mathematical quantification of the functions is 
not yet possible, but semi-empirical analysis using extensive flume data supports 
the theory and has allowed the development of a basic theoretical resistance — 
equation. Development of a more ‘general depends on detailed ‘Studies 
‘The writers thank those Students and staff ¢ of Colorado State ‘University 


per |gence; of | Campbell Kellerhals Hey” 
2.81 2.07 2.28 | 2.23 
334] 283 | 268 2920 
CronkleyA | 1.37 | 1.02 | 1.60 | 1.52 1386 
_| 400 | 136 250 | 216] 210 =| 213 | 247 
— 165 | 3.18 | 265) 262 | 244 295 
4.00 | 1.79 | 3.33 | 286] 2.86 29 
2.17 | 4.22 | 3.38 3.23 3.880 
&g 
| 
| 
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which supported his stay at Colorado State University. & fo 
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MEASUREMENT OF IN FLOWING 


By Alan F. M. “ASCE an and nd Henry C. it) 


engineering ‘studies, in 1 both the ‘field and the laboratory. In some instances, 
- knowledge of the bubble sizes constituting the entrained air is also important. 
_ The writers have been interested in developing instrumentation to measure both 
oe ‘ air concentration (void fraction), and bubble size to assist in assessing the — 
of structural geometry on dissolved air air downstream fi from | spillways. ae — 
~. Martin (3) and Schrock (5) used gamma- Tay attenuation methods to measure 
void fraction with considerable accuracy, but those methods are not suitable 
for measuring bubble size. Much of the development of the gamma-ray method = 
has taken place in those disciplines involved with the measurement of soil moisture 
content. The work of Gardner and Calissendorff (2), for example, provides’ 
useful information for r designing a gamma-ray system. In this method, a portion 


of the gamma- ray beam is absorbed in passing through the medium. The density 


A hot-film anemometer 1 method has been used to measure both void fraction 


& bubble size. _ Delhaye (1) studied the interaction of ascending bubbles with 
a conical- shaped film probe, and Resch et al. (4) studied the two- o-phase flow 
of a hydraulic jump, also using a conical hot-film probe and a digital computer. 
In these studies, a threshold signal is used to detect the presence of a bubble. 
A disadvantage with the hot-film method, however has been its dependence 
on empirical constants determined for conditions that frequently differ from 
those to which the method is to be applied. 
major objective of this paper is to introduce and describe a s substantial 
~ refinement of the previously used hot-film method. The important “feature of 
“this refinement is use of the gamma-ray method for determining constants in 
r the more versatile hot-film procedure. A second major objective is to describe 
= movement and distribution of entrained air in a hydraulic jump which was 


affects the amount of absorption and is determined by measuring the residual 
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used to assist in the development of the hot- film a, wars ay oe 
_ Many aspects of the study have been treated seaoges but | this paper ag 


In addition, the continuity equation is used in conjunction with the void fraction -_ 
and mean velocity data to compute the spatial distribution of air flow entering © 
and leaving the hydraulic jump. The mechanics of air entrainment in a hydraulic 
jump have been en studied Thandaveswara and Rao (6) using photographic 


Apparatus For FoR HYDRAULIC, 


has A To 46 cm wide was used to generate the hydraulic jump. The test poll 

= located at the upstream end of the flume and consisted of three 73 cm 
long glass" panels on each | side. Water was supplied from a constant head d tank 

and measured with a ‘magnetic flow meter before entering the flume. The flow 
was maintained with 1% of the predetermined 0.0566 m*/s which produced 
a stable hydraulic jump ata Froude numberof6. = = 


Arrewuaron 


“Fig. 1 shows ‘the experimental system used for measuring the attenuation 
of a gamma-ray beam transmitted through an air-water mixture. A carriage 
straddling the flume was designed to move the source and detector together 


in vertical and longitudinal motion. 


The source was Cs'” which emitted gamma rays with an energy level “4 
0. 662 Mev with a half-life of 30 yr. The size of the source was 100 mc which 7 
' produced a count rate of 3,600 cpm through the 46 cm wide flume filled with — 

_ Water, and a count rate of 192,000 cpm through the flume when it was empty. 


Cs" *” was used because its gamma 11 rays are e easily absorbed by water and ore 


The source was housed in a lead shield 7.9 cm thick with a collimator hole — 
J 


0.437 cm in diameter and 9 cm long. As the size of the collimator hole is _ 


decreased, the resolution is increased, but this also reduces the count and increases 

_ Radiation was detected with a thallium-activated, sodium iodide crystal and 


photomultiplier tube. The crystal detector was mounted in a shield of lead 5. 1 

cm thick and ‘received the beam ‘through « a collimator 0.635 cm in diameter. 
The two collimators were leveled and adjusted optically for alignment. The — 
source was positioned behind the collimator to achieve maximum count. 
_ The crystal detector emitted a voltage pulse that was proportional to the _ 
energy of the radiation pulse. The discriminator \ was used to ) select only the 


scattered lower energy radiation. The discriminator had two due eliminated 
| One selected the energy level, and the other selected the energy band width 


at that energy level. A pulse must have an energy within this band width to 


_Gamma- -ray emissions are to follow a Poisson distribution which | 


| 


void fraction itself. and the time. a a high void ‘fracti 
z (10%) and five minute counts produce a . percentage standard deviation in — 
fraction of 2%. At low void fractions (0.1%), four hour counts are needed 
to reduce the percentage standard deviation in the void fraction below 40%. > 
_ Other errors can be reduced or eliminated. Due to variations in the flume 
- width and the thickness of the glass walls, counts were made with the flume — 
full of water at every lo location. _ Corrections were also made for the d ous time 
of the counting equipment and the natural radiation count. ems? ad — 
_ Another important error involves measurement near the free water surface 


igh Voltage 
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or 
where the count rate rate increases. This error is believed to be caused by y photons 
_ scattering above the surface which are not being attenuated by the water. | 
As a check on the system, the attentuation coefficient for pure water was 
4 measured and found to be within 0.6% of the ee Nhe teizeheces 


tm 
= 1(a).—Glass-Walled Flume Showing Hydraulic Jump (Carriage for Gamma-Ray 
| 
_ The rate of heat lost by a hot-film anemometer is proportional to 
a ___difference between the sensor and the fluid and the square root of the mass 
| hemometer measures a combination of 
ymperature and velocity are held constant, 


or their variation is known, the change of phase in a liquid can be detected, 
Of the two types of heated sensors, the hot-film probe is preferable for 
measuring air concentration and bubble size in an air-water mixture for two 

~ basic reasons: (1) It is stronger and more durable than the wire probe; and 

(2) it can be electrically insulated from the ambient conducting liquid (water), 

_ whereas the wire probe cannot. This prevents decomposition by electrolysis 


& and also the formation of carbonate, which can cause a shift in the signal - 


(1). A conical film ‘Probe (DISA 55R42) was used in this study. It has a nickle ' 


a quartz 
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ne.  3.—Block ‘Block Diagram of Components Use Used 


insulation coating 2 thick. The maximum-rated velocity is 25 msec for this 
probe, and the upper frequency response is approx 30 kHz. The conical probe — 
was selected because its shape provides the least distortion of available probes - 
2g bubbles moving past it, and it is able to be operated at higher velocities. — 
&s this study a DISA 55D01 anemometer was used in the constant temperature > 
a When a bubble comes in contact with the probe, there is a significant drop 
in the bridge voltage, as shown in Fig % Delhaye (1) has found that when 
the bubble first comes in contact with the probe, time is taken for the wale 
to leave the sensor. As the bubble leaves the probe, a meniscus — forms _ 


| 
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void fractions given by the two methods were equal. 


bubble has completed its passage by the probe. Both of these egoagpenrene tend 
- to shorten the length of time the probe senses it is in a bubble, 1,, compared — 
with the time the bubble would have taken to pass the same point without — 


must be added to to arrive at the duration of the bubble, _. These 


‘The duration factor was suggested by Dethaye (1) to be 3.3 msec for a conical 

probe. A shorter duration factor has been found in this study using higher 

velocities. It is shown later that the duration factor is a fonction of 
The total time that the — is in air is given git 


in which n= the total of b bubbles The fraction a is given 


4 which rT. = the total time the iia: is in water. The unknown duratio 
was calibrated using information from the gamma-ray ‘method s that 
it was ‘was found by Delhaye ( () ‘that t the water velocity next to the bubble is 
a. equal to that of the bubble. By calibrating the hot film probe for velocity, ‘ 
the velocity of the bubbles can be found. If the assumption is made that the 
_ bubbles are spherical and are not deformed appreciably by the probe (4), the 
bubble size can be “estimated knowing the bubble Ds aca and the time for 
_ The equipment used to record ‘the hot- film ‘signal is s shown schematically in 
Fig. 3. An oscilloscope was used to monitor the signal as it came from the > 
-anemometer, and a DC digital voltmeter was used to take velocity readings. : 
In order to slow the signal for recording on a strip chart recorder, a magnetic 
— recorder (Precision Instrument Model PI 6200) was used. The : signal was 
"recorded in the FM mode and played back at either a 10: 1 or 100:1 speed 
reduction to a Brush oscillograph. Using this method , 10s sec of data produced © 
- 25m of oscillograph tape. A hybrid computer was used to expedite data reduction. 
_ The EAI 860 analog computer accepted the signal directly from the magnetic i 
_ tape and controlled the sampling rate. The EAI 693 analog to digital linkage ~ 
4 converted the signal to digital form and passed it to the EAI 640 digital computer | 
for data reduction. The digital computer had 8k of storage and received programs 
through a high speed paper tape reader. Data and computational results were 
ia order to determine the optimum time length of record for the hot-film 
probe, a 30 sec record was recorded on the ‘oscillograph. The void ratio 
_ calculated by hand for various time lengths of record from 2- -20 sec. The average { : 


= void fraction, standard deviation, and percentage standard deviation of the void 


fraction \ were calculated for various time ' lengths of | record. In order to obtain 


the probe ic almoct instantaneously back in the water even before the 
q 
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7 ast standard deviation less than 10% of the n mean n void on a minimum record 
the data directly ‘from ‘the magnetic tape. The ‘program ‘required 17, min and 
printed out the void fraction, average velocity, total number of bubbles, average — 
bubble sins, ond cine 
Compurarionat Proceoures For Hot-Fiim Anemometer MetHoo 
Med of rg, ach 
Method. —Because « of limited computer si storage, each piece was 


FIG. 4.—Photograph of Scope Showing Signal from Computer Indicating Bubbles 


from the tape every 10 msec, the t time for computer to complete 
the longest loop. The magnetic tape recorder was played at a 100:1 reduction — 
in speed which, in effect, gave a data sample each 0.1 msec from the anemometer 


The method for detecting bubbles was based on the fact that there is a larger 


drop in signal voltage during the presence of a bubble than for turbulent velocity — 
If the voltage drop exceeds a predetermined value, hereafter referred 
to as the signal threshold (Fig. 2), a bubble is assumed to have been detected. 
= end of the bubble was detected by noting an an increase in voltage exceeding | 


«IG. 5.—Example of Start of False Bubble at End of True Bubble ty ia 
| 


w 


the threshold. The true bubble duration, was determined in 
_ with Eq. 1 by adding the time from the beginning of the bubble to the end 
of the bubble, ¢,, to the bubble duration factor, f,. All the bubble durations n 
_ were summed for the void fraction calculation in accordance with Eq. 2. exe 
‘The development of the program required. careful tuning so that bubble initiation — ‘ 
and termination were detected properly. To assist in this development, an 7 
. oscilloscope was used to view concurrently the hot-film signal and a signal — 
ine from the computer that indicated when the computer detected a bubble. Fig. 
4 shows the program working properly. When the broken, horizontal line at 7 
. the bottom of the photograph in Fig. 4 is in a high positi m, the computer — 


i 


FIG. 6 of | False on ral Decrease i in Flow Velocity 
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et onal 7.—Example of End of True Bubble Detection Teo Late 


_ is detecting a bubble. Fig. 5 shows a bubble having a pronounced end point — 
_ which started a false bubble at the end of the true bubble. The program would 
‘ end this false bubble at the end of the next true bubble, suggesting a bubble 
of a very large size. The computer was programmed to disregard these large 
L bubble sizes when they exceeded a predetermined value. Fig. 6 shows a false 
bubble that was , detected because of a general decrease in velocity. The false 
‘bubbles. in Figs. 5 and 6 could have been eliminated by increasing the signal — 
threshold. But when this was done, the computer would miss many of the 
smaller bubbles not shown in these figures. The center bubble in Fig. 7 is 
an example own the program did ‘not sense the end of the bubble soon 


bin got od) 
7 
1 
_ 
| | 


enough. In this ‘case the high point | in the signal did not occur at the the end of 
the bubble. The bubble in Fig. 8 is an example of where the program stopped — 

_ sensing the bubble too soon. A small peak was found on the slowly rising — 
signal before the end of the bubble was reached. 2” haere ae 
| The program was designed to eliminate the error types previously analyzed 
_ by using the average of a number of the most recent data points, instead of 
one point, to find the beginning of the bubble. The maximum number of data - 
points used, called an array, was specified at the beginning of the program. — 
__A bubble is detected when the difference between the array average and 
a the latest data point is greatcr than the specified threshold. This is the bubble 
E detection point. The end of the bubble is detected either when the difference 

between the average of the array and the signal is again less than the — 


q FIG. 8— of True cay 


generator. the as a visual check on the detection of 
bubbles by the computer, the number of errors was tabulated for a specific 
— length of run using different factors for both the threshold and the length of 
array. The minimum number of errors was found using an array length of 15 

and a threshold of 0.72 v. For these values, from a total of 677 bubbles detected P 
on the oscilloscope, 21 of the smaller bubbles, , or 3.1%, were missed by the ; 
- computer program. In addition, 15 bubbles, or 2.2%, “were not started or ended ron 

a correctly. This latter condition occurred when the bubbles were spaced very © 

_ together. The insertion of a program statement signaling the end of a 

bubble at a sharp rise in signal eliminated this problem. = vidded 
Duration Factor, Function of Velocity.—The value of the duration factor does 

7 not affect the detection of the bubble and thus was adjusted after the optimum © 
_ values of threshold and array length were found. The duration factor i is a function — 
of how the bubble passes the probe and was used to correct the length of — 
time the computer sensed the was in bubble size 


a 
| 
: The last condition was designed to end the bubble on the characteristic sharp ; 
[_ —s‘ie in signal, even if the signal never rose far enough to satisfy the first condition. — ; : 

a The beginning and end of the bubble were calculated at the intersections of 

&§ - the falling and rising slopes of the ee of the array (Fig. 

es 


and void fraction are e dependent on ‘the duration factor, which was adjusted 4 
_ to achieve a void fraction equal to that found by the gamma-ray method. — 
duration factor was found to be a function of the velocity of flow. = | : 


Data w were taken at seven | different locations in the hydraulic | jump, with = 


& plotted against average velocity i in Fig. 9. The following, exponential equations A 
the plotted data from this study: 


in which fa= factor in msec; and V= thes average flow 3. 
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FACTOR, 


DURATION 


oN 


Factor as of Velocity 


9. This value is approximately double the value predicted from Eq. 4, bietn : 
the importance of the character of the flow field on the duration factor. There = 
are two possible reasons for this difference. The first is the large difference. 
in turbulence levels between the column flow and a hydraulic jump. The second 
reason is the greater significance of buoyancy in Dalhaye’s study. The. resultant 
slip between bubble and flow would tend to distort the bubble shape. — al tci Bo 
_ Bubble Sizes. —The estimation of bubble size depends ona number of simpli- 
fications that result in a prediction method which might appear to provide only “4 
approximate values. It is possible to overcome some of these simplifications _ 
by choosing | the duration factor ‘so that the hot- ‘film method ) graces the same if 


ab 
| 
_ duration factor, determined using high-speed photography by Delhaye 


ot in establishing the duration factor. This 1 may affect to s some 

an the distribution of measured bubble sizes, even though they atin 7 
the same void fraction measured with the gamma-ray method. 28) 
‘The assumptions made are that a bubble is a sphere and the probe always | a 
intercepts a bubble along its diameter. Photographs indicate that the bubbles 
in the highly turbulent hydraulic jump are not always spherical, and the probe, 
_ of course, does not always intercept the bubble through its center. The gamma- ray 

method, however, furnishes the calibration that permits these difficulties to 
: te partially overcome. The calibrated hot-film method is able to produce average 

- bubble sizes within the accuracy of the gamma-ray method and a bubble size 

‘distribution that will be at least qualitatively 


the sampling time itself to end the bubble. This is expressed as oe ee 5 


d, the minimum bubble diameter; V= at the * 


DATA EVERY Olms 
- (USED IN THIS STUDY) 


DATA EVERY 0.0ims 


= the data sampling time. The results are 
10 as a function The effect of various sampling 
times is also shows in Fig. 10. A sampling time of 0.1 msec was used in this 


‘The hydraulic jump u used in the development of the emeniiamaai procedures 

y had a Froude number of 6 with undeveloped inflow. The total flow was 0. 0566 
m’/s which, in the 460 mm wide flume, produced a flow rate per unit width 
of 0.123 m ?/s and an upstream velocity of 3.51 | m/s. The depth just upstream 
of the jump was 35 mm and the downstream depth 1 was 280 mm. Resch, et 

| 7 al. (4) also used a jump with a Froude number of 6, but the upstream depth 

| was 12.2 mm and the upstream velocity was2.07m/s. = 

- Fig. 11 shows the distribution of void fraction determined by the gamma-ray 


| 
| | 
Minimum bubble sizes were computed using Fig. 9 and the assumption that 
would take a minimum of twice the samnlino time to detect a hubhle 
| 
FIG. 10.—Relationship Between Velocity and Minimum Bubble Size 
Ps 


method: Due to the surface e1 effe fect reviewed earlier, it was not possible to measure 
void fraction near the surface. The locations where air enters and leaves nay 
By gory jump can be identified by using the void fraction values from Fig. 

11 and and | mean velocities from Fig. 12 in the continuity equation for the air a a 


i‘ “component. Considering the finite control volume shown in Fig. 13, and assuming 
that the velocity of the air bubbles is equal to the adjacent water velocity, 


the state continuity equation for the air which expresses 


FIG. 11 —Void Fraction i in seit Jump 


Locarion, fea a 
FIG. 12.—Distribution of Dimensionioss Mean Velocity, 


‘that the net flux of air t the control surfaces -inflow) is zero, 


ol surfaces 


— in which V. = ia average velocity Nee air entering the jump over Ax; V, 


the temporal average velocity in the x-direction at a point; a = = the corre 


4 
ay 
Air Component 
nding 


average p = the. density y of air in the on and varies 

g depth; p, = the density of air at atmospheric pressure; y = the vertical 


distance from the channel invert; x = the horizontal distance referenced to | 


the upstream end of the jump; B = the slope of the water surface; and 5 
=the width ofthejump. = | 


of Air Through Vertical Cross Section 


© Computed trom veld ratio end meen velocity date 


Be 


hie Entrance Velocity, 


FIG. 15.—Dimensionless Air Entrance Velocity “4 


From Eq. 6, 


the ter temporal average velocity, of air the | at 


A 


ad 


in which A V <pady = the net t flux of air through the control surfaces 


| 

| 
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‘The integral in Eq. 8, representing the flux of air through a ney cross 
section, is plotted in Fig. 14. The slope of the curve in Fig. 14 is equal to 
a (1/cos B)(V,/V,) which, when multiplied by cos B, gives the dimensionless _ 
_ velocity V./V,, ploted in Fig. AS. The maximum value measured for B, the 
= surface slope, was Ue af which yields a minimum cosine of 0. 95. Thus, 


surface of the jump upstream from a location x is given 


in which Q, = the air flow r rate. upstream frot from the location. x; and Q, 
h fl from th eee 
the water flow rate upstream fromthe jump. 


_ The accumulated area under the carve in Fig. b poets closely the a “a 
the hydraulic jump) equal to 1.2, indicating that, on the average, air enters 
the jump upstream from this station and leaves the jump downstream from _ 
; this station. The equality of positive and negative areas shown in Fig. 15 provides a 
a check on the validity of the measurements, as there is negligible air in the 
¢ flow at the end of the test section. As a further check, air was collected 
b, volumetrically as it escaped through the water surface. These measurements, £ 
also plotted in Fig. 15, show good agreement with the values Someones from 
the void fraction and mean velocity data near the end of the jump. ~~ 
Bubble sizes for nine locations were calculated from the hot- film anemometer 
-: data. To assist in the analysis of the bubble s size ¢ distribution, a density function, 


a range of bubble size, , Ad,, divided by the range, ne... and multiplied by 
bh void fraction a. This is expressed mathematically ae ibid 
volume of air of Ad, 


A cumulative distribution function, F (d,), giving the air 


i 


| 
a the water 


1981 


that = a. al 
The cumulative Seton functions for cack of the nine locations a are bal 
_ in Fig. 16, which also lists the major characteristics of each curve. Of particular 
interest is the bubble size, (d,.), for which 50% of the air is contained in 
ye bubble sizes (@,., _ = bubble diameter for F (d,) = (a/2) 100). It varied 
ie a maximum of 6 mm at Stations | and 3 near the upstream end of the 
- jump, to a minimum of 3.8 mm at Station 9 near the downstream end of the 
= Fig. 16 also shows the importance of location on the size distribution 
2 the bubbles. The larger bubble sizes occur near the upstream end sien 
jump (Locations 1, 2, 3, and 4) where the air is entrained, and also toward | 
the the upper reaches of the flow (Location 3 a 3 compared with 2). The larger bubbles — 
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BUBBLE DIAMETER, mm 
1G. Distribution Fun netion of Bubble 
evidently are broken up by the turbulence and intense shear, and also have 
greater rise velocities than the smaller bubbles. These factors produce the | 
Variations s shown in Fig. 16. 
Suacestions For F Furtuer S 
: A more exact relationship nonin the duration factor and velocity could a 
z developed using controlled conditions designed for this purpose. The ol 
of a cylindrical probe, which has a much shorter duration factor than for the 
conical probe, should be explored. This would allow the: detection of smaller 


_ The computer program eo be modified to allow for a duration en - 


varied with instantaneous bubble velocity. The ability to record data in — 


| 
| | 
| | 
| 
| 
| 3.38] 2.61] | 0.092] 4.8 | 8.0 
5 | 4.47/ 3.48] 178 | 0.061 | 4.2] 6.0 
6 | 5.55|4.35| 76 | 0.036| 4.7| 7.0 4 
| 
pi 


restricted by a data sample interval. This would allow the computer to work . 
much faster and would allow for the use of more time- -consuming calculations, _ 
such as a velocity calibration on and a function that would correlate the 
duration factor with velocity. B 
f Procedures to determine the void fraction and bubble sizes ‘in an air-v 
mixture were developed using a hot-film probe and digital computer. A 
method for measuring void fraction was used to calibrate the duration factor 
for the hot-film method. The duration factor was found to be a function of a 
velocity and varied from about 1.5 msec for a velocity of 0.5 m/s, to about — 
_ The hot-film method is described in detail by identifying possible | sources 
_ of error and then outlining procedures for their elimination. These errors primarily 
_ involved the identification by the computer of a bubble that did not exist, and 
the failure to properly detect the end of all bubbles. Concurrent display on 
a recording oscilloscope of both the anemometer signal, and an indication of 
_ whether or not the computer was detecting air or water, was used to optimize 
The minimum size of bubble detectable by this method 
4 


a function of the size of the duration factor and, therefore, a function of velocity 
_ Measurements of void ratio and bubble size in a hydraulic jump of Froude — 
number 6 were used to both assist in the development of the _ experimental 
_ methods, and also to develop an understanding of the mo ovement | : of air in and 
out of the hydraulic jump. Use of the void fraction and mean velocity measure- — 
ments with the continuity equation showed that air was entrained upstream — 
_ from x/y, equal to 1.2. Air escaped from the jump downstream from this location. . 
_ The maximum air velocity entering the jump upstream from x/y, equal to 1. 2 = 


+ was estimated to be constant at 10. A% of the water velocity upstream from 


the jump. maximum | measured air velocity leaving the jump was equal 


_ Bubble size measurements indicated that larger bubble sizes (maximum of 
10 mm) are located near the ‘upstream end of the jump. The longevity of lenge : 
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width ° 
diameter; 
minimum bubble diameter; 
_ bubble size for which 50% of air is contained in smaller bubble 


bubble duration factor in msec; 


fa.) = density function for bubble sizes; 
- number of bubbles in test; 


“total time probe i isinair; 


- total time probe is in water; 


> time ¢ probe senses it is in bubble; 


| 


pe = average velocity of air entering jump; indi 
average velocity of air entering jump over Ax; 
= temporal average velocity in x-direction; 
= velocity immediately upstream from jump; 
x 
= vertical distance from channel invert; 
pee 


w 
: B= slope of water surface; 


&; = mass density of air in bubble; and ie 
mass density of air at atmospheric 
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INTRODUCTION 


Gachenge flow, and a simultaneous sink flow within a. discretely si stratified medium, 
are considered under conditions of radial symmetry about a vertical axis. a 
_ The motivation for the study of this problem derives from an analysis « of fo 
7 _ the operating characteristics and environmental perturbations induced by ocean 
thermal energy conversion (OTEC) plants. OTEC is a method of generating | 
a useful power using the vertical temperature gradient of the tropical ocean | as 
a source for energy availability. The upper layer of the ocean collects heat _ 
by absorbing solar radiation. ‘The underlying water is colder due to the return 
flow from polar regions which occurs within the global ocean circulation. Thus, | 
a temperature stratified — is Sores, as shown in Fig. l(a), for different 


of temperature. In total, a ‘stable. ta stratified system (Fig. 1(5)) is seems 
- consisting of a relatively thin upper layer (of the order of 50-100 m) and a 
deep lower layer, separated by a thermocline- -pycnocline with strong — 
_ An OTEC plant is basically a conventional thermodynamic engine operating 
- on the Rankine - principle. Water from the warm upper layer and the cold lower 
d layer is taken in and passes through an evaporator and a condenser, respectively. 
= heat exchangers evaporate and condense into an internal working Mid 
(e.g., ammonia) which, in turn, drives a turbine generator. The major distinction 
from conventional power plants is its very low thermal efficiency, of the order 7 nd 
of 2%, because of the si ‘small ocean temperature difference, of the order of — 
20° c; Thus, in order to , produce ‘electric power ata significant level, an OTEC 
sie must use very large water flow rates to exploit its low-grade energy resource. _ 
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For example, to produce 100 0 MW, a total flow rate of 1,000 m 
(8), which is large by any standard of engineering design. For a potential l, 000 
MW plant, the flow rate is correspondingly higher and enprenaye about half 
of the Mississippi River’s mean discharge! 
- Various OTEC design options are being pursued within the research and — 
: "development programs which ar are in progress in t the U.S. Ss. and abroad (5,11). 
One of the critical problems, in view of the large flow rates involved, is the 
potential for recirculation. Recirculation is defined as the repeated intake of 
diluted water which has passed through the evaporator or condenser. Recirculation 4 
will have a detrimental effect on plant efficiency due to a further reduction — 
of the available temperature difference. Recirculation could also lead to adverse 


(g/em?) bibs 


‘ ‘oe impact due to the local accumulation of control chemicals (for 

{ ‘biofouling and in plant operation. ATO aA 

_ The actual ocean density stratification (Fig. dD is idealized as a two-layered 

sys ‘system with densities p, and p,, and upper layer thickness H (Fig. 2(a)). This 
idealization preserves the major characteristics of the actual ambient potential 
(buoyant) energy distribution, while introducing a considerable simplification - 
into the analysis. The ambient ocean is assumed to be stagnant. This is considered: 

to represent the most restrictive condition for plant recirculation. In fact, 


7 observations on current conditions at several candidate sites for OTEC installation _ 
(e.g., Hawaii, Puerto Rico) show periods of very low velocity or of complete a 
_ The OTEC plant is schematized as a simple | cylindrical structure with balan - 
at the top (evaporator flow) and the bottom (condenser flow). Two types of 
discharge schemes are considered. te 
= Separate Discharge Scheme.—The evaporator and condenser flow are discharged | 
through separate radial slots above and below the thermocline, respectively. : 
As recirculation is much more critical in the shallow upper layer, the evaporator — 


_ discharge is located as deep as possible, ble, just above the thermocline. (In principle, 


| | 
q 
| 
| 
q 


the evaporator would also be below the ‘thermocline; however, this 
appears to be very complicated and costly ae one standpoint of internal plant 
design.) The evaporator discharge density is p, =p, + Ap in which Ap aa 
the density change corresponding tot the evaporator | flow temperature change, > 

which is small of 2° Thus, as Ap << the turbulent discharge 


bas 
Schematic density. structure b) Seporate discharge 


FIG. 2 eidiliadibidaees of Density Structure. and | OTEC Dis Discharge Schemes : 


TAA 


‘FIG. 3. 3.—Single Layer Simulation in Laboratory 


flow will ex experience strong buoyant damping v with respect to the — layer, ' 
her only weak damping with respectto the upper layer, 
Mixed Discharge Scheme.—The evaporator and condenser flow are combined | 
and discharged through a radial slot at the interface level. In practice, the flows — 


may be within the Plant or or they may mix the discharge 


| 
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¢ openings a are close t to each other—by v virtue of the Coanda effect. The adiechengs 
= 
n the current OTEC development program. 
In addition to its recirculation potential, each scheme has certain 
_ and drawbacks, as far as plant design, operation, and environmental impacts _ 
are concerned. The alternate possibility of a separate port design instead of . 


the radial set geometry has also been studied @). Its behavior has been shown 


_ Three distinct fluid zones can be expected during plant operation: (1) 1 
= buoyant jet zone, in which turbulent entrainment from both upper and lower © 
layer occurs leading to an increase in jet thickness; (2) an intermediate buoyant — 
- layer at larger distances which undergoes a transient spreading motion with © 
little additional mixing; and (3) an intake zone in which flow toward the = 
: intakes occurs. The interaction of these zones, in particular the magnitude of 
the flow rates, the degree of turbulent mixing and of buoyant damping, and 
the thickness of the upper layer, will determine whether recirculation occurs. 
This interaction is described in the following analysis. 
_ A dynamically similar problem, whose geometry is much more convenient — 
for laboratory study, is depicted in Fig. 3. Here attention is restricted to the 
_ upper layer of the stratified ocean which is simply inverted so that the discharge, 7 
Po = Pp, — Ap, occurs at the surface which forms a shear-feee boundary. An 
analytical model for this case, which closely resembles the separate Gacharge 
scheme (Fig. 2(b)), is developed first and verified with experimental data. Then 
an analytical model for the mixed discharge scheme (Fig 2c)) is presented. 


Finally, applications to prototype design are examined. eis See 


Equations.—In ‘the following, the source-t sink interaction is 
modeled by simply zones | and 3 and | the 


FIG. 4.—Structure of Turbulent Surface Jet Flow with » Potential Return Flow i in Single = 


Raoiat Surrace Jet-Sinx Fiow in Layer 
7 
| 
me 
— 
| 
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dynamics of the intermediate | buoyant layer. 1 This is justified in hindsight by 


the theoretical and experimental results which indicate a considerable horizontal _ 
extent of the jet zone, which then controls any potential recirculation.  —_— aan 
The flow in the shallow layer surrounding the source-sink is made up (Fig. 
4) of a turbulent zone with active momentum and heat diffusion, and a potential i 
return flow zone (neglecting any thin boundary layers which would develop 
under conditions of Fig. 3) which supplies ambient fluid for the intake and 
_ for the turbulent entrainment. In the radially : symmetric coordinate system fixed 
tothe still-water level, we the f ing of motion and of buoyancy 


in which r, z coordinates; u, w = mean (ime- J 
mean buoyancy = 4 = gravity; p = total density; Ap = mean 
difference; p = mean pressure; = Reynolds stresses; and j,, j, 
= turbulent buoyancy fluxes. Eqs. 1-4 assume » shania incompressible flow and 4 
nan all molecular transport processes. Further, neglecting the normal stresses 


relative to the shear stresses, and making a w 


= 


+w—]= 


Thus, Eq. 6a assumes a 4 purely hydrostatic pressure distribution. In the extreme | 
. of a recirculating flow, this omission of the dynamic pressure should remain 
_ valid as long as the horizontal extent of gc recirculating motion is larger than . 


a equation for the streamline between. rand oo: 


| 
| 
| 


in which u, = the local, "vertically uniform potential flow velocity. Thus, the _ uN 
gradient is in Ea. 5, the turbulent h, 


_ Integral equations are formed by vertically integrating over the turbulent layer 
a and observing the kinematic and dynamic conditions at the boundaries. This : 
gives, after multiplying by use of the Boussinesq and 
manipulation, the following conservation 


= 
flux: 


Buoyancy f = 


‘Buoyancy fax: der 


in which h = the ayer thickness, hy — h,; u, and g’ = maximum values 7 
_ of velocity and buoyancy, respectively; and w, = the local vertical entrainment — : 


velocity which incorporates nonturbulent fluid into the turbulent zone. Inherent 


£ in the above integration Process” is the assumption of self-similarity, i.e., an 


valid in the present problem which involves several length and dynamic scales, 
but should not be of serious consequence due to the bulk nature of the mage 


approach. Assumed self-similarity implies constant values | of the he integral « coeffi- : 
cients which are defined as 


in which = sand f and = = te velocity and | 


| 
= 
q 


an alternative conservation equation is obtained 


Up win 


“The of | Eq. 15 ‘hes in as 
_ flow force is an approximately conserved quantity, while other parameters (e.g., 
— M, u., h) can undergo jump-like transitions. This conservative aspect of M, 
is readily seen by Cone in Eq. 20 the limiting two-dimensional case, ta 
_ Thus, Eqs . 11, 13, oy 20 are the governing differential equations of the 
problem, Eq. 13 being trivial. Several algebraic equations relate the auxiliary 
variables to the three main These are Egs. from 


ag. 


Furthermore, continuity over the . entire layer implies equal return flow and © 


a: Entrainment Rate.—The only remaining unknown in the equation system is 


_ the entrainment velocity, w,. Following the usual convention in free hosp 


shear flows, it is taken as proportional t to | the jet velocity 


entrainment coefficient which describes ‘the ‘effects of radial 
geometry and of buoyant damping of turbulence. 

7 ee the nonbuoyant case and a deep stagnant ambient fluid the momentum, id 
12, indicates conservation, dM/dr = = 0. The mixing length arguments for 
free 


turbulent motions given by ‘Schlichting (13) can be readily extended to 


Using the volume flux conservation dQ/dr = a, u,., together with dM/dr 


one 0 and the appropriate definitions, gives the following relation for the nonbuoyant 


the radial indicating a Constent nt (linear) ) spreading 1 


re | 
a 
— 
ig 
i 
— 
— 
q 


} 

_ Thus, the entrainment coefficient will depend on on initial conditions r, and — 
h, and attains an asymptotic condition &, = J, k for r > », (h/r—> dh/dr). 
Available experimental data (15, 7, sad, alte extensively, 17) indeed support — 

constant spreading rate at sufficiently large distances when similarity has 
been attained. Witze and Dwyer’s data give a gradient kos = 0.106 for the - 
half-velocity width for r,/h, = 20 in which r, = the initial radius, and h, 
the slot half-height. “Somewhat larger gradients are found for r,/h, < 20, but 

- discounted here as og appear to be caused by efflux nonuniformities 


tite 


: Thus, the ‘total width is 2.27 times larger ‘than the \ velocity half-width, so that 
the experimental radial jet spreading gradient is k = 0.241. Consequently, a, 
_ = 0.054 + 0.225 h/r and is limited to the aforementioned profile. 
Turbulent entrainment is strongly suppressed if buoyancy acts on the shear 
: layer, as has been shown | by the seminal experiments of Ellison and La 
a (6). Different empirical damping factors, all dependent | on local Froude or bulk P= 
Richardson numbers, have been (19,1 12, 14). An relationship 
in which F? = u?/(g’h) is used in study. Other relationships, including 
those which give a complete “cutoff, a* = 0, at some critical Froude a 
(10), yielded practically identical results. 
Nondimensionalization.—For the present case of Reynolds number, Bous- 
sineeq- type flow, we expect any jet poopenty ® (e. .. velocity, jet depth, volume 


and seometry r, 


= Jor 


in which v2 2, J, = (assuming top-hat efflux 

profiles); and = 2nr, h, taking M, and J, as the two normalizing 

- variables because of their fundamental destabilizing and stabilizing roles, respec- _ . 

= in buoyant jet behavior, we obtain the following scales and their relationship _ 
to the non-integral initial conditions (U,,g’,Z): (1) Length scale M?/*/J!/? 


( 


(2) volume flux scale = Q, F,; 3) velocity scale 


a 
4 
q 
1 


the denotes the densimetric Froude 


tnd 


= 


Three he radial surface. jet jet problem: (1) The Froude amber 


4 AG. 5- —Radial Surface Jet as Line Source 


respectively; note that, h, J, (27 ; (2) the radius 2 
- and (3) the layer depth A. Special cases are: Radial line source, F, > 0; radial z: 


in Water 


point s source, F, @, 0; and the two-dimensional line source, F,—> 0, __ 


furthermore, deep | layer co conditions, A — in 1 which the jet dynamics 
- are unaffected by the layer depth; the incipient condition for recirculation denoted _ 

by a critical depth value H.; and shallow layer condition, with significant 

a The governing equations are put into nondimensional ‘fonts and solved with 

a fourth-order Runge-Kutta algorithm, Eqs. . 11 and 20, _ and a Newton- ap on 

_ Model Results. —The results for a radial line source in nei water (7, = 0.5, F bewir 
= 100,A = 100) are summarized in Fig. 5. This includes all integral (0 = 
_M, = M,/M,) and non-integral variables (A = h J'/?/M2/4 


F,). Three major regions can be distinguished: 


q 
= 
| 
od=f 


An initial high entrainment region in which the jet | behaves as largely , 


nonbucyant (FL= 2,M, = 1), and a significant depth increase occurs. — ied 


A gradually expanding, super-critical with damped entrain- 


be 

FIG. 6. —Radial Surface Jet: Effect of Shallowness es and Finite Sour 

the ‘entraining surface. The jet momentum function M, increases doe t to the 


contribution of the pressure term, P/r, in Eq. 20. 


-_ 3. Asubcritical expanding region in which entrainment becomes rapidly damped — 


dj The transition from super to subcritical jet flow is characterized | by rapid 


changes in @,, f, and F,. The critical Froude number if F.=V 21,/1, = - 
te 08 which is found readily es the equation for dh/dr and setting 


| 
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the denominator to zero. The solution plotted in Fig. 5 can be considered as 
“ately upstream controlled. Yet, the super-subcritical flow transition in stratified 
a flow, just as in free surface flow, should be downstream controlled as well. va 
; In the laboratory, such a control can be created - through a radial weir or — 
an unsteady layer build-up in a finite basin. In that case, and depending on 
the strength of the control, a jump transition would occur at a distance less 
than that in Fig. 5 (~40). Mathematically, this is given by choosing the alternate, 
larger root for A in Eq. 22. If the subcritical region is reached earlier in that 
_ manner, then the overall flow rate @ is limited. Thus, Fig. 5 can be considered 
as the ‘“‘maximum entraining”’ Solution to the radial buoyant jet equations. — 
a of this type are assumed i in the sequel since the downstream build-up 2 


effect was avoided in the experiment, since the OTEC field conditions probably 
lack a control due to some residual ocean currents which are always active, = 


7 


Unstable 


Radial 


an 


and since the dynamic instabilities which lead to recirculation in shallow A Zz 

5 ; conditions occur usually at much smaller distances within region a. For compari- 4 

_ - son, it may be noted that solutions for the ‘‘maximum entraining type’ ae are 

4 not possible beyond the critical section in the two-dimensional case, 7, > o, 

Be to absence of the radial expansion term. This is evident hanardipesin work — 

of Wilkinson and Wood (16) and Chu and Vanvari (4). ceaniany love ed cue 

¢ The effect of shallowness and finite discharge flows is indicated in Fig. 6, 

‘now restricting attention to 7 < 10. Thus, F, = 4 can be approximated as line 4 

a source flow, and A = 3 represents practically deep conditions. Further decrease ~ 
in depth causes a distinct region of maximum thickness, Ainax at * between 

a (2-3, and subsequent thinning of the jet until the subcritical region is: is reached © 5 


— 

| | 

| 


au! in Fig. 6a). ‘Finally, if A falls ‘below a critical value A AA. 
0.44 for 7, = 0.5), then a real solution of Eq. 22 for A cannot be found at — 

a certain position ?, which is of the order of 2-3 for A < A. and becomes 
emall as H << A... Thus, a buoyant jet flow with a given momentum function 
oan be sustained within a shallow ambient fluid layer. Physically, this condition ey 
can be —s as leading to instabilities and recirculation. It may be noted 
that A,,,, ~ 1/2 A under incipient recirculation. 

_ Finally, Fig. 7 gives a comprete radial buoyant jet stability diagram asa 


a 


x 
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wre. 
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Exp. 33 dota range) 


FIG. 8.—Radial Surface Jet: 

_Non-Reci 


Comparison of ‘Theory and Experimenta Data = 


function of radius 7, as obtained by repeated ales application. For 7, < 0. s 
the behavior seems independent of 7, giving the radial point source regime = 


A. = = 0.42; while the range 0.1 < 7, < 10 is characteristic of radial line sources 
e- some dependence on?,.A different behavior (A. ~ F>'/? due to differences 


in the line source quantities) is ir ag to result for large 7, > 10 approaching — 


| 


| 
| 
| 
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device could be installed either in the basin c center or as 
 half-model at the basin boundary which then can form a symmetry line for 
¢ flow field. A window along the basin boundary allowed ine ied 
observation. The jet was formed by discharging heated water through a radial iy 


slot at the water surface. An equal amount of water was withdrawn from near 

the bottom. A series of 28 experiments was performed with stagnant layer » 

conditions, covering the following parameter range: H = 0.24-2.5; 7, = 0.08-0.83; mi 
_F, = 2.0-30.5. The discharge Reynolds number, based on slot height and velocity, 
varied from 3,000-12,000. Temperature measurements were taken with 100 
thermistor probes (1.0 sec time constant), mounted along five radial sections 


on a movable frame so that they could be positioned at any elevation. The 


cel 


© Exp. 


56 


oat a 


FIG. 9.—Radial Jet: n of Data for 
readings of these probes were automatically recorded and radially averaged, 
3 thereby producing the mean temperature field and, thus, buoyancy field field g’(r,z). a 

- Other measurements included: vertical temperature profiles via a traversing fast 

thermistor probe, velocity profile determinations from dye streaks released from | 
_ dye coated vertical threads, flow visualization through dye release into the 

; discharge water, and direct recirculation measurements by monitoring the 
temperature or fluorescent dye concentration in the intake flow. Details of the 
—— set-up and procedure are given in a separate report (8). tek 

Two experimental error sources are of concern here, namely radial nonunifo or- 

; penn and unsteady temperature build-up. The radially expanding flow within 
the model tended to instabilities and thus unevenness in the exit flow distribution. . q 
This effect was minimized through use of internal baffles and flow dissipators. J 

Repeated checks of the surface temperature distribution showed satisfactory 
radial symmetry, although minor deficiencies (of the order of +20%) persisted. 
— The Sate basin size limited the steady-state experimental period to about 20-30 


Experimental Study.—Experiments were conducted in a 12.2 m x 18.3 m 4 
bes 7 0.35 m shallow basin at the R. M. Parsons Laboratory for Water Resources _ 
| 
| 


oe. . After that, the mixed water buildup in the basin “‘far-field” started 
to move in toward the near-field with hydraulic jump-like conditions and ultimate 
near-field flooding. The data taking was limited to the steady-state period. —__ 
All of the experiments were within the shallowness affected range, A “a 
3. Fig. 8 summarizes experimental data and model predictions for several stable — 
(i.e., non-recirculating) experiments. Temperature half-depth, hos ; (= 0.63 A for a 
the chosen profile), normalized centerline buoyancy, g’/g’,, and, in some _ 
instances, total layer depth f, as obtained from local velocity profiles, is given. 
The data range obtained in the detailed investigation of Experiment 33 is typical 
for all runs and a probable measure of residual radial nonuniformity. Except — 
= a slight over prediction of the initial non-buoyant jet growth, f, which 
is due to the direct adoption of Witze and Dwyer’s ond rate, the agreement - 

Fig. 9 shows two experiments in which A was obeerved to occur. 

This is coincident with an instability occuring in the mathematical model, = 

A< A.. The model predictions are terminated at the instability. Experiment — 

6 (A = 0.35) is a marginally unstable case. A recirculating motion could readily in 
be seen through the lateral observation window. However, no recirculation was 
- measured in the intake flow. This is due to the fact that the inward return 
- flow does not become uniformly mixed, but has in itself a layered structure: 

_ the upper portion consists of the recirculating flow which gets re- ~entrained — 

' into the buoyant jet, while the lower portion consists of ambient flow which 

enters: the intake. Experiment 8 (A= 0. 24) is a more extreme case of recirculation. — > 


; rapid i increase of the experimental jet depth, wt At the same time, the entire 


water depth, A, was observed to be occupied by a recirculating eddy and a 
direct temperature rise was in the intake flow. 


InterFace Jet-Sink Frow in Laverep Ocean 


‘4 “consisting of a finite upper sit and an infinite lower layer. Thus, return velocities _ 
and egret deviations in the lower layer are negligible. ‘The 1 radial momentum 


before, and the same global dynamical interaction for the surface jet is expected 
_ to govern mixing and stability of the interface jet. Assuming vertical symmetry 
for the turbulent zone, i.e., A, = —h, h, = +h, the integral equations, in 
conservation form for the upper jet portion, are the volume flux eqeatien, 


of the flux equation, Eq. ‘The | global th and 
_M, are defined over the half depth A, and all auxiliary variables, including 
entrainment rates, are defined as before with exception of h which is now 


the solution of — 


i 
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Using M, and g the n normalizing the 
length scale = volume flux scale M*/°/g’'”? = 


and velocity scale M 4 Fi! in which the interface jet Froude 


“layer 


ent 


10.—Structure of Turbu u 
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- CO a FIG. 11.—Radial Interface Jet as Line Source with Deep Upper Layer | tat _ 


Eq. 31. The nondimensionalization of Eq. 35 i Sis 


Model Results. —Fig. ll gives the results for a radial line source at the 
_ with deep upper layer (7, = 0.5,F, = 100, = 100). The same general features: 
_ for the surface jet (Fig. 5) occur. The major difference is a distinct maximum 4 


value of the jet depth A at 7 ~ 2.5, with a subsequent decrease in depth. 
This is related t to the considerably larger buoyant damping effect. For the same 


F~ 18) 
ing 
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point Radial line source» 


12.—Radial Interface Jet Stability Diagram: Condit tion o Recirculation 


is possible and the instability condition is reached. Fig. 12 is the interface jet 
Stability diagram, again indicating the point source range (7, < 0.1, A. = 0.50), a ; 
<j and the line source range | for mere 7,. Beyond 7, = 10 the fully two- dimensional 
line source regime with A, ~ is approached. 
The theory is now applied to typical OTEC applications in tropical ocean 7 
latitudes. First, the mixing characteristics of a 100 MW OTEC plant are examined. 7 
_ survey of ocean data and proposed designs (8) give flow rates Q, = 500 
m’/s for both evaporator and condenser flows, a conservatively low discharge pe: ; 
velocity U, = 2 m/s, a plant _ = 25 m, and a yearly average value sind 


forthe mixed layerdepthH=70m. 


Separate Discharge Scheme .—The flow is into the upper 


layer and a ‘ 
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to om Its discharge buoyancy relative to the lower layer 

_ is considerably higher (factor of 5-10) so that mixing with lower layer water 

is neglected. The “inverted surface jet’”’ variables are, therefore, H=0. oo, 

‘ fF, =0.22, 2, and iF, = =1. 12 indicating line source behavior with significant shallowness — 


a Mixed Discharge Scheme.—The combined flow is emitted at approximately — 
_ the thermocline, which is characterized by a total buoyancy difference 2 g’ t 
= 0.02 m/ s* (corresponding to the density change over the | upper 150 m Nad 
the ocean (Fig. 1(5)). Thus, the interface jet parameters are H = 1. me Se 
0.54, and F, = 5.0, giving a relatively deep line source regime. The siniad ‘ 
predictions for jet depth at an undistorted scale, and total induced radial flow _ 
in the upper layer, are summarized in ‘Fig. 13° for a near-field region of 500 
m surrounding the OTEC plant. Thus, the separate discharge scheme tends 
a4 to induce considerably more vertical perturbations and mixing within the ba 7 
ocean structure than the more damped mixed scheme. 


_ Secondly, the maximum OTEC plant size can be calculated under the — 


evaporator discharge depth 


et boundaries for 


wi m Mixed Layer Depth 


of recirculation. In view of this of 

_ OTEC power delivery to the available temperature difference, this is, indeed, 

a reasonable objective. Also, since recirculation is tantamount to a full agitation 
of the entire upper layer, it is further undesirable as it destroys any small 
temperature gradient within the upper layer | and precludes utilization of the 
warmest near-surface waters. Despite strong data campaigns in recent years, 
the environmental data base for OTEC is still very limited, particularly for : 
extreme conditions. Yet it appears that, for many OTEC localities under current P 
consideration, the minimum mixed layer depth is of the order of 50 m. It seems 
prudent to reduce this by: another 20% to allow for short- term weather related 


Discharge Scheme.—The point source recirculation ‘condition, = 
42, can abe expressed using the earlier definitions (M, = Q, U.): 
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With the design from above, the evaporator | = 250 m m/s 
which corresponds to a 50 MW plant size. 

Mixed Discharge Scheme.—The recirculation , A, = 0. 50, translates 


= (4) — ( ) 


This is evaluated as 0, = 2,560 m */s which stands for one half of the 
‘mixed discharge flow so “that the corresponding plant size is about 500 MW. 7 
Thus, the separate discharge of the evaporator flow in the upper layer is extremely | 
prone to recirculation. Mixed discharge schemes seem preferable. In any case, . 
the results suggest that the external fluid mechanics of OTEC plants are A 
‘ definite limiting factor on plant size and seem to preclude very large plants — 
* the 1,000 MW range) which have been suggested i in the literature. ia cel 


In summary, the present analysis appears to give reliable first-order predictions 


on the interaction of a large momentum source and a sink flow within a layered, ol 


- Stratified fluid system. The model appears well verified under experimental 
5 conditions which are dynamically similar to those of the tropical ocean. The > 
_ model predicts the amount of mixing under normal operating conditions and, 
7 ‘if the system is stressed to the extreme when emitting large momentum sources, 


it indicates a complete breakdown of the orderly layered structure of discharge _ 
withdrawal zones, 
_ The present study focused on the schematic essentials of the source- sink 
interaction. Numerous other complexities exist in the actual OTEC problem — 
_ and are the subject of more detailed experimental and theoretical work, e. “2 ‘ 
[ the effect of continuous stratification (2), the effect of ambient cross flow (8), 
- design details in form of discrete, rather than radial, discharge ports (2,8), and — 
transition to the far-field (9). 
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The following symbols are used in this paper 
interface jet Froude : number; er; 
local Froude number; 


surface jet Froude number; Reyicn, which 
gravity; 


= initial (surface jet) or total (interface jet) Saiasildan: ~ 
= critical layer depth for incipient recirculation; 
integral constants; 
= total radial buoyancy flux; 
= turbulent buoyancy flux; 


| 
; 
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a radial and vertical distance, respectively; 


= sink flow return velocity; 
non-buoyant radial entrainment coefficient; wil 


buoyant radial entrainment 


‘a 


Ss 


tl 


1,2 = upper, lower layer; 
“centerline (maximum) ré 


normalized ire jet Property. 
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Verticat Rounp BUOYANT JET IN SuALLow WATER 


By Joseph H. W. Lee,’ A. M. ASCE and Gerhard Jirka,? 


a A simple we way to dispose of the large « quantities ‘of waste heat, resulting from 


- steam electric power generation, is to discharge the heated condenser water 
_ through a submerged round outfall located at the bottom of the receiving water. 
An understanding of the induced excess temperature distribution, in relationship 
to the hydrodynamic characteristics of the buoyant discharge, is essential for 
: a sound design which utilizes to its fullest the cooling capacity of the ambient — 


7 waterbody. The information also provides a basis for an assessment of the 
impact of the thermaleffluent. = 

In recent years, high velocity injection as 
; have increasingly been used as effective mixing devices to dissipate the heat 


-‘Tapidly | within a small area. . Many once- through condenser cooling systems — 


discharging either vertically or horizontally. The San Onofre Macher Plant and 
- Redondo Beach Fossil fueled plant on the Pacific Coast, and the Littlebrook _ 
and Sizewell power stations on the Thames estuary, United Kingdom, are two 
cases in point. Frequently, these discharges are located in waters only a few a 
port diameters deep. Significant surface and bottom interaction may then dominate al 
mixing characteristics of such depth-restricted thermal discharges. 
_ This paper reports an analytical and experimental investigation of the stability. % 
and mixing characteristics of a round turbulent buoyant jet discharging vertically 
_ into an otherwise stagnant uniform waterbody of large horizontal extent. The © 
- discharge configuration is termed stable if the discharged fluid is not re-entrained _ 
‘into the > buoyant jet (Fig. 1@)). In such a case, four flow regimes of distinct — 
which the discharge mixes with the ambient fluid by ‘turbulent entrainment as 
it rises to the surface; (2) a Surface Impingement Region, within which the — 
- buoyant mixed flow is turned into a radially outward moving supercritical flow; © 
_ (3) Radial Internal Hydraulic Jump Region, transforming the flow into a subcritical — 
 - ‘Asst. Prof., Dept. of Civ. Engrg., Univ. of Delaware, Newark, Del. 19711. 
Associate Prof., ‘School of Civ. and Environmental Engrg., Cornell Univ., Ithaca, — 


- Note.—Discussion open until May ‘1, 1982. To extend the closing date one month, 
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ASCE. Manuscript was submitted for review for possible publication on October 17, 

1980. This paper is part of the Journal of the Hydraulics Division, Proceedings of the 
suelo Society of Civil Engineers, @ASCE, Vol. 107, No. HY12, December, (1981. 
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_ region downstream with an an abrupt change | of up upper layer thickness and loss" J 
of energy; (4) Stratified Counterflow Region, with negligible entrainment across 
_ the interface. In general, stability can be maintained for deep water high buoyancy — 
pseeonanes that are characteristic of sewage applications. For these cases, the — 
_ minimum dilution at the surface in the vicinity of the discharge can be estimated 
using various buoyant jet models that have tacitly assumed the existence of 
‘ On the other hand, under high momentum, low buoyancy, and shallow water 
conditions, the energy input by the buoyant jet into the receiving water may 
be so strong that the buoyancy contained in the discharge is not sufficient | 
to establish a stable stratified flow. Recirculation cells are set up near the jet — 
efflux, leading to heat re-entrainment into the discharge and thus an impairment 


in mixing saith Such a discharge configuration is termed d unstable Fig. 


FIG. 1.—Vertical Round Buoyant Discharge in cnatow Water: (a) Stable — 


jet in agua water has been studied in considerable detail by Jirka and Harleman 

_ (7). The same line of approach has been adopted in the present work: ~The — 

: salient features of the four regions are separately analyzed and the solutions 


the theoretical predictions. Practica! implications of the results are discussed. " . 


Consider a vertical axisymmetric buo: 


1652 
= 
| 
| 
| series of laboratory experiments carried out with a half jet are compared with — - 
| 
: a ant jet of fluid (discharge velocity u,), — 


and the field T(z, r) generated by the buoyant: discharge 
_ (temperature T,, density p,) in a cylindrical coordinate system. a. 
In the following sections, a steady-state formulation of the problem, which 
incorporates the dynamic interaction of three near-field regions—the buoyant 
4 jet region, surface impingement region, and the radial internal hydraulic jump 
region—will be presented. The buoyant jet region is analyzed with the assumption | 
of a constant jet spreading angle. the impingement and 


REGION OF FLOW 
es 


he att 

ior pare 


jump regions are determined on the basis of. control 
A coupling of the solutions of the three regions yields ; a prediction of the near-field = 
stability, as well as the bulk mixing characteristics of the buoyant discharge. _ ; 
a Buoyant Jet Region.—The mechanics of a pure axisymmetric turbulent free 
det or plume, in an infinite field, has been extensively studied (2,17,19). Thus, 
the detailed physical description of jet or plume behavior will not be belabored 
upon here. Suffice it to say, in the general case of a buoyant jet—a source 
- with finite initial volume. momentum, and buoyancy fluxes—the governing 
dimensionless pasamater is the Froude — 


issuing from a source of fi id 
il 


which Ap, =p, — p,/p. = the discharge relative density difference; 
i and g = the acceleration due to gravity. The extreme values of F, (the ratio _ 
4 of inertial to buoyancy force per unit mass) define the two limiting cases of | 
a ‘pure source of momentum, F, => @, and a pure source of buoyancy, F, a : 
'% Fig. 2 shows a schematization of the flow structure in the plane of symmetry 
of the buoyant jet. This consists of a region of flow establishment near the 
source, a = z,, in which the turbulent shear layer has not penetrated to the 
- jet axis, , and a region of established flow, z> z,, in which velocity and excess 
temperature (density deficiency) profiles are similar at all levels. The following 
Pers layer type approximations customary with buoyant jet analysis can > 
be introduced into an integral model, which accounts for the effect of buoyancy ~ 
in the zone of flow establishment: (1) The piezometric head is constant across 
the jet, and the transverse mixing scale is much less than the longitudinal scale, : 
<< i, which is the nominal jet width. Turbulent transfers in the 
Z z direction are neglected; (2) the Boussinesq assumption is such that small density 
_ differences can be neglected except when multiplied by the gravity term. For 
small density differences induced by temperature differences, a linearized 
; equation: of state relates excess temperatures AT = T — T, to density deficiency 
Ap = .- Pp, Ap = BAT, in which B = the coefficient of thermal expansion 
at T,; (3) the vertical velocity, u (2, r), and density deficiency profiles, Ap(z, ry, 
are assumed to be similar. In the region of flow establishment, a linear decrease 
of width of the potential core, b’, is ‘The velocity and density 


u AZ r)=u, and r)= dp. for 


in which \ measures the difference in diffusivity of momentum and heat. 7 rs 
In the of established flow, the profiles : are assumed to be fully ¢ Gaussian: 
in which the subscript ¢ center ter line these assumptions, 


be integrated across the buoyant jet from 1 r= 0 tor=@ to yield a set of 
governing equations in terms of the volume flux | Q, momentum flux M, and 


Continuity; — =— 


_ Momentum: —=— | utardr= \ rdr=f 


fj 


HALLOW WATER JET 
| Apu,2ardr= 
bail ow ot al gaits, avewied nad 
in which Q, = = the rate of entrainment of ambient fluid into the buoyant jet. 
quantities, Q, M, J, can be related to the centerline 


In this type of integral formulation, the closure problem in turbulent diffusion 
7 im to relating the entrainment flux to the local characteristic quantities. 
In many buoyant jet applications, (9,6), this problem is circumvented by an 
entrainment hypothesis first proposed by Morton, Taylor, and Turner (17) to 
describe the broad mechanics of the flow rather than the details of the actual 
mixing mechanism. The entrainment flux, Q., is assumed to be proportional 
to the local centerline velocity and width as Q0.= 2nabu, in which a = an > 
experimentally determined entrainment coefficient. Invoking this assumption 
in Eqs. 4-6 yields a set of three equations in u, , Ap., and ~ 
S. _ Based on purely dimensional grounds, and confirmed by experiments (6, 12), 
the entrainment coefficient, however, can be oom ¢ to be dependent on the 
densimetric Froude number: 


and thus is variable in the _ case of a — jet. ‘The exact whee 
a ‘depends on the shape of the profile assumed. For Gaussian profiles it 
is equal to a, ,= 0.057 for a pure momentum jet, and a, = 0.085 for a pure 
plume. The considerable variability of a necessitates a numerical solution of 
= 4-6, if some form of dependence of the entrainment coefficient on F, 
As an alternative that is deemed more _ appropriate for the present work, “4 
the entrainment hypothesis can be replaced by an assumption of a constant 
_- spreading angle, first proposed by Abraham (1). hot diffusion layer, b, 


is assumed to spread linearly; ad} mdi 

db 


The cor continuity equation then can be bypassed and an analytical solution can — 

be obtained by solving Eqs. 5 and 6. The coefficient «, of course, has to be 
experimentally determined. Jirka et al. (6) have shown that, in order for the . 
‘ two approaches to be consistent, certain relations have to hold between a and o 


7 €. For the pure a momentum jet and plume, it can be shown that that «, 


q 
| 
| 
| 
i 
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0. ond: = 6/5 a, at 106 independent of the profile assumptions. Thus, 
‘wheress a varies considerably as a function of F,, there is a difference of — 
less than 10% between the jet spreading angle in the two limiting cases. Throughout ia 
- the subsequent analysis, « and \ are taken to be 0.109 and Ld 14, , Tespectively. ; 
represent intermediate values in the jet-plume range. 
Following Abraham (1), the solution of the characteristic quantities in the | 


region of established flow can be found: 
2 


D 


in which subscript e= = the ontana at the ‘end of the he region of Sow establishment. a 
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FIG. 3.—Length of Region of Flow Establishment as a asa Function of Source 


a The momentum flux within the region of flow establishment i is given aad integrating 


wae 


ae 


= 


| 
| 
‘ 
{ 
4 
Defining Z = z/D, Eq. 12 can be eval = 7 


length of establishment, Z,, then can be evaluated ated by invoking Ap/ Ap, 
= 2, in Eq. 11 and combining with h Eq. 13: od) ta 

The of Z is down in Fig. 3. Z, increases rapidly 
i from zero for F, = 0.0, to an asymptotic value of 5.74 for F, beyond 25.0; 
this is somewhat smaller than the value of Z, = 6.2 determined by Albertson 
et al. (2) for the three dimensional momentum jet, indicating that heat diffuses — 
faster | than momentum. (Eq. 14isa slight it modification o of Abraham m’s ’s solution 
(1) for rZ .» the effect of buoyancy being properly a accounted for in the evaluation 
of M.). In shallow water, the length of flow establishment can constitute a 
_ significant portion of the water depth, and its inclusion in the analysis is important. 
Given Z,, the characteristic quantities Ap., u. can be obtained from Eqs. 10 


and for z > z,. For z <z,, te 


= ‘Surface Impingement Region 


When the buoyant jet impinges on ‘the free the pressure 

- causes a surface hump, which gives rise toa radial pressure gradient and hoviontal 
_ spreading of the mixed discharge (Fig. 4). Intense turbulent mixing occurs in 

- - this region. Instead of attempting to describe the complicated flow details, a 
_ ‘simplified control volume analysis is given in this section to relate the properties 

, of the vertical buoyant jet at the entrance of the surface impingement region, 
the characteristic parameters of the horizontal axisymmetric surface 


@ 


"jet at the exit of this 1 region. For brevity, ‘the flow is assumed to to be fully 
established at the entrance section i, with center line velocity u,, density deficiency 
‘Ap, and jet width b,. The case when the entrance flow is not fully established 
can be treated in an entirely analogous manner. The following assumptions - 

; are made: (1) The outer radius of the region, r = r,, corresponds to the location 
in the buoyant jet where the vertical _ Velocity ‘is 5% of the center 
line value, ice., r= V3 b,; (2) at the exit section, the vertical profiles of 7 
the radial velocity and Geasity deficiency are assumed to be —_ we 


z in which h, = characteristic length scale of the 1 upper layer. An appropriate 


= of the upper layer thickness A, can then be defined as the point in 
q Neglecting turbulent entrainment into the the governing 
_eeations f for the control 


Conservation: .. 


In In_the mechanical — y conservation 


Ap,z refer, r 


Kei is a te loss coefficient, accounting for the turbulent energy dissipation 
in the transition. K, = 0.2 is adopted as appropriate for a 90° bend and a ; 
wide range of (5). The results are not sensitive to K, 
The preceding analysis is applicable only if z,> z, . When entrance 
_ flow into the surface impingement region is not fully developed, and thus the 
_velocity and density deficiency profiles are not Gaussian. This has to be accounted 
fori in the analysis because: (1) Very often practical submerged discharge designs snail 
‘fall within the small ai. range; and od z, is also a function of F, (Eq. 14). 
AY based on the same principles just ; 
“outlined, can be y , 
Coupling Buoyant Jet Surface Impingement Region.— n.—Noting that the “2 
A ey condition relating entrance section i and the exit section J is z, = 2 
3 h,, and b, = €z,, Eqs. 17-19 can be rearranged to yield Benes 


| 

i 

18) 

19) 

§ nd 
| especiively, to the proiile-averaged velocity and 

be shown that 

1 

ry 
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s the local densimetric Fr number of the jet; this can 


related to the source conditions by ¥ wht 


or Surface Impingement Region Upper Layer matin? 


4 FIG. _ Solution for Surface Impingement Region, Densimetric Froude Number ; 
a. 


= the analytical solution 7 
@ 
7 | 


(Eqs. 10 and for any source conditions The atites for the 

7 dimensionless upper layer thickness h, _/H, and the densimetric Froude number 

defined by F, = u /Vghp,/p,h,are shown in Figs. 5(a) and 5(b), respectively. 

7 Upper Lew Thickness. —In the deep water range, i.e., when z, wh J/A 

is quite insensitive to F, or H/D, and is equal to about 0.08 or 1/12. Sbiaseee, 

as the shallow water range is approached, h,/H increases to about 0.2-0.4 
in the low H/D range. It is interesting to note that for H/D < 4, h,/D ~ . 
0.4, indicating that D replaces H as a characteristic length scale in the Shallow 

_ water range. Given the ; Solution of h,/H, all the characteristic properties at 

the entrance section, u,, » Op, F,, ca can be evaluated at : Z,. The near field bulk 
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G. 6.—Radial Internal Jump 

dilution i is s then en related ‘to the buoyant center line concentration as 

+? Ap, 

Ap, 2d? Ap, 


Furthermore, the upper densimetic number can be obtained as 


« 


_ Upper Layer Densimetric Froude Number, F,.—The pure plume is cheracterined 

a constant asymptotic local densimetric Froude number | FF, = 3A7/2€ Aw 

4.45. Regardless of the discharge densimetric Froude number, F ,, F ,approaches 

_ for su sufficiently large values of H/D showing that all buoyant jets tend 
to exhibit plume behavior far away from the source. Within the range of F,, 
im in Fig. 5(b), this trend is more apparent for small F,. When F, > F,, 7 
+t | decreases as the deep water range is approached. As a contrast, if F, < => 
"ae then as the submergence H/Di increases, F, will first increase and then ‘de- 
crease, indicating fluid To when F, is ‘below the asymptotic value F, 


As the supercritical upper basis flow moves radially outwards from the surface - 


| 
— 
| 
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SHALLOW WATER JET 
impingement region, it entrains fluid from. the lower - layer spreads like 


an axisymmetric surface jet. At sc some me radial distance, 1 r=P,, the high momentum, 


flow, with a thickening of the upper layer depth and accompanied by a loss 
of specific energy (Fig. 6). The characteristic properties of the ‘upper layer 
at the toe of the jump, r = "sad can be coupled to those at the « exit of the 


‘surface impingement region. A bulk momentum analysis yields the conjugate 
jump heights of the respective layers. iginds ne bas 
_ Coupling Surface Impingement and Internal Jump Region.—The governing 
equations for the upper layer flow prior to the jump are (assuming half-Gaussian © 


¢ velocity and and density deficiency profiles as before (3,8); wollads ei is 


a 

d 

in Q, = the volume, momentum, and density pope fluxes 
of the radial surface jet; and a = an entrainment coefficient; Q= 
de. M = dz; and J = 
"Substituting the assumed profiles into Eq. 26, and imposing the known conditions 
at the end of the surface impingement region, r =r,= 3 b, (~0.2 A), 


which = dn. The entrainment coefficient for the Tadial 


Bh ds jet: a = 0. 057 for Gaussian profiles. The ‘buoyancy induced radial 
3 pressure gradient that has been ignored in Eq. 26 has been shown to introduce 
a anegligibleerrorforr/H= 
sy Defining the upper layer thickness at r=r,ash, 5 3 h,, the average 
upper layer velocity u,, and density deficiency ‘Ap, at the jump toe, can be 
evaluated as u, = Q,/20 ryh,; ;Ap, =Ap,/V2 which subscript J = properties 
Internal Ju Jump Analysis. —The flow in the internal jump Tegion is treated as 
a two-layer system | m with a distinct interface, , and with constant densities p,, - 
Pa — Ap. The upper and lower layer depth at r r= r,are ae h, respectively, 
and their conjugate jump heights at the end of the jump, r = r,, areh|,hj. 

_ Whereas radial internal jumps heve hitherto received little attention, both 
the two-dimensional internal jump, and the free surface radial hydraulic jump, | 
4 have been studied extensively (7, 15, 16, 22 An momen- 

= by Yih ont Guha (23), and the solution to > angen system 
of equations sought by various investigators. Using a Boussinesq approximation, _ 
the two- dimensional | internal jump \ was also ‘Studied d by J Jirka ¢ and Harleman 2) e 


4 
a 
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in the context of the stability of plone vertical buoyant jets in cusfined depth. 
On the other hand, the radial free surface hydraulic jump has been studied 
7 experimentally and theoretically by Sadler and Higgins (20) in connection with 
s design of channel transitions for free surface flow. In this case, since the 
discharge per unit depth varies with radial position, both the location of the 
jump toe and the length of the jump, r, — r,, may be important determinants | 
of the resulting conjugate jump heights. By imposing a downstream control, — 
- assuming an abrupt vertical jump of zero length, Mehrota also presented a 
a solution of the inviscid circular free surface hydraulic jump (16). pape) 
The present problem, i.e., the internal jump created by a vertical round buoyant 
“jet i in shallow water of large horizontal extent, differs from the previous works 
_in two major respects. Firstly, there is no downstream control in the near field 
of the vertical buoyant discharge. The present investigation reveals 
the jump occurs at at approximately » = 6. 6H 


the conjugate jump height, 
of the jump toe, r, ~ H, it can a shown Ql) that the length of the jump 
cannot be neglected in the momentum analysis. Whereas there is no basic _ 
experimental data on the jump length of radial internal jumps, the little data 
available for radial free surface jumps al allow a crude approximation it in relating 


interfacial and bottom shear stresses, assuming 


pressure and uniform velocity distributions, and an average head of 1/2 (h, 
+h‘) over the jump length, a force balance of the two layers gives the governing 
for the jump heights h{ , of the respective layers 
2 2 


in which F*, Fr = = the free surface Froude numbers defined 
in which Q,, = discharges in the upper and lower respectively: 


= aa aes = Pa — Ap. The layer depths before and after the jump represent two 
ible dynamic states for the same given flow force. By equating alternate — 
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expressions of hy (29 and 30, and introd 
the Boussinesq approximation, Ap/ p << l(a condition satisfied a 
submerged buoyant discharge applications), to the resulting expressions, the 
governing equations for the internal jump become 
1+—}l1-— 
1+ 


b 


1+ 


(33) 
a in which Fi = Ft?/Ap/p; F? 2 = F? For typical values of F,, F,, and 
=< Ap/p, the solution to the approximate Eqs. 32 and 33 differ from that of the 


exact governing Eqs. 29 and 30 by, at most, 1%. It can be shown that - 
yy 32 and 33 imply a constant total water depth h, +h,=h{ + Ay which results 


3 in the equation governing the upper layer change: 


a 


1+ 


Defining the critical state as that which is conjugate to itself, h’ <- h,, hs 

5) 
The flow is supercritical upstream of the jump (compare Fig. 5), Fj + F3 


7g 


general, the governing for the conjugate jump height, 
two solutions that lie within the physically feasible range 0 < hi < H. One > 
of these can be discarded by energy considerations. However, for high densimetric 
Froude numbers of the upper layer F, (or large values of h,/h,), Eq. 34 does 
-— yield a solution. The nonexistence of a solution can be interpreted Physically 
asa hydrodynamically unstable situation. Under these h high momentum conditions, 


the buoyancy contained in the upper layer flow is not sufficient to “establish 2 


s and vemmens impingement — (Eq. 21) can be coupled to the internal an 
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SUBMERGED 
THEORETICAL 
PREDICTION 


oo UNSTASLE NEAR FIELD 


—Stability Criterion « of Vertical Round Buoyant Discharge in Shallow Water | 


_ jump region via the surface og solution (Eq. 27) to determine the flow conditions { 
prior to the jump at r = LA predictive criterion of the near field = 
is then provided by the ‘egietion for upper layer change, Eq. 34. In terms 
of the discharge parameters, this stability criterion—the line delimiting the stable 
Fig. 7 shows that for a given H/D, the domain is papproached 
by increasing the source densimetric Foude number F,. The departure of the 


stability criterion from Eq. 38 for low H/D indicates initial size effects. It 


conjugate downstream Conditions, leadin oa ak 

flow and heat re-entrainment into the bu 

| 
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WATER JET 

is is interesting to note . that the pure pone (low F,, 


determined solely by the dynamics va the jet, surfa 


and internal j jump regions. * a 
A total of 33 experiments were carried out with a vertical semicircular buoyant | : 


wil jet (i. €., a jet in lieu of a full round Bom ina x Bx q 


Fi.OW METERS | 
@ SYSTEM OF 


SYMMETRY 


WOODEN FRAME 


SLOTTED WEIR ps SUPPORTED ON 


MOUNTED AT 
SAME LEVEL 


ADJUSTABLE _ 
WOODEN FRAME 


1-43 
sucks 


x 5.5 x 0.3. } m) hydraulic scale model basin which had a false floor made > 
of styrofoam insulation material and concrete blocks (Fig. 8). A 16 x 34 ft 
ai (4.9 x 10.4 m) area was formed by constructing a partition which formed a 
a plane of symmetry for the axisymmetric jet along the whole length of the basin. 
a _ The center portion of the partition was made of 0.5 in. (13 mm) thick Plexiglas, 
partial flow visualization possible at an oblique angle. 
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(14.7) 


(15.8) 
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The flow injection device half jet was a ‘Plexiglas box 
ith a semicircular fitting sealed onto its top (Fig. 8); different jet diameters — 
could be obtained by using different fittings. Flow entered the box at one end 
; and exited upwards through a semicircular opening. To avoid flow separation, | 
the exit section of the half jet was rounded off smoothly. ‘he a 
Temperatures were measured with 44 thermistor probes (YSI 701, time constant q = 
= 9 sec, accuracy 0.3° F) which were individually calibrated and mounted at a3) ' 
the same level on a wooden platform supported on four screw jacks. Soon 
after the start of the experiment, dye was injected to observe the flow . The 
‘first scan of the surface temperatures was started when the dye front had gone 
"past a substantial area. After two or ‘three ‘surface s scans were taken, the » wooden 
platform was then lowered to record vertical temperature profiles. To ensure 
_ that some kind of quasi-steady situation was reached in the experiment, a surface 
scan was always taken at the end of the experiment. In all the runs, the temperature 
_ recordings of the last surface scan in the near field were very ¢ close to those ~ 
of the first few scans. The majority of the runs took about 20 min. The character _ 
of the near-field flow in the plane of symmetry was visually observed by due 
injection. Pertinent flow features such as the upper layer thickness, the location 
=... the jump toe (for the stable discharge configuration) were also noted. AF 
| summary oft run parameters and t results are in Table 


10.2 448 76.0 0.78 — 
(14.6) (24.6) (5.68) 
(4.35) 
| 
: 


‘SHALLOW WATER JET 


Continued | 


Unstable 
Unstable 

Unstable 


ony 


Unstable 


In general, | visual observations of stability and other features of the flow 
- in every way consistent with the temperature data. Fig. 9 shows typical © 
radial temperature transects corresponding to the possible near-field conditions. 
Vertical excess temperature profiles, A T(z,r), averaged in the azimuth, (6), — 
direction, are shown at various radial positions r/H from the s source. The vertical _ 
- homogeneity of the near-field temperatures, in the case of an instability, contrasts 
sharply with the vertical temperature profile of a stably stratified system in 
the In addition to the stable and unstable near-field 
conditions, a weak instability was also sometimes observed as a transition | 


by a an n internal jump with a conjugate “depth that exceeded ‘the water depth 
Bere, jump), causing slight re-entrainment of the upper layer mixed water. 


‘Comparison OF THEORY AND Dara 


4 Near Field Stability 
_ From visual observations of the upper ter flow (and confirmed mw the 
- temperature data) over a wide range of experimental conditions (F, ~ 10-80; 
10-30), the location of the jump toe is determined to be r,~ 0.6 


19.0 2.6 | 16,100 | Unstable 
a2 | 40 | 6,800 | stable | 25 4 
| 
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The stability (Fig. 7) is well grodicted by the stability 


“criterion. Under deep water, stable discharge conditions, the upper layer thickness — 


q in the surface impingement region is about 0. 0.1 H. This compares favorably a 
F with the theoretical prediction (Fig. 5), and also Rawn’s data (18) for ctr 


=] _ Under stable discharge conditions, th the th theoretical near field dilution is defined 


‘nahety as the ratio of the upper layer flow at the exit of the surface > 
impingement region, r= Tht to the discharge flow, = Q, /Q.- Experimentally, 


to the maximum observed excess ; temperature in n the upper layer, Sexe = 
AT,/AT,,,,- Predicted dilution contours are shown in Fig. 10 and compared 
to the observed experimental results. In general, good agreement is obtained, — 
the observed dilution being somewhat higher than predicted. This could a 
at attributed to the additional entrainment in the surface impingement region, and > 


og 


ope 


sgiedgal G. -——Comperison of Observed and Predicted Near Near Field Dilution 


ee the fact that in most experiments even the temperature probe closest to } 
the discharge may lie outside the surface impingement region. In the case of | 
an unstable discharge configuration, an eddying zone is observed close to the 
_ jet, , leading to considerable heat re-entrainment. The extent of this mixing zone 
can vary substantially depending on experimental conditions. For the range» 
of ' experimental conditions reported herein, the outer edge of this mixing zone, 
- r=r,., is located at 2-3.5 H. At the edge, the interfacial height changes sharply, 


a and a stably stratified counter-flow system is observed beyond this point. a. , 


7 _ The dilution capacity of the discharge is severely impaired if the discharge — 

parameters fall within the unstable domain. For this case it is reasonable to. 

postulate that the buoyancy driven exchange flow Q. of the counterflow system 
depends only on source flow Q,, B, =Q, g Ap, 


— 
—e 
| 


| 

Noting that = aD? and the dimensions B, 


Thus, ines of constant dilution i in the unstable domain follow pon of 2/ 5 Sil 
(H/D ~ ). A best fit of the dilution deta gives 


_ The preceding presentation demonstrates that the mixing and stability snie 


teristics of large diameter outfalls discharging buoyant water vertically into a 
is dependent on two dimensional parameters: F ,, a measure of are 
TABLE | of Design Parameters for Three Plant Sites 


Condenser 


flow, 


discharge 


7 (high water) | 7.9 Unstable 

0.9 (low water) 

LS water) 


Onofre 
Littlebrook 


: 

momentum and and a measure ure of the of the > 

water. An evaluation of typical conditions at existing cooling water discharge 
- sites (the 450 MW unit of the San Onofre Plant in California, and the Littlebrook © #. 

. and Sizewell Power Stations (500 MW each) on the Thames River estuary) 4 

is given in Table 2. Using the - stability criterion of Fig. 7, it would ‘appear — i 

ii? that all these « existing discharge designs operate in the unstable discharge domain. we ; 

_ Some field information (14) seems to support this fact. On the other hand, onal 

design details, such as the height of the port (which is one-half of the total ef 


depth at San Onofre), may alter the near field behavior, = 


sf 
Ms _The major conclusion of this paper lies in the observation that the shallowness a 
of a site has a pronounced effect on the dynamics of a discharge. Accordingly, - 


entirely different predictive techniques have to be employed in the design ae: : 
an unstable discharge, compared to the simple buoyant jet theories which can ie 
A secondary question that arises, namely, which discharge configuration ee: tg 
more desirable in actual design practice, is considerably more difficult to answer. : 4 
Given a heat rejection rate AT,,, a condenser flow and a chosen 


oat 
¥ 
a 


_ similar dilutions, S, are attained. Thus, no atten design strategy exists if 
_ maximization of near-field dilution is the only objective. However, other 
- ; considerations may favor the choice of a stable discharge design: For r example, Pid 
the prevention of a recirculating eddy and the associated large benthic temperature 7 
rises from an ecological viewpoint, or the fact that a stable discharge experiences 
additional dilution, beyond the near-field value, S, in the internal hydraulic | 
_ _Insummary, it is clear from the preceding arguments that shallow site conditions 
ne limit the dilution that can be obtained with single part discharges. Moving the 
| 4 outfall to sites with larger depths, or using a multiport diffuser design, are 
‘the obvious alternatives to avoid that constraint. Possible extensions of the 
present work to tidal situations and non-vertical discharges in shallow water 
4 The bulk of the research presented herein was “supported under the Waste 
Heat Management Program, 1974-1977, administered through the Energy Labora- _ 
; tory of the Massachusetts Institute of Technology, and was reported in a echnical 
- report (10). Additional work, notably i in the treatment of the surface impingement 
>. region and the radial al internal jump, _ and leading to a slightly revised criterion, q 
38 or Fig. was by the of Civil Engineering at 
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_ The following symbols are used in this paper: aad 


i jet width or diffusion thickness; 
ape F ie width of jet potential core in region n of flow « establishment 
densimetric Froude number; ge, tow 
4 Fe = free surface Froude mumber, 
= diffuser length; 
M = momentum flux; - 


temperature, ratio length to jump oc 


i 
| 
ta 
AT=T-1 
wien 


entrainment coefficient; 


ik ratio of temperature diffusion thickness to momentum diffu- a 

sion thickness; = } bal 


density; and 


center line values, value at critical section; 
entrainment, end of region of flow establishment; 
a = entrance to surface impingement region; 
= exit of surface ‘Surface impingement | region; forces te 


on 


om ie Te 


= 


q 


“Tes 
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"arg 
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Flows in alluvial ‘channel: bends" are among the the more difficult encountered 


of the channel bed. The balance between the centrifugal forces due to the 
vertical nonuniformity of the velocity profile, combined with the flow curvature, 
the shear stresses, and the radial pressure gradients, tends to produce strong 
_ secondary currents and a transverse slope of the channel bed, with an increase 

in the depth toward the outer (concave) bank. As a result, the outer bank 

_ This paper examines and compares existing - models for the prediction of the 4 
- steady- State transverse bed slope in channel bends, using laboratory and field 
data, and Proposes an alternative model. While traditional models are based 2 


3 helical flow pattern in n the bend, and its interaction with ‘the So 


g based on a a theory for the critical conditions for sediment ‘catrainment. The 4 
model gives good agreement between measured and predicted bed slopes. The - 
analysis herein may lead to: (1) An improved understanding of the problems 
involved in the description of bend bed topography; and (2) improved analytical > 
tools for development of channel-protection measures that are effective, environ- 

mentally acceptable, and ro justifiable. 
a Several analytical models have been developed for the ‘seiandtien of the 
0 -state transverse bed slope in a curved, erodible bed channel. Among — 
them are the models by van Bendegom (19), NEDECO (13), Yen (21), El- Khudairy 


—@), Ikeda (11), Allen (1), Ragiund (8), Kikkawa, Ikeda, and Kitagawa (12), 


a radially inclined bed. The equilibrium transverse bed slope calculated is that 
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hich the radial-plane comp -plane components, to the bed, o bed, of fluid 


in which sin B = radial or transverse slope of the bed s surface. The radial-plane 


_ force diagram is shown in Fig. 1. The radial plane is the plane normal to the 
Stina of the he mean flow. It is principally in in the calculation | of F, that the 

a Van Bendegom (19) determined F, by balancing the forces acting on. a fluid 
volume traveling at steady state through the bend. He considered the forces 
7 due to the longitudinal and transverse slopes of the water surface, and the 
centrifugal force. By assuming that: () The shear stress behaves according 
to Prandtl’s mixing length hypothesis; (2) the direction ‘of the shear stress 
corresponds to that of the velocity vector; (3) the radial and vertical components 
of the velocity vector are small composes with the longitudinal component; 


w 


alo 


... 
in which 7 = - depth-av averaged ‘mean velocity; p = density of the fluid; d a 
local depth of flow; r = local radius of channel curvature; and A= projected 
surface area of a sediment pe. oe Substituting Eq. 2 into > Eq. lyields 


al 


in which F7,. F, is the particle Froude number defined 


= law, van an Bendegom obtained — 


| 
fo 
F 
a 
| 
a 
ponents of Forces on a Sediment Particle 
y profile may be described by Strickler’sspower 
| 
I 


= mean ty of sediment par Ea. 2, van Bendegom assumed that: 
| Ay (1) The longitudinal drag, F,, on a sediment particle is determined by F, = a 
= _ Apg V*/C’, in which C = Chezy’s coefficient; (2) C = 50m'/*/s; and (3) 
: j the direction of F, deviates from the mean flow direction by an angie, 8, determined 
Van Bendegom’s s results were later used by NEDECO (13) with sin B replaced — s 
by the derivative d(d)/ dr to. ‘yield a differential equation which was solved 
6 obtain an expression for the radial profile of the bed. Yen (21), in his analysis, _ 
related F, in Eq. 1 to the radial component of the flow velocity and obtained 
; : = relationship between sin B and the maximum radial velocity. El-Khudairy’s i 
(7), and Ikeda’s (11) analyses are very similar to the analysis by Yen (21), 
except that the radial velocities used in the calculation | of F, were re computed 


int of vi van n Bendegom’ s ; results, oni a fo orce e balance : similar to that of v van n Bendegom — a 
(19). He obtained a relation between: (1) Transverse slope, sin B; (2) deviation — 
a angle 8; (3) bed shear stress, +,; and (4) mean particle diameter, D. if the ~ 


(8), Kikkawa, Ikeda, and that the 
longitudinal drag on a sediment particle is balanced by the friction force on a 
the particle from the channel bed: 

oP in which p = tan ®; @ = friction angle, and F,= lift force on the particle. 


Engelund (8) related F, and by the deviation: angle, using his 
experimental finding that tan $= 7 d/r, Kikkawa, et al. (12) related F, and in 
F., by assuming that the ratio between them equals the ratio between the radial, — 
V,, and longitudinal, V,, velocity components of the flow at bed level. V, — a i 
was determined from Fa bere F, Teplaced by a relationship between a 


anda a logarithmic profile for the mean flow velocity, 1 modified by the forced-vortex on . 
distribution, f,. Both Engelund (8), and Kikkawa, et al. (12) obtained Eq. 2; nae 
however, K is essentially a in Engelund’s while it is 
function of F,, and / V, in which t,/p ,in Kikkawa, et al. Engelund’s 
model no explicit on grain size or flow velocity. Bridge 
(5) interprets Engelund’s » as Bagnold’s (3) dynamic friction angle which relates 
to normal and tangential collision forces arising between sheared grains. According 
to Allen (2), it is then theoretically possible for the transverse slope to become > 
ine on grain size, at least under certain circumstances. Engelund’s 2 
: model was used by Engelund (8), and Gottlieb (10) to calculate the bed topography 
. co? in a given curved channel with given discharge and grain size. ee 
Zimmermann and Kennedy (22) related in Eq. | directly to ‘stream-flow 
“a characteristics by equating the moment, about the channel axis, of the radial — q 
= boundary shear stress (averaged around the wetted perimeter) to the moment > 
"produced by the nonuniformity (over de depth) of the centrifugal force. The balance 
is a simplified expression for the conservation of the flux of moment of momentum. 
: = Falcon’s (9) approach is an extension of Zimmermann and Kennedy’s at 


ry mH ed ctream _tunctionn tormulatinn a conctant ed 
| 


of the radial shear stress. The pressure gradient i in n his relationship i is determined a 
by equating: (1) The integral over depth of the moment of the pressure gradient 
- with respect to a longitudinal axis; and (2) the integral over depth of the moment é 


d 


in which V = point velocity at distance y from the bed; and n = dimensionless 
factor (which will be termed velocity- “profile exponent). The radial shear stress, 


TABLE 1. .—Summary of Key Relations in Reviewed Models 


an Bed-friction © _ Slope to depth / 

“Author, ‘or authors parameter 


1/2 


Ikeda, Kitagawa (12) = (2.75 — 4.35 v, 
Zimmermann, » (22) ve = (n + 1)/[n(m + 2)) (a/2) NF, 
Note: 1 ft = 0.305 m; a = ratio of projected surface area to volume for a a cote 
Ee divided by that for a sphere 


radial drag, F,, byt multiplying t, at bed level by the er" surface area 


According to Zimmermann and Kennedy (22), n° = the inverse of the overall 
‘friction factor f. In Falcon’s (9) analysis, n i is determined from the longitudinal 
slope, of the water surface, using the Darcy-Weisbach formula “modified 


Falcon aie Both and Kenney (22), and Falcon 0) obtained 
Eq. 2, with 


which c= constant. Zimmermann and ‘Kennedy’: s analysis, the 


_ c is one-third that in Falcon’s analysis, because Falcon neglected the moment a 


contribution of vertical shear stresses. (Table | is a ened of key relations =s 


+4 
He assumed that F, is 
eon me tO the radial nressure sradient and the vertical pradien 
| 
| — 
2 
= 
| 
4 
| 


_ The preceding models are based on — assumption that he sediment particles © 
shaping the bed are moving or being dragged over the bed in the longitudinal _ 
direction, with no velocity component in the radial or transverse direction. The — 
relations for the transverse bed slope are obtained from the expression of no 
net radial force on these particles. The grain size, D, used in these relations — 
2 is that of the moving particles. Considering that natural channel beds usually 
“consist of nonuniform sediments, and often have an armoring of coarser sediments 
at the bed surface, there may be a _ significant difference between the grain — 
‘size distribution of the bed load and that of the bed surface. To be consistent 


ot theory, estimates of the bed slope by the above models must be based a) 


on bed-load characteristics rather than the characteristics of the bed surface. _ 
The writer’s approach is an attempt to relate the transverse bed wae to a 
-_bed- surface characteristics. It is assumed that, in the state of equilibrium, the 4 
sediment particles shaping the bed surface are statically stable in both radial 


\ POINT OF SUPPORT. 


That is, the size “of the bed- is that the are 
just about to move in the longitudinal direction be rolling about their — 


about the Point —* ™ 


“of the submerged sediment. This analysis is s similar to. the one proposed by ah 
White (20). White, however, did not include the effect of lift force. Acconfing — 
to Raudkivi (14), Chepil (6) was one of the first researchers to consider lift 
force. Based on experimental results, Chepil concluded that for nearly tee 
‘grains F, = C, F,, in which c, = a constant of the order of 0. 85. al 


Chepil’ s result into Eq. 10, a and assuming that cos B = 1, yields 


2.—Vertical Longitudinal-Plane omponents TParucie 


Eq. 11 may also written 


in which +, = critical bed shear stress; D. = diameter of the particle whose 


is impending; and = dimensionless factor. by White (20), 


“function of a critical critical boundary Reynolds number, defined as be 


— inwhich V,.= critical shear velocity defined by V,. -/p ;and v = kinematic 
viscosity of the fluic. Shields’ results, which as a plot of 1./[(p, 
— p)g D.) versus R acs Were approximated by a curve (the Shields curve) by ‘ 
i. Rouse (15). The curve is seen in Ref. 17. The main features of Shields’ curve 
are a laminar-boundary layer region, < ~ ~2, in which @ varies as 
a transition region, and a region of fully developed turbulence, R,. > 400, 
_ in which 6 is substantially constant at a value of approximately 0.06. The ‘transition — ke 
_ curve has a minimum at R,. = 10, where @ = 0.03. Shields’ curve is valid 
_ only for flat beds; therefore, T. in Eq. 13 is the critical flat-bed shear stress. _ 
_ The critical shear stress, t., may also be related to a critical mean flow 


_ velocity, F, ~» by means of the Darcy- -Weisbach friction factor, if  (flat- -bed friction 


off,= 8K? and Eqs. 13 and 15, yields 


g D.— — 
i in which k= Karman’ s constant. Eq. 16 may be reduced to 


e particle Froude number based on ond D.; 


that Falcon’s relationship, Eq. 1, holds, Eq. 1 ee 
P,-p)g¥ sing . 
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which Fp. = th 
| 


n which - A = volume and projected surface area, respectively, of a 1 sediment 
tha-vid 


@) 


in which a = 1 = ratio of projected. surface area to volume for a a sediment particle 


(23) 
| a which is similar in form to Falcon’s relation, * 9. Thus, while Falcon’ s K- value. 
7 is a function of the overall bed resistance and the size of the bed-load particles, — 


the writer’s K- value is a function of the flat- -bed Tesistance and the size of 


oF Prepictive Capasiuty 


25 

in which S,;= d(d)/dr a at the center 29 


Comparison by Laboratory Data. _—Initially, the .e performance of the models 
was tested by means of Zimmermann and Kennedy’s (22) experimental results. oe 
_ Fig. 3 (a, 5, c, d, and e) shows a comparison of measured and predicted slopes, — 
the ordinate being the measured values of S,/(d./r.), and the abscissa the 
values determined by the formula for K, using measured values of d, D, and 
3(a) shows the van Bendegom model 


by of 2.7%) it Pigs, ged] 
in which K gee x n7 N, 
| transverse bed slope to 
portional to the ratio d/r: 

oportionality, K. Therefore, 
_—— =} a comparison of their performance may be limited to a comparison of their — 
capably of predicting the local slope, which may be taken at the center line — 
of the cross section, in which r = r. and d = d.. Eq. 24 is then expressed a 
= 
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between measured slopes and slopes predicted by this model (calibrated 
by a factor of 3.55) is quite good at smaller values of S,. Systematic deviation 
of the measured values of S; from the predicted is seen to occur at e larger 3 
AS seen in ‘Fig. 3(b), ‘the correlation between observed slopes, and ‘slopes 7 


pootictes by Engelund’ s (8) model, is fairly small. Since the dynamic friction 


Mage 


= _3.—Comparison of Tran 


Kikkawa, Ikeda, and 8; id) Zimmermann and or Falcon’ ‘and 


coefficient, 4, is more or less constant, Engelund’s model always predicts 
the same value of S,/(d./r.). The reason for the poor agreement may be that 
_ Engelund’s model cannot be applied to Zimmermann and Kennedy’s data. The 
_ flow in the test section was a fully developed bend flow, and Engelund’s model — 
was never validated (or intended) for such a flow. -_Engelund’s s model seems ~ 

to perform very well for weakly meandering flows, in which the secondary 7 


| 
erse Bed Slopes Measured by Zimmermann and 7 


_ To apply Kikkawa, ¢ et al. (12) model, the function M (Table 1) must be evaluated a 
at the center line of the channel. Using: (1) A value of one for the forced © 
; vortex distribution (which is the value to be expected at the center line); (2) 
pe Vi/V=avVf/8 = 1/(nV8); and (3) a sheltering factor of 0.592, Kikkawa, 
etal., K-value is reduced to F ,,(1.628-0.912 /n). The value of n for each experiment | 
_ was determined experimentally by Zimmermann and Kennedy, using the side- -wall 
; a correction procedure of Vanoni and Brooks (17). As seen in Fig. 3(c), 
agreement between observed slopes and slopes predicted by Kikkawa, et al., 
- (calibrated by a factor of 0.76) is very good. Fig. 3(d) compares the observed — 
_ Slopes and the slopes predicted by Zimmermann and Kennedy’s / Falcon’s model. oy 
_ The agreement is seen to be somewhat smaller than when Kikkawa, et al., 


“current is in its ant and ti by residual currents 
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» 

4- 
To apply the author’ s model, my must be be evaluated "using 17. ou 
is a function of R,, which, in turn, is a function of n, (at least, at lower zs 
values of R,), the “value of n, has to be determined by iteration. The value 
is found to be approximately four i in Zimmermann and Kennedy’s s (22) 


to Shields’ curve, these values of R, ‘correspond to values of 6 “of 0.04 and ri 
0.03, respectively. As shown in Fig. 3(e), the agreement between observed slopes 
and slopes predicted by the writer’s model is good. The writer’s model was 

= without introduction of a calibration factor. MECH) 


is against field data. A ‘measerement: ‘program v was carried out by ‘the S. 


_ Geological Survey, Sacramento, on a bend of the Sacramento Kiver be vet 
river miles 187.5 and 189.5. Fig. 4 shows a schematic plan view of the site. 


| Six cross sections of the bend were surveyed i in April and May 1979, and 
- March 1980 to define the channel geometry at two different times for two 
different flow rates. Also, detailed velocity measurements and bed-material 
samples were taken at each | cross | section to determine the variations in the 
velocity profiles and the size of the bed sediment. The direction and — 
of the flow velocity were measured over 4-10 verticals across each section, 


using a Marsh- eye flow meter. In each vertical, the pow st = 
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FIG. 5.—Cross Sections No. 3 and 4 of the Sacramento River Study Ben Bend at Low | 


9.11 100 Cubic Feet per Second (260 Cubic Cubic Meters rs per Secot Second) 


each of the vertical. Based on these measurements, a depth-averaged 

velocity, V, was determined for each vertical. The cross-sectional average 

velocity, V,, was determined by dividing the known (gaged) flow rate by the 

‘measured cross-sectional area. The bed samples were taken by a “Digging 
“sampler” (USBM-54), which samples to a depth of 2-3 in. below the bed surface. 

_ The n- -values were determined on the basis of the longitudinal slope, S, of - 

the water surface, using the Darcy-Weisbach formula modified by Falcon 2 


| The procedure for the determination of n may be seen in Falcon (9). The slope 
th the was measured and found to be of the order of magnitude 


| 
off 
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of 0. 0003. According to Brice (4), the average — of the bed between river - 


the seca of the best fitting line. The cross-sectional average value, | Nyas Wi was 

determined as the arithmetic mean of the n,-values in the cross section. The - 

transverse bed slope was measured oo oo the channel center line; however, the 
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oh 
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Kitagawa’ 8; (d) Zimmermann and Kennedy’ or Falcon’ 8; ‘and (e) the Writer’ 


= sein to be practically constant across most sections. (A summary 
7 of the data and computational results is presented in Table 2. Fig. 5 shows 

Cross Sections No.3 and4atlowflow.) 
4 Fig. 6 (a, b, ¢, d, and e) shows a comparison of measured and pasties. 
4 ordinate | is the measured value of S,/(d, /r. and the abscissa. 


is the value calculated by formula measured cross-sectional 


i | 

| 


sont 


gilt. TABLE 2 —Secramento River 


Flow rate, jin Cross-sectional average velocity, Cross-sectional 


cubic feet area, in V,, in feet — depth, 
per second square feet per second in feet 


ae D,,, and S, respectively. As will be seen in Fig. 6(c), Kikkawa, et al. 4 
- model (calibrated by a factor of 2.0) gives good agreement between i 
and predicted slopes. The agreement is also good at smaller | values of S, when 
Zimmermann and Kennedy’s (calibrated } by a factor of 5.1), or Falcon’s ‘model 
(calibrated by a factor of 1.7), is used, as will be seen in Fig. 6(d), oF when 
van Bendegom’s or Engelund’s model is used, as will be seen in Figs. (a) 
a and 6(b), respectively. Systematic deviation of the measured values of S, from 
the predicted is seen to occur at the larger slopes. Although Falcon (9) recommends 
- using D,, for the calculation of F,,, Falcon’s model was found to give better im 
‘agreement between and slopes when Deo was used. ‘Using 


a ‘surprisingly good poll of the Ny, Was “calculated by Eqs. 20 
and 17, using 6 = 0.06, x = 0.4, and the measured values of F, for Foes 
The calibration factor (2.43) may reflect the difference between F, and F,., - 
_ or the fact that « and 6 are not Constant in movable bed channels, or "both. Po 
“ So far the ‘comparison of the models has been ‘made on the basis « of their 
~ capability of predicting the transverse bed slope at the center line of the cross 
sections. Strictly, an evaluation of the models should also include a comparison 
of the complete solutions in the form recommended by the respective authors. — 
== (8), Kikkawa, Ikeda, and Kitagawa (12), and Falcon (9) =n ol 
writer recommend ‘using constant slopes across the sections. However, it is” 
easily verified that the trend seen in Fig. 6 remains unchanged whether -etiedhel a 
solutions or constant slopes are used with the present data. At small and moderate 
slopes, there is little difference between the predictions. In all all — the 


Section 
 numbe 
24,900 6510 3.82 
28,4000 7 _ 
9,100 | 4,000 | 2.28 ‘15.20 
6H | 26800 =| 670 
Note: 1ft= 0.305 m;lin.= 254mm 
| 
it 


_ TRANSVERSE BED SLOPE 


 Cross- 


diameter, Froude | Radiusof | Transverse Profile 
Exponent 
“number, ‘4 curvature | bed 


ta 


key” parameter is the whether it be that of the bed load ‘that 
of the bed surface. As will be scen in Table 2, the variation of D,, throughout — 
the bend is very significant. No technique or model is yet available that cm 
simulate a variation like that in a given bend, given the flow data. Instead, 
a rough estimate of Dog ‘may be obtained on the basis of a combination of : 
3 Manning-Strickler type of equation and Eq. 17, giving a relationship between 
_ the Froude number, F = V,/V gd,, and the relative roughness, D,,/d,. Fig. 7 
7 shows a plot of F versus D,,/d, based on the data from the Sacramento 
River. Although the theoretical curve approximating the points cannot be a 


straight line, due to the variation of n,, through the bend, it is seen that a 


straight line of slope 3.9 may | be si sufficiently accurate for a rough, tentative | 


es timate of so, given the values of V, and d,. The line is given byt by - pean 


Note that this $ equation is based on only the data from the Sacramento River 


Application. .—The procedure for the application of the writer’s - model 
summarized as follows: Given the cross-sectional average values of depth, d,, 
and velocity, V,, the cross-sectional average value of D,. may be estimated 
using Eq. 26. (Eq. 26 is tentative and may apply only to the Sacramento River 
data.) The cross-sectional average values of n, and N, may then be determined 
: by Egs. 17 and 20, respectively. Knowing the radius of curvature at the — 
section in question, sin B is obtained by substituting n,, and Nye into Eq. 21. 
- Based on the Sacramento aed bend data, K in Eq. 21 should be multiplied 


Bend Flow Data bux 
< 
(11) 
(2.81 3920 0.010 | 4.0 | 63 
| 1,800 0073 22 | 96 : 
| 
- 
: _ The differential equation for the transverse bed slope, which is obtained by » 
: _ replacing sin B in Eq. 21 by d(d)/dr, is not readily solved since both n, and } ; 


4 are functions of me mean velocity and grain size, e, both of re vary in a : 
yet unknown manner across the channel sections. However, the field dataindicate 


that a ‘reasonable ‘first approximation may be be sin B 
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FIG. 7. —Relative Roughness as Function of Froude satian for Sacramento River r - 
able approximation over the central section of the channel —_ bank effects 
we small; and pad the local longitudian! slope v varies as il in n which 


of mean over the central channel section may ay be estimated by 
ad.) \r 
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in which (d/d.) may estimated by Eq. 21; and V.= mean flow 
igs study has: (1) Reviewed existing models for the "prediction of the 
steady-state transverse bed profile i in channel bends; (2) formulated an alternative 
_ model, and (3) compared the models, analytically and by laboratory and i 
Existing ing models are all based on the equations of motion for the sediment 1 
a og in the bed load. The writer’s model is based on the assumption that — 7 
the bed | topography is determined by the critical conditions for sediment = 
entrainment. All models that relate the local transverse bed slope to mean tow 
eet. which includes the writer’s model, are alike in the sense = 
they all predict the local transverse slope to be proportional to the ratio between 
tess d, and radius of curvature, r. The Pon on between the models io 
_ Engelund’s model predicts K to be essentially c constant throughout the channel ail 
_ bend. In van Bendegom’s model, K is a function of the bed-load particle Froude - 
number, F,,. The models by Kikkawa, et al., Zimmermann and Kennedy, and © 
Falcon, predict K to be a function of both the velocity- profile exponent oa 
(corresponding to the Darcy-Weisbach friction factor /), -andF,,. While, Kikkawa, 
etal. model predicts K to be linearly dependent on F,, Zimmermann and | Kennedy 
> and Falcon’s models predict K to vary with F,, to the second power. In the 
_ writer’s model, K depends on only the flat-bed velocity-profile exponent, n, 
4 (corresponding to the flat-bed friction factor f,). Since, according to the writer 
theory, n, is proportional to the bed-surface particle Froude number, F,., at 
fully developed t turbulent flow, the writer's K-value may be expressed 
in terms of Fics thereby reducing the | problem 1 of predicting the . transverse bed | 
; slope to, essentially, one of predicting the overall grain size distribution of > 


the bed — given n the flow rate. The writer’s K-value is largely proportional | 5 


aa _ Using proper calibration factors, all the models give reasonable agreement 

between measured and predicted slopes at smaller transverse slopes. Systematic 7 
- deviation | of measured slopes from the predicted occurs at larger slopes, exc ‘ 

when the models by Kikkawa, et al. (12), and the writer, are used. _— a 


: . a The study has shown that the key parameter in all available models i is ; grain 


- size. The prediction of the complete bed topography of a given bend for given 
flow data requires knowledge of the spatial variation of the grain size. Unfortu- _ 
nately, ‘fo technique or model is yet avaliable’ that this veriation 

Support for this study was provided by the U.S. Army Corps of Engineers’ 
Sacramento District, California, under Contract No. DACW0S5-80-C-0083. 
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om The rhe following symbols ar are neat in in this Paper: 


(= projected surface area of sediment > 


= lift force coefficient; ; 


constant; sss 

coefficient; 
= particle diameter; 
= mean sediment diameter; 


ng and 
An of suck 2 Hew. iv without) | 
= flow depth at channel center rline; genic al srikingly 


Froude number = = gd; 
iM Froude number V/{ [@, - 
F longitudinal drag component; 


radial drag component; curren envy 
om 


factor; 
= filat- bed friction factor; 
« They 
acceleration of gravity; _ 
Tatio of transverse slope. to d/r; 
functions of friction factors; 


the 
al eniry | the to? 
radius of curvature; radial coordinate; . 
mean channel radius; 
= longitudinal channel slope; 


depth- the flow velocity; nit The yw 

Projected area- -volume for normalized by that for 


6 = Shields’ parameter = 1,/ [(p, — D 

= fluid density; 
particle density; 

shear stress; 
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angle of repose of submerged sodiment; ont 


"related to to Darcy- -Weisbach’s flat- -bed friction 
radial direction; an ond 
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SHEAR ‘STR ESS DISTRIBUTION AT 

By Bishnu P. Das,’ M. ASCE and Ronald D. Townsend 


A channel constriction results i in contracting g and ae flow through the 
_ gap. The flow is predominantly curvilinear, accelerating and decelerating while — 
- negotiating the constriction. An analysis of such a flow is incomplete without _ 
reference to the flow through a bend where the general pattern is strikingly 


_ Flow Through Bends. —For uniform rectilinear flow in a prismatic channel, % 
the boundary shear stress distribution has been studied in the laboratory. Both 
the longstream isovels and the shear exhibit a symmetrical uniform pattern. 

if the flow negotiates a bend, generation of secondary currents before entry 
N the bend create asymmetry both for the velocity and shear distribution. a 
Numerous researchers (7,15,16,18) have studied the flow through a single bend — 


‘for channels having both rectangular and trapezoidal- shape cross sections. They | 

_ have observed specific locations, such as close to the inner bank at entry and 
outer bank at exit to monitor maximum bed shear stress. The locus of the 
/ maximum shear filament shifts from the inner bank at entry | to the outer bank > 
at exit, more or less following the locus of the maximum velocity filament 
“through the bend. It is postulated that the secondary flow generated | by the the : 


stream curvature produces the distorted shear pattern. 


<~a _ Nouh and Townsend (15) have advanced the study by « developing for a a specific 
Ss bend geometry, nondimensional contours of shear maxima treating the uniform 
shear stress in the straight reach, upstream of the bend, as the reference value. — 
Erosion of the outer bank is generally considered to be a consequence of 
the fh high shear stress occurring there. Yet, although the inner bank at entry e 
experiences high shear, deposition occurs forming a point bar. This is explained 
by Yen (19) as follows: ‘‘. . . that deposition and scour are not directly related 
to the magnitude of the bed shear. Rather the secondary flow directed inward 
- near the bottom which carries the sediment particles from the outer side of 
a bend to the inner side of the following bend determines the bed ss __l ig 


‘Superintending Engr., (Civ.) Dept. of Irrigation, Orissa, India. 
Prof, , Dept. of Civ. Engrg., Univ. of Ottawa, Ottawa, Ontario, 
~ Note.—Discussion open until May 1, 1982. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, — 
ASCE. Manuscript was submitted for review for possible publication on September 26, 
1980. This paper is part of the Journal of the Hydraulics Division, Proceedings of the 

(ete Society of Civil Engineers, ©ASCE, Vol. 107, No. HY12, December, 1981. 
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presen’ have been conducted extensively, actual measurement of shear stress | 
along the constriction boundary and downstream has not been reported. The 
analytical derivation of maximum scour depth, however, has been made by 


— GRIGINAL, Beo Lever 


ened (c) m= 0653 § 
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FIG. 1 —Velocity Distribution wi with h Contraction 


of frequent necessity, for siting a hydraulic structure. bridge on an ‘alluvial 
= is generally provided with a constricted waterway keeping the well-defined - 
~ bankful Stage within the bridge vents, and closing the overbank flow by the = 
~ embankment approaches. Even in deep valleys, abutments frequently protrude 

* the main stream causing partial constriction of the deep channel. Often — 

| Soedams built with loose fill are pushed into a stream from one bank to 
_ a work area which can obstruct up to 20-30% of the waterway. Studies 


by Lane (11), Kindsvater and Carter (10), Valentine (18), Izbash and Lebdev 


- Flow Through Constrictions.— Although laboratory studies on maximum depth a 
| 
one or both banks. Such closures, which result in an abrupt constriction, are 4 : 
| 
— 
= 
| 


e been ‘directed mostly towards derivation of the backwater rise and scour 

- depth at contractions. Studies of scour at spur dikes, which are abrupt constric- 
tions, have been made by Awazu (1), Garde et al. (5), and Gill (6). Although © 

_ these studies did not report boundary shear measurements, Das (2,3,4) has 

_ experimentally studied and analyzed for end-dump constrictions the stability 
of closure material forming the body of the dam, and the associated depth 
of maximum scour downstream, by relating them to the maximum ~ 

consequent to that imposed constriction. 


fi. 


4 


This University of Alberta study (3) also and 
related boundary shear stress distributuions on both the advancing face of an 
end-dump constriction and on the channel bed in the vicinity of the constriction. 
Sample velocity and 1 boundary shear stress distributions from the study are 
To examine further the applicability of such laboratory studies to field problems, 


Similar investigation was performed at the University of 
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a deep central section flanked by relatively shallow berms. The redistribution 
; of shear stress occurs as a result of flow on the berms suddenly being channelled 
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FIG. 3.—Variation of Nondimensional Shear Maxima with Contractio Ratio 
‘examined in the study. A more detailed discussion of these 1 ton ts is containe 


A uniform flow implies constant energy y expenditure in the direction of 1 of th the 


flow and, accordingly, the uniform shear stress on the bed 7, is given as _ a 


boundary shear 
in a channel formed in cand and having 2 cection concictine o 4 
; s Fig. 5 shows the basic shape of the cross section and also presents scour E 
patterns relating to two flow stages, and three different contraction ratios 
MELOCITY PROFILE 
| 


is which 9 = unit cunts’ of water; h = uniform flow depth; and s = slope 
a as the same for bed and water surface. $= va 
Introduction of a constriction will result in the flow contracting sharply through — 
the g gap. The flow against the constriction face will be curvilinear and will 
Ps accelerate up to the vena contracta. It would then decelerate, the live stream 

Pr seen beyond the vena contracta to diffuse gradually into the full downstream — ; 
width of the channel (Fig-6), 
vs In Fig. 6, a channel of width B is constricted by end dumping to an average 
width 5, and the flow further contracts to a width b. at the vena contracta. 
The contraction ratio m is defined as (B — b)/B. by 
Iti is obvious that the boundary shear through the constriction would be motifies 
depending on the geometry of the constriction, the flow characteristics, the — 


“fluid characteristics, and the bed material characteristics. 


De Tere Oe 


three of pebbles of (6 640 mun, 1900 id 

‘The followit ng parameters define the problem: B = approach channel mean 
width; b = = mean gap width for end-dump constrictions; V = approach channel a 

- mean velocity; h = approach channel mean depth; p,, = specific gravity of 

= water; g = acceleration due to gravity; v = kinematic viscosity; d = median 

size of bed material; o,, = geometric standard deviation of the bed material; 

a p, = specific gravity of bed material; and ee max» = Maximum local shear 
against the advancing face and downstream bed. 
As the flow against the closure material will be fully kinedintic 
4 viscosity v is considered unimportant. Similarly, p, and p,,, being constants — 
_ for noncohesive sediment and water, are ignored in the analysis. — te pcos 


Dimensional analysis indicates that ihe nondimensional increase in stress 


| 
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_ The principal nondimensional terms are considered to be m, the contraction — 
ratio, and V/V gh, the Froude number of the approach flow. The terms d/h = 
and B/h are not significant as determined from the experimental findings. a 
_ Observations by Ippen and Drinker Ms § and Replogle and Chow (17) also attest 
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FIG. 5-—Nondimensionl She Shear Distribution Overbank Channel Constrietic ns 7 
Awazu (1) studied scour around spurdikes concluded that with an increase 
in the contraction ratio, T,,,, increases ia and the nondimensional 


shear maxima is related to m in the form att 


in which j, and j, : = coefficients primarily | dependent onF, 
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In view of the | limitations of past investigations 4 


_ in this area, it was decided to proceed with a laboratory study in which actual 
ma shear stress distributions in a model channel constriction were deter- 2 
mined from velocity measurements. The study used Froude numbers and ; 
contraction ratios | in the range 0. . 10-0. 50 and 0. 10-0. 90, respectively. palo A ie 


Flume, Material in Bed and Dam. were carried out 
in a flume 60 ft (18.3 m) long, 7.5 ft (2.3 m) wide by 4 ft (1.2 m) deep with © ; 
7 a horizontal bottom resting on the laboratory floor. Of the 4 ft (1.2 m) depth 
available in the flume, the bottom 1.5 ft (0.46 > m) filled with the bed material © a 
over a 50 ft (15.25 m) length. A head tank, an upstream stilling pool, and a 
. a drop box downstream of the tailgate occupied the remaining 10 ft (3.05 m) 4 
length of flume. The beginning of the 50 ft (15.25 m) long ‘‘sediment’’ section 
was designated as station ‘‘0,’’ with downstream stations (at | ft (0.31 m) spacing) i 
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FIG. 6. —~Definition Diagram for Flow and 

A 30 ft (9.2 m) stretch from Station 13-43, where uniform water surface — 
slope could be obtained, was selected as the test reach. End- -dump dams with 
three sizes of pebbles of D,, equal to 6.60 mm, 18.00 mm, and 25 mm were . 
built symmetrically into the flow at Station 24. Three sizes of sand of 1.20 4 
am, 0. 60 mm, and 0.25 mm median diameter were used as representative ae 


Measuring Devices. —The comprise the fol- 

eee point gage for measuring the water surface elevation and the static | 

An electronically-sensed propeller meter for The 


| 
| 
| 
| 
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could we taken ‘the normal to the 
curved advancing face). The propeller meter was connected to an autograph 7 ] 
recorder which permitted measurements of velocity to the nearest 0.01 fps (03 
3. An electronic bed probe (working on the optical principle) was connected 
to an autograph recorder to measure bed elevations to the nearest 0.001 ft 
~All the « equipment was meted: on an instrument carriage that moved on 
7 rails located on the side of the flume. 
Test Series and Flow Parameters.—Since the Froude number of flows in 
mild- sloping natural chenaets rarely exceeds 0.5, even vor flood conditions, 


carried ~ For each 1 run, the boundary shear maxima and pattern were obtained 


TABLE 1 of Experimental Program 


fication numbers | in millimeters 
and sand 1.2 mm Both side closure 
a sand 1.2 mm 7 Both side closure 
sand 1.2 mm | Both side closure 
sand 1.2mm One side 
sand 0.60 
0.25 mm | 


of actual closures 


coal and large 
pebbles 


Procedure. —Equilibrium and steady flow conditions were produced by fine 
adjustment of the tailgate control. Closure material dumping from both sides — =~ 
symmetrically, or one side, was arbitrarily stopped to correspond to a particular — 
- contraction ratio. . The water surface profiles, the live stream boundary and 
* velocity profiles, normal to the boundary, were then recorded starting at the = 
section of maximum backwater upstream, and ending where the flow widened 
nor or less uniformly to the parent channel width. The readings were taken 


fairly close in plan to the end-dump face and at | ft (0.3 m) intervals elsewhere. 


The bed profiles were recorded by the ‘‘optical’’ probe at suitable intervals 

i. as to complete the bed configuration u upon which the water surface configuration a 
and velocity pattern could be superimposed. This correlated data provided a 


detailed description of the flow w field. 
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= shear. 


a Contraction and Expansion of Live Stream ast —The large- scale 2 
- ‘pattern was visualized by dye injection and also by dusting confetti on the — : 

water surface. The local angle of expansion, ,,, is smaller than the local 

angle of contraction, ,, at the vena contracta (Fig. 6). Ina 
constriction the contracting stream maintains a more or less vertical edge. For 
the end-dump constriction, the contracting stream edge is more or less vertical _ 

_ immediately upstream and downstream of the dam, but hugs to the end-dump 

face in between. It is therefore obvious that the geometry of the constriction _ 

has a decisive influence on the coefficient of contraction, as weil as s the flow ‘ 


A nondimensional plot of contraction and expansion of the livestream is shown 


in Fig. 7, which also presents, for the purpose of comperieen, the relationship © 


Velocity Distribution through Constriction. _Fig. typical velocity 
a - distribution through the gap at increasing: contraction ratios, and Fig. 8 shows 


unrestricted channel flow had a mean due of 1.15 fps (0.35 m/s), wllstens 7 
the flow negotiating the constriction attained a maximum point velocity of 1.81 
fps (0.55 m/s) at a contraction ratio of 0.25, increasing to 2.05 fps (0.63 m/s) 
at a contraction ratio of 0.775. The interesting feature, however, is that although © 
the contraction ratio increased from 0.25 to 0.77, the overall mean horizontal - 
velocity through the gap remained around 1.50 fps (0.46 m/s) because of 
considerable local bed scour (Fig. 1). This phenomenon of mean a 
velocity in the ‘‘gap’’ region remaining more or less constant is, , to a large 
extent, supported by field measurements of flow w through bends, where 4 ‘ve 


The maximum velocity filament hugged to the livestream boundary, akhough 
4 the velocity decreased with time along with the progression of bed scour. It 


is is pertinent t to mention that neither the mean velocity through tl the gap nor the 


4 of scour) had satisfactory correlation with either the stability. of pape ns ol 


On. Advancing Face. —The contracting flow accelerating between Sections: 1 
| and 3 (Fig. 6) was a case of boundary layer flow under an exceedingly favorable 


pressure gradient. The shear stress on the boundary was, therefore, determined © 


97093 
g ch experimental Gata ODtained for test series B, C, and Df 
, nts the measured variables pertinent to flow description and = J 
7 the longitudinal! velocity distribution at one contraction ratio. This latter pattern 
| 
q 


Maxie | 2d) 
velocity | | = 
| 
Length ben traction: 
ofex- | 
tracta, | traction,| pansion) 
in feet | in feet | in feet on face on bed 


320 | 5.90] 067 
0840 


2.35 | 1.80 


| 4so | 100 | | 2.50 | 


lll 
TABLE 2—Basic Experimental Data 
rt 
mwiverr 
gap 
width, 
| in feet 
| in feet 
. ___@ Series B, Q = 1.56 cfs, B = 7.5 ft, h = 0.18 ft, V= 1.15 fps,F=0.50 
40 190 3.80 0.90 | 3.50 | 6.00 | 430] 0.47 
| 3.80 195 | 250 | 150 | 5.00 | 7.80 | 5.70] 
8609 |: 3.40 | 2.00 | 260 | 1.70 4.20 | 620 | 480] 054 
m | 2.40 | «2.06 | 1.90 | 2.00 | 6.00 | 13.40 | 9.70] 068 
5.65 | 1.80 | 5.50 | 080 | 3.00 | 215 | 160] 025 
m 12 | 400 | 185 | 3.70 | 1.00 3.50 
| 13 | 330 | 195 | 3.20 | 050 | 3.60 | 9.00 | 750] 050 — 
614 | 2.60 | 1.00 4.10 | 12.50 | 8.90 | 0.65 
(220 2.07 | 4.80 | 13.60 “10.30 O71 
| 1.70 | 230 | 130 | 180 | 15.80 | 12.10 | 4 
| 450 | 190 | 420 | 080 | 280 | 860 | 3.10 040 
| 3.10 | 190 | 250 | 090 | 390 | 670 565] 059 
| 230 | 155 | 110 | 460 | 13.15 | 8.60 | 0.69 
| 1.80, 2.23 2120 | 1.80 | 470 | 15.60 | 1210] 076 
| | 243 | 100 | 1.80 | 5.70 | 17.50 | 1250] 080 
——____@) Series C, Q = 1.00 cfs, B = 7.5 ft, h = 0.38 ft, V=0.35 fps,F=0.10 
080 | 3.80 | 1.00 5.00 | 2.80 | 2001 043 
| 4.30 | 0.80 | 3.80 1.00 500 | 200] 200) 043 
| 2.56 2.20 090 | 6.00 | 6.20 (6.50 | 0.660 
203 | 140 | 200 | 230 | 5.00 | 13.50 | 10.20 073 
5.05 | 063 | 400 | 1.10 5.00 | 185° 033 
4 28: | 3.90 | 086 | 3.00 | 150 | 5.00 | 475 | 329 | 
«| 2.45 | 140 | 200 | 1.70 | 
30 | 1.20 2.07 | 7.00 
. | 2.20 | 085 | 220 | — | 17.500 0.87 
| 450 | 073 | 3.80 | 1.10 | 600 | 2.65 | 1.95 | 0.40 
| 3.20 | 2.40 | 6.00 | 9.00 | 
34 | «1.70 | 153 1.50 1.40 | 7.00 | 1410 | 1065 | 0.77 
35 | 090 217 0.88 1.50 | 650 | 18.20 | 13.50] O88 
0.75 2.27 0.70 | — | — | 30.00 | 22.00] 090 
_ (ce) Series D,, Q = 1.70 cfs, B = 7.5 ft, h = 0.26 ft, V = 0.87 fps, F = 0.29 Sie _ 
| 1.00 | 550 | 3.30 | 210] 030 
4 40 415 | 147 | 3.60. | 1.00) | 6.00 6.15 | 420] 045 
4 | 7.50 | 9.70 | 5.80] 068 
— Series D,, Q = 1.70 cfs, B = 7.5 ft, h = 0.26 ft, V=0.87 fps, F=0.29 


CHANNEL ANNEL CONSTRUCTIONS 


160 | 650 | 8. 6. 0.73 


and eliminating constants _K, — 

“manipulation gives 


- 


was nondimensionalized, with reference to ‘the normal boundary shear stress, . : 
of the unconstricted channel, and plotted (Figs. 2 and 3). 
On Bed of Channel Around Constriction.—Boundary shear on the bed ne’ 
determined using the same two- -point velocity measurement procedure and was Ss) 
plotted, nondimensionally, as contours of +,/7, extending from the maximum | 
backwater section (upstream | bal the constriction), to the of the downstream 
The locus of maximum shear was found to have a shape similar to the livestream — 


3 boundary. The shear stress in the backwater reach, and even on the bed a 


upstream of the gap, was less than the normal shear value; whereas in the 


_ many times the normal value. The maximum point shear occurs on the advancing» 
face and decreases gradually as the vena contracta is approached and near 
“Cee the maximum bed scour occurs. In contrast to the highly Se 
“(ares -dimensional) flow against the face, the flow field at the vena contracta 

_ was two-dimensional, though accelerating faster. Apparently, because the ae 


severely contracting flow region downstream, the stress very quickly rose to : 


at the point of maximum contraction tends to become rectilinear, the shear 
Stress intensity decreases downstream from same. These features parallel ob- , 
servations ns by many Tesearchers investigating flow through bends, where the 


a The nondimensionalized shear Has Aad , both on the embankment face and 


at the vena contracta, are plotted against the contraction ratio m, (Fig. 2) from 


| 
45 2.05 
q 46 | 1.60 
47 230 1.10 | 650 | 965 | 5 
48 | 130 8.00 | 13.20 | 9 
: 49 | 160 | 223 | 080 | 200 | 7.60 | 1610 | 12 _ 
:1ft=0. 
Note: 1 = 0.305 m | 
| 
| 
which +, = boundary shear stre 
 6=—S__sy, and y, measured, vertically, abov 
| 
q 
4 


= 


Re 


in? 


_ 


4 £06 | 


the shear increase exceeded 15 times the normal at contra 


1706 
| 
tm 
| 
. 7.—Nondimensional Profiles of Livestream Boundary 
_ which an exponential relationship is seen. As will be noted from Table 2, although 
_ the maximum velocity against face was less than double the mean velocity,™ 
ction ratios exceeding 


— 


sien, Ts can be 
Pe box wi ooh 


Se 


wo 


Note* Maximum poi point “Velocity: Contour sie? 


wth. £0. of ff 


AG. G. 8. —Typical Longitudinal Velocity Distribution 


a As mentioned earlier the preceding sections refer to a major study conducted — 
at at the University of Alberta (3). A second study, employing a 4 m long physical 
model of a typical river channel, recently was completed by the senior writer 
at the University of Ottawa. wd Now ‘ot 
The model channel in this recent study was moulded in a 10 cm thick sand 
ue (ds0 = 0.7 mm) i in a ce flume. As shown in Fig. 5, the ‘fiver ¢ cross 


é 
A 
| 
| 
| 
| 


‘Flows, to depths 15 30 mm on the berins, were 
rao initiating any bed movement. These two stages (I and IT) were considered © 
as representative of normal flow conditions in the channel. Boundary shear 
distribution in the unconstricted channel for the two- “Stage ¢ conditions are presented 
in Fig. 5. It was noted that the segmental overbank flow was 18% of the total 
for Stage I, and 26% for Stage. 
- Constrictions of 20%, 30%, and 40% in ‘terms of overall channel width were 
a introduced, symmetrically, from both sides of the channel (the constriction units 
_ were 4 mm thick rigid vertical walls projecting, laterally, into the flow field). 
- This type o of constriction — the most general case of bridge crossing 
TABLE 3.—Test Results of Flow Constrictions on Natural Channels. 


Nondimensional 
Maximum | Shear Maxima 
Velocity Upon | | Depth, in 
Constricting _ centimeters 
per second 


_ tracta, 
in centi- 
meters Over 

| 
| per | central 
= 


‘second channel 


(10) 
(a) Series I. _ Normal width B = 30 cm; Depth of flow over berm :- 


flow 1.5 cm, Over center 
i of channel = 4.9 cm; Mean velocity over berm = 0.10 m/sec; Over yet tral channel 


_ (b) Series I]. Normal with B = 30 cm; Depth of flow over berm = 3.0 cm; Over center | 
of channel = 6.4 cm; Mean velocity over berm = 0.14 m/sec; Over central channel 


| 
on alluvial channels with adjoining flood plains, where the uninterrupted portion 
of the waterway is, generally, marginally in excess of the central channel width. 
_ The velocity redistributions, immediately upon introduction of the constriction | 
units, as well as under equilibrium scour conditions, were measured and from A 


4 


‘The measured redistribution of the boundary in the region 
followed, in general, me pattern a.ready observed for end-dump constrictions 


section compri ‘flood plains’ or 
1 
«| meters | trac- | On 
Run tion | berm | Onc 
num-| second, | ratio, | above _ ter 
ber | | m | face | channe 
q 
2 21 =| 0.30 0.18 029 ie 16.20 1.75 | 5.10 | 1.75— 
|} 18 | 040] 025 | 033 | 13.45 | 1.95 | 5.20 2.80 3.20 
whieh the Shear pattern fOr the NOW Was SUDSequcnuy 
| 
“Ra, : 


in channels of u depth. At sattion ratios ©. 30° and 0. 40) 
pronounced high shear zones were observed against the constrietion face, whereas al 
the shear increase was nominal at 20% contraction. This can be explained. 
For the Stage II case, the overbank flow that was redirected towards the midstream 
zone was 18% and 22% of the total flow at 30% and 40% contraction ratios, 
respectively, while it was only 9% for the 20% contraction. The deflection — 
of the overbank flow into the midstream was thus substantial at higher contraction - 
Tatios; consequently, the live stream profile against the constriction face seome 
‘more curvilinear with the flow accelerating faster. 5 ma 
_ The significant feature is that the neastienddelonel shear maxima for the 40% _ 
‘contraction attained values of 2.88 and 3.40 for the two stages of flow, whereas 
the corresponding longstream velocity increase was only 1.40 and 1.50, respec- 
_ tively. Another pertinent feature was that the high shear zone conformed to 
- the path of the high velocity filament from the backwater section to the vena | 
_ contracta section. Scour holes developed adjacent to the downstream face of 


immediately following introduction of the constriction. Maximum ¢ depth « of scour . 
_ was observed to be approximately 28% of the local depth. (The observed data 
velocity and shear distribution are shown in Table 3.) 


- _ These observations of flow pattern and shear distribution through overbank 
constrictions can n only be considered as $ qualitative. Nevertheless, they do: confirm 


of the main channel flow, and the former is suddenly confined, by abrapt side 
constrictions, to the main channel zone, it is obvious that a substantial increase 
in boundary shear stress eres occur in the main channel which will result in 


- ‘The flow characteristics, in general, and the boundary shear distribution, in 
particular, ‘through end-dump constrictions and overbank flow constrictions of 
‘natural: channels in alluvium, been presented. The following 
conclusions can be drawn for flow through channel constrictions: aoe. iby 
 * The velocity distribution through a constriction is highly nonuniform; the _ 
_ flow accelerates from the maximum backwater section upstream to the vena 
contracta and decelerates gradually thereafter to that section downstream where 
expanding flow field extends, finally, over the entire width of channel. 
maximum longitudinal velocity filament, whose path i is highly y curvilinear, 
clings to the constriction face and remains close to the livestream boundary _ 


2. thei increase, both in mean flow velocity anew the gap and local velocity — 


corresponding increase 


3. Because of the curvilinear and accelerating flow field at the channel 
constriction boundary, shear maxima, f,,,,, 0n the upstream face of the constric- 


tion, and on the channel bed in th 1 the e region of the vena contracta, may Teach’ 


a 
increase for constricted flow situations, which also is the case for flow 


values many times higher than the corresponding uniform bed shear ae 


=,, in the unconstricted reach. Furthermore, experimental evidence suggests 
that the shear maxima increase exponentially with increase in the contraction 
_ As is generally the case for flow in channel bends, the locus of boundary _ 
shear stress maxima through the constriction follows, closely, the path of the re. 
-* For a given contraction ratio, the stable size of closure material required _ 
to form an end-dump dam can be evaluated satisfactorily using a shear stress an 
analysis approach. Furthermore, this approach can quantify the extent of related | ‘ 
a local scouring, and can assist also in establishing criteria for any bed protection — 
scheme that might be considered. Nondimensional curves of shear maxima versus - 
4 contraction ratio have Seen | developed with this puspess in n mind. The conventional 


the case of curvilinses: ond accelerating flow fields. ear 
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The following symbols a are used in this paper: — 
width of approach channel: 
width of opening at dump line; — 
width of livestream at vena contracta 


acceleration due to gravity; _ 
= depth of normal flow; 


= length of contraction from dumpline to vena 


contracton ratio = (B — b)/B; 


bed and water surface slope f for normal flow; 
of tow at at depths y, and y, above boundary; 
approach ‘channel r mean 
coordinates along flow direction for contracting and 


flow depths above bed “at which velocity us u, are 
measured; 


coordinates to flow direction for contracting 


4 
and expanding (Fig. 5); 
unit weight of water; 


pull 
kinematic viscosity; at 


specific gravity of bed material and water; 


geometric standard deviation of bed material; We sree 
local intensity of shear stress on constriction region, af 
_ average intensity of shear stress in unconstricted channel; — 

on ne = maximum intensity of shear stress on boundary, on 
channel bed, downstream of constriction, on face of 

die = overall angle of contraction, local angle of contraction 


Si. 


‘ pot = overall angle of expansion, local angle of expansion (Fig. 7 
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si hod MAPPING F FOR CHANNEL 


UNCTION Fiow | 


3 : The analysis oft flow behavior at channel junctions where two eggs combine 
to form a third is a problem which, so far, has received little attention in 
situations and their analysis is of importance engineers e 
in the design of irrigation systems, storm sewers, power installations, and allied _ 
projects. The complex flow pattern which results from the meeting of the two 
streams makes a complete analysis of the problem impossible. However, by 
a adopting conformal mapping techniques, a two-dimensional approximation of _ 


d the fi flow pattern at rectangular channel junctions may be carried out. This method — 
: of analysis, though it involves oversimplifying assumptions, is quite useful. 


Similar analytical solutions already have been obtained, with the help of conformal 
mapping, for other problems in the field of open channel hydraulics, such “ 


flow past a hump or a step in the bed of a channel, dividing flow at channel — 


‘The conformal ‘mapping approach was first adopted for the problem of | channel 
junctions with combining flow by Greated (2). By using , the Schwarz-Christoffel 
_ transformation, Greated developed expressions for locating the stagnation point 
; (where the dividing streamline meets the boundary) for the particular case of P 
: rectangular channel junction having a junction angle a equal to 90° and all 
_ the channels of identical bed-width. Further, by using the free-streamline concept, 
. = also obtained a streamline pattern for a 90° rectangular channel junction 


_ The method proposed by Greated is not applicable to junction angles other 
g 90°. It is also not applicable for channels of different bed-widths. Thus 3 
_ it lacks a generality of approach. A more generalized analysis of the sections ; 


4 ‘Reader, Dept. of Civ. Engrg., Malaviya Regional Engrg. Coll., Jaipur, India. 4 =o 
Lect., Dept. of Maths., Malaviya Regional Engrg. Coll., Jaipur, India. we je 
Prof. Dept. of Civ. Engrg., Malaviya Regional Engrg. Coll., Jaipur, India. oll 

~ Note. Discussion open until December 1, 1982. To extend the closing date one month, ~ 
a written request must be filed with the Manager of Technical and Professional Publications, — 
ASCE. Manuscript was submitted for review for possible publication on September , 
1980. This paper is part of the Journal of the Hydraulics Division, eer 4 of the 
_ American Society of Civil Engineers, ©ASCE, Vol. 
ISSN 0044- 796K /81 /O012- 
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“any junction angle and different bed-widths. 


At dimensional junction of rectangular channels with bed 
. oma as shown in Fig. l(a). The main assumptions in this method all > 


fluid is viscous effects are absent. sent. 


3. The loss of energy is absent 0 or, if present, it takes place downstream 


“such an n analysis does j give a good insight into the flow pattern at channel junc- 5 


= 


- Te two distinct cases which need to be considered are: 


on the method of mapping for both these cases 


using ‘the theorem of Schwarz-Christoffel, the various of 
channel junction, as represented in the z-plane [Fig. 1(a)], may be transformed — 
3 the real axis of a (-plane [Fig. 1(b)], and the points within the region of — 
flow between the boundaries are transformed to the corresponding points in _ 


the upper half of the t plane . Referring to —— the expression for the 


But, as m may ~ verified by placing ‘the following - fractions over a common 


presented in this naner Iti bcable to rectangular channel junctions having 
| 
\ 
+ 
| 
q = 
| a With F as origin in the z-plane (Fig. 1(a)] and using the standard results 


CHANNEL JUNCTION ‘FLOW 


2 


a4 


of in integration given by Cc ‘the integration 2 and 


in which bed-widths of the channels 1, 2, and 3, respectively. 
K = K, + and separating the real and imaginary quantities, — 


weve rier foi wofl} woll cosa 
+6, 6 ,) 60s a] = 


and (b,-b,)=K 


wt 
La +ayl+b) (+ aXb-a) 


COseC a 
be Poster, Eq. 2 separately in the limits a ond a + for 


point D, and b — « and b + « for the point E in (-plane a ‘which € = infinitesimal — 
distance on Lemead side of points D and E on real axis of {plane, such that 


__ FIG. 1.—Transformation Layout for Locating Stagnation Point = 
| 
— 
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( the parts only, w we e obtain 


Eqs. and by eliminating K,, one 


+ a) (1+ by 
have 
6 and 7 into Eq. 4 yields 


Dividing Eq. 9 by $ yields 


_ The preceding expressions, which correspond to the most general case, may 


be used for evaluating the unknown quantities, namely, a, b, K,, and K, for — 
ithe known values of b,, b,, b,;, and a. For a particular case of the channel E] 


_ junction with b, = b,, Eqs. 8-10 were programmed for solution on an IBM 
1130. The values of a, b, and K obtained for a equal to 30°, 60°, and 90°, 
b,/b, equal to 0.5 and 1.0, are given in Table | sections 
» As ‘thown in Fig. 1(b), if ‘in the {-plane two sources, one each at points 
C and D of strength V,b _/m and V,b ,/™, respectively, : and a sink at point 
E of strength (V,b, + V,b,)/m are considered, then the 
in the (-plane takes the form a. > 


in which V, and V, = the mean velocities of flow in n channels 1 and 2, so, . 
Differentiating both sides of Eq. 11 yields 4 


from w which the express for the at any is s obtained as 


in which V, and V, = the velocity components in the x and y directions, _ 


1 


4 


| 
q 
| 
| 
| 
| 
J &g 


4 


TABLE 1 Parameters for Location of Stagnation Point 


— 


TTT 
| 


‘CHANNEL JUNCTION FLOW 


by 


| 


~0.63789 
~0.32791 
~0.05956 
0.17503 
0.38185 
0.56555 
0.72981 
0.87755 
~0.78942 
—0.59008 
~0.40110 
0.22171 
—0.05118 
0.26578 
0.41332 
0.55423 


Distance per Unit Width 
of the Main Channel, of 


Stagnation Point 
from A Along 


ao | 


0.09911 
0.02266 


| 0.25758 x 107° 


0.07117 

0.01903 

0.08849 x 10~* 
— ¢ 


0.01453 
0.06934 
0.17110 
0.37132 


point be the 0.30 | —0.03198 
| | 62.s120 |125.270 |0.10 | 0.79008 0.50532 | 
13250] 3.6691 | 4.9785 ]0.10 = 
6 | 1.00 0.68894 | 6.6786 | 12.9680 | 0. 0.82832 | 
Ye 12.9680 | 0.10 mt) — 
1.6404] 2.1250]0.10] | 01 
| 1.00 | 1.00 | 0.38196 | 2.618 0.90| 0.54270 
.00 | 1.00 | 0.38196 | 2.6180 0.10 | 0.78940 0.50854 | 
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0.70} 0.09257 
0.80 4 19701 
0.90| 0.29315 
It is evident that somewhere along the boundary | CAD there must be a point 
of stagnation where the dividing streamline, ¥Y = V, ms meets the boundary. 
J 


_ a(1 V,b,(b — a) 
(b — a) + V,b,(1 + bd) 
Let the discharge ratio n, be defined as the 1 ratio in 


o the discharge ii in channel 3, which, , since the flow is two-dimensional, is =~ | 


“ap 
Ea. 1s into Eq. ‘14, 14, we obt 


From Eq. 16, ‘the ‘position of the stagnation point may be located in 
‘Galen. It may be seen from Eq. 16 that for any value of a a, the value of | 1 
4 corresponding to the stagnation point depends on the value of the discharge _ 
ratio, n,. In other words, as the value of n, varies, the stagnation point moves | 
along the boundary CAD from C to D and, correspondingly, the value of s 
: varies from —1 to a. Thus, when the stagnation point lies along the boundary © 
_ CA, ¢ varies from -1 ‘to 0, and when it lies along the boundary AD, C varies 7 


F from 0 to a. The stagnation po point will coincide with corner A when l= = 0, 


corresponding to which the value of Na = be ices as critical, _ 
n,_» which is obtained from = as 


| 


4 


= or 
5 
| 
~ 
{ 
— 
| 
may be integrated 
along AC and AD, whic ‘eo | 
in which -1 <{ <a. by 


In order to obtain the distance of the stagnation point from the corner A » 
along the boundary AC or AD, the integral on R.H.S. of Eq. 18 has been > 
evaluated numerically by using Simpson’s rule on an IBM 1130. The values 

ra of the distance of the stagnation point so obtained for the different values | ; 

n,,a, and b,/b, are givenin ‘nondimensional form in Table 1. The corresponding 
values of ¢ as given by Eq. 16, and the values of n, |. as given by Eq. 17, 

ate also gweninTablel, ‘ont 24) al SG 
Case (B) Fuow Separatine at Bounpary: Free-Streamune THEORY al 


‘e According to ‘the idealized theory described in the previous case, with no 
‘Separation at any point on the boundary, the ge sary at the downstream comer — 


corner has to separate from the seendiins at the corner, thereby developing 
a finite tee a oh at the corner. In order to account for this actual Legal 


FREE- E-STREAMLINE 


CHANNEL 3 
Z-PLANE 
— 


@-PLANE 


fe) PLANE 


| streamline, which separates the fluid in motion from the fluid at rest, is called - 
a 
i 
FIG. 2.—Transformation Layout for Locating Free-Streamline when Stagnation Point 
SCoincides with Corner 4 OR ho ¢ al 


DECEMBER 
a free-streamline. The velocity as ; well a as the | pressure remains ‘constant along = 
a free-streamline. Accordingly, as shown in Fig. 2(a), the streamline BE, starting 
from the downstream corner B of the junction, may be considered to be a 
: free-streamline. The fluid between the streamline BE and the channel boundary a 
_ BE is assumed to be stagnant. Consequently, the flow on the downstream of = 
the junction is confined between the free-streamline BE and the channel boundary _ 
_DE. In the study of the flow pattern at the junction, it is, therefore, essential z 
to locate the free-streamline starting from the downstream corner B of the — 
junction. In the following Paragraphs the method for the location of the 
_ The analytic function, Q = = {[{In (V;/v)] + 16} = In [V, (dz/dw)) (in 
y = the modulus of V. — iV,), transforms the boundaries of the z-plane o 
= Q-plane, in which ‘the bodogrephs or the velocity — assume one of © a 


2 along the boundary The free- streamline f ame cases 


be 


FREE-STREAMLINE 


_ 


AG. 3.—Transformation 
lo 


tion 

q S Lies in in Channel 1 segduneamaidht AD 

g 


point S being coincident with corner A; (2) the stagnation point S lying : 
— 
4 
— 
| 
| 
Streamline when Stagnation Point 
| 
| 


s Stagnation Point sc Coinciding with Corner A: :—Using the Schwarz- Christoffel — 
aa theorem, the region included in the semi-infinite strip of the QD -plane of Fig. 


: 2(b) can be transformed ot onto the upper half of the A- plane, a: as shown in ‘ae 
The corresponding for the transformation is 
wh 


20) 


a which C = constant of integration, which, along with K, ay be ‘evaluated 


by using the boundary conditions as follows. __ 
4 For point E: 4 = 0, Q = im; and for point B: oie 
thus C = in, and K = Eq. 20then becomes 
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the 
Since se the | flow near ihe sins streamline is of prime interest, a relation between 
‘¥ and \ must be found which maps the interior of the infinite strip in w-plane 
" a E, F, DandC shown in Figs. 2(d) and 2(e). The corresponding transformation 
| which . K’: = a complex constant. 
Integration 0 of 2B yields 


w= K’ ind) 4 + 
» may 


_ in which C’ = constant of integration which, a with K’ 


by using the boundary conditions as follows. 
i V,; thus 


_ For point F: A = 1, w = 0; and for point B: dh 1, w = 
=0,and K’ =¥ in which = the stream function of streamline 


«Eq. 24 then becomes 


which, on differentiation with respect to to yields 


Further, oe ‘into Eq. 22, we obtain 


in which \ = —£; dd = - sin’ , or arc sin \ 

By 


d taki 
a as 


sin + isin arc sin V ) 


 —_ 
| 
4 
meg 
4 nd E yields 
7 


Separating the ‘real | imaginary quantities ‘gives equations for the 


_free- -streamline in the parametric form: 


dé: 


sin 


(Bai arc sin VE 


- 


Further, introduction of an auxiliary variable t=7, 


€ preceding equations to be as 
sin 
ow, 
| 
Vs 
|. sini —arcsing 


— arc sin y 


7 


n 


2n 


permits 


dt = 


| 


_ ‘Eq. 31 represents the coordinates of the points on the free-streamline 


( 


t 

v,= 

— arc siny 


tye 


cos 


sin (=. arc sin 


san) 


(> 


Point S Lying in Channel 1 Along Boundary AD.—Using the 
the shaded region of Q-plane of Fig. 3(b) can 


Christoffel theorem, 


in 


V,b,, Eq. 31 may be expressed in ‘nondimensional form as was 


4 be transformed c onto o the > Upper half of the ) - -plane, as ‘shown in a Fig. 3(¢). TI The — 


“an 


Integration of Eq. 33 between 0 to —), ee is, from B towards E, yields 


ad end 4 


g 
| 
——— 
| | 
| —— 
| 
ad 
 Q-i(r+a)= 
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es point E: A = thus = 1/2; and Q = in. Introducing these values 
— 


> 


integration of Eq. 33 in the b - €, rer the | point 


(such that — 0), 


From Egs. 36 and 37, K 
= 
Vb(l+b) 
all Further, tear Eq. 33 between points B and C we ate! 
| 


in which \ = ee 


which m = a)/V b(l + b) 


a From Eq. 38 41, ‘for given values of n,, and b a/ by, the values 


b(1 +a) 


the parameters a and b can be obtained forthiscase. = 


_ Since the flow near the free-streamline is of prime interest, a eal between — 
w and a has to be found which — the interior of the infinite st 


points B, E, F, D and C as shown in Fig. 3(d) and 3(e). ‘The dag 


of the expression in the limits, - + « and -1 -«6 


&g 
37) 
38) 
Pao 
39) 
ration,, since b, V,/b, 


Zz 
TABLE 2. —Dimensioniess Coordinates of on with F as Origin 
‘when Stagnation Point S Coincides with Corner 


of 


 @a= 30, @, /b,) = 1.0; (b,/b,) = 1.0; Critical n, = 0.50156 (Limiting value of 


a = 60°; (b,/b,) = = Critical n, = 0.53101 value of 


neo —0.0444 0.9499 
70.1027 
-0. 1745 


—0.5046 
—0.6826 


(c) a = 90°; (6,/b,) = 1.0; (6,/b,) = 1.0; Critical mn, = 0.61803 (Limiting value of . 


~0.8356 
~1.2720 


| 

| | 
| 

| 

a 


; TABLE 3.— 3. —Dimensioniess Coordinates of Points on Free-Streamline with F as — 


when Stagnation Point Lies Along AD in 


m= 006 on, = 0.8 
(@a= 30; (b, /b,) = 1.0; (6,/b,) = 1.0 (Limiting value of (y/b, ) at n, = 0.6-0 
is 0.7660, 0.7583, 0.7424, and 0.7117, 


0.0356 0.0355 —0.0352 


2 oy 
0.4196 0.4157 

—0.6580 
9954 0.9954 0.9953 

0.8767 


0.9825 
na 
0.8413 


(b) a = 60°; (b, /b,) = ‘1.0; (b,/b,) = 1.0 (Limiting value of (y/6,) at n, = 0.6-0.9 
is 0.5161, 0.4934, 0.4628, and 0.4207, respectively). 
(x/b,)—0.0054 | 
0.0610 -0.0598 = 0.0554 
—0.1194 0.1159 | 0.1104 
0.3606] 
(y/b,) 0.9919 0.9917 


9674 
9272 
8718 

0.8020 


de 


0.9095 


(ce) a = 90°; (b,/b,) = 1.0; (, roe 1.0 (Limiting value of (y/b,) at n, = 0.7-0 9, 
is 0.2662, 0.2373, and 0.2028, respectively). 


~0.0239 = = 


| 
— 
ree ~0.4597 0.4427 9/5799 
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= it 


Toast J 
9 
for dd/dw 44 yields 


in which = —sin and = sin n cos dé. On integration, t the 
expression Yields: 


fs 


__ 364 


exp Qtand db 


exp lice +f)] = —ex 
a- arc tan tan b> 
= db 
Separating the real imaginary quantities me the 


x=-— 


: Eqs. 52 and 53 ‘represent nt the coordinates of the | points on the free- tialies. 4 
inthe zplane with Fas origin, 


Further, since V, b,V;, Eqs. 52 and 53 may be exp expressed i in nondimensional 


o 


“4 
7 
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TABLE 4.—Dimensionless Coordinates of Points on Free- Streamline with Fas 
when Point Lies AC i in Channel 


—0.0361 


0359 
“a 0. 


—0.0358 
—0.0843 - 
0.1588 
0.2671 -0. 2665 
0.4241 | -0. 4233 
-0.6637 | -0.6627 | 
—1,0927 


(x/b,)—0.0087 0. 0087 —0.0086 


0.9955 | 0.9955 | 0.9955 
0.9839 | 0.9834 0.9831 9830 0.9830 
0.9671 | 0.9654 | 0.9646 | 0.9641 | 0.9639 
0. 9410 0.9401 | 0.9396 

0.9249 «0917s | 0.9137 | 
0.9021 | 0.8898 8836 0.8806 | 0.8793 
0.8790 0.8607 | o8sis | 0.8471 | 0.8451 

wae | 08308 | 0.8181 | 0.8120 | 0.8094 


(6) a = 60°; (b,/b,) = 1.0; (b,/b,) = 1.0 (Limiting value of (y/b,) at n, = 01-05, 
= 0.6586, 0.6216, 0.5960, and 0.5785, respectively), 


 (x/b,)-0.0058 -0.0057 | -0.0056 | -0.0056 | -0.0056 
0.0267 0.0261 | -0. ~0.0252 
-0.0667 | -0. -0.0645 | -0.0638 
—0.1336 -0.1276 
ey —0.2357 2314 | -0. 0.2258 
0.3962. —0.3883 | -—0.3822_ | 0.3742 
70.6348 | -0.6250 | -0.6173. at 0.6072 
) 0.9922 0.9921 | 0. 0.9920 
0.9709 0.9704 97 9697 | 0.9695 | 
9399 | 0.9379 | 0.9365 | 0.9355 | 0.9348 
9029 0.8978 | 0.8916 | 0.888 | 
0.8631. 0.8527 8455 | 0.8404 


0.7073 | 0.6795 | 0.6 


(c) ia = 90": (b, = 1.0; (b, /b,) = 1.0 (Limiting value of (y/b,) at n, = 0.1-0.6, 
0.5443, 0.5290, 0.4824, 0.4470, 0. and 0. 3986, respectively), 


-0.0015 | -0.0014 | -0.0013 
| | —0.0102 


-0.0392 —0.0346 
0.0933. =0.0832 | 0.0832. 
—0.1673 —0.1673 
0.3408 -0.3292 0.3194 | -0.3038 | 3038 
___-0.5795 -0.5608 | -0.5482 -0.5281_ 5281 


| 
j-imiting value of (y/b,) at n, = 0.105, 
spectively). 
| -0.0086 | -0.0086 [| 
0.0841 
0.2702 0.2663 
0.4285 0.4257 0.4230 | 
| 
| 
| 
j 
—0.2979 
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S Coincides Corner 


i. Stagnation Point in ‘Channel 2 Along Boundary Ac. —Using the 


| 9.9488 | -0.9400 
09903 | 0993 
gai 0.9184 (0.9175 } 0.9167 | 0.9161 09156 
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. tions when Stagnation Point 


Schwarz-Christoffel theorem, the shaded region of Q- -plane of F 


be transformed onto the upper half of the A-plane as shown in Fig ri 


for the transformation | « 


in which A = cos” 6: dy = —2 cos sin do; and taking -2= 


integration of Eq. 56 from B 


Further, it may be own that for tk this cas case also, the parameters ‘Kia a, 


7 Again, by adopting the same same procedure as —— in section 2, the integration — 


of Eq. 56 between points D and E leds to” i 


FIG. 6.—Free-Streamlines for 30°, 60°, and 90° Junctions when Stagnation Point 
Lies in Channel 1 on on Boundary AD, Discharge R Ratio n, = 0. 
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+a 


+ a) al + 
le 
in which m = — - + b = =(1 wore from Eq. 38. 
The factor V,/V, is representative of discharge ratio n since b, V;/ b, 


From Egs. 38 and 59, for given values of a, b / bs, the values of 
_ the parameters a and b can be obtained for this case. re 
Next, for this case also, a relationship between w and ) may be found which | . 
maps the interior of the infinite strip in the w-plane onto the upper half of — 
>. the A-plane, with correspondence between the points B, E, F, D, and C as 


_ shown in Figs. 4(d) and 4(e). Further, adopting the same procedure as in section 


spy 


of To sole zai] iniog 

+t. 

ary bas 

FIG. 7. —Free-Streamlines for 30°, 60°, and 90° Junctions when Cugnaten 2 a 

Lies in Channel 2on AC, Discharge Ratio n, = 0. 
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in which arc tan tan (62) 
=... Further, since ¥, = V, b,, Eqs. 60 and 61 may be expressed in nondimensional 


x 2 i 


In order to obtain the coordinates of the on the free-streamline, 


Z integrals on R.H.S. of Eqs. 32, 54, 55, and 63 have been evaluated numerically 
_ by using Simpson’s rule on an oe The values of the coordinates obtained 
for b,/b, and b,/b, = to l and a = = to 30°, 6, and 90°, for each of the 
three cases, are given in Tables r ap 8 and 4, respectively. Further, Figs. 5- a 


An solution has been obtained, the help of mapping, 
for the problem of channel junction with combining flow. The analysis indicates  __ 
that for smaller values of the discharge ratio n, (corresponding to which the 
discharge in channel 2 is less than that in channel 1), the stagnation point is 

located along the boundary AC of channel 2. With the increase in the Value 
of n,, the stagnation point approaches the upstream corner A of the junction. 
‘For es each junction there exists a unique value of n, termed as n,_, corresponding — 
to we the stagnation point coincides with the upstream corner A. The value 
n, 


pre lies along the boundary AD of channel 1. ¥ 

- For different locations of the stagnation points in nae ‘channel junctions, the 

_ solutions of the equations for the free-streamline have been presented in 


dimensionless form. The zone of the stagnant fluid between the free-streamline , 


Ps Such an analytical approach i is helpful in n developing insight about ‘the flow 
behavior at the junction. It also presents quantitative conclusions which could - 

f obe compared with model studies in order to understand the effect of the 


q 
| 
| 
cases for the three values of the junction angle. 
i 
| 
7 and the channel boundary BE increases with the increase in the junction angle _ i 
assumptions made in the analytical Gevelopmemt. Pr 
| 
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ar i ts paper 
e following symbols are use in this pape: 
= distances of points on on real axes of C and athicina 
= _ bed-widths of rectangular channels 1,2,and3; 
= constant which could be complex number = K, + iK,; 
4 = discharge ratio= 
= complex variable = [log (V,/v) + ie seiaialiilaih 
formation plane for boundaries of channel 
002.0 = discharges in channels 1, 2, and 3 (Q, + Q, = Q;); 
= mean velocities of flow in channels 1, 2, and 3; - 
velocity components in x and y directions; 
: filament velocity = modulus of V, — iV,; 
= complex variable representing rectangular flow field; 


= + infiehocmel distance on either side of points on = a axis of 


= complex variable representing transformation plane for 
aries of channel junction; 


aries of channel junction; — 
auxiliary variable; and 
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Journal or on any paper presented at any Specialty Conference or other meeting, the Proceedings 
ot i who has significant comments or questions regarding the content of the paper/technical _ . 
note. Discussions are accepted for a period of 4 months following the date of publication 
_of a paper/technical note and they should be sent to the Manager of Technical and Professional 
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Loss From SIDEARMS oF Coouinc Lakes" LAKES 
teats och Kennedy,°M. ASCE 
‘The writers are gratified by Jain’s (11) use of ‘their r numerical re results a and 
‘the verification he obtained of their proposed correlations for surface heat loss 
f from cooling-lake sidearms. The co contention by Jai Jain 3 and Yildirim that the boundary 
condition, (T),-0 = T,, is incorrect requires some > clarification and comment. 
One of the principal advantages of the Karman-Pohlhausen approach to numerical _ 
analysis of equations of the boundary-layer type, which was used by the writers, § 
is that it is unnecessary to know precisely the vertical distributions of the variables _ 
of interest (temperature and velocity in the present instance), but only to reproduce 
eae essential features. The © experiments 's conducted | by th the writers were used 
as a guide in selecting an _ acceptable profile shape, shown in Fig. 2, for the 
vertical temperature distribution, which then was used in the Karman-Pohlhausen 
analysis of the equations. Although the experimental results exhibited some > 
smaii longitudinal variation in the tepenatuse at y = 0, the experiments = 


that the condition (T),.0 = = T, is approximately correct, and no more than 
ti was claimed by the writers. An ordering analysis of the continuity equation | 
and the equations of motion clearly shows that a representative vertical velocity 
_Vis of order 8 U, where 8 = D/L and U is a representative horizontal vilodiy 
= (10). The authors imposed the restriction L >> D or 8 << | in their analysis, 
z which case the vertical velocity is much smaller than the horizontal velocity. — 
: Now in the main reach of the sidearm, where the flow is nearly parallel and — 
_ for the case of A, = 2 it is apparent from Fig. 3 that the temperature of — 
the fluid convected from the inflow current downward to the outflow cerrent 
is nearly equal to 7, at the point y = A, in which the vertical velocity attains 
_ its maximum but is stiii much smaller than the characteristic velocity U.— 
~~ vertical turbulent diffusion downward into the outflow current Jf 
from the inflow current is inhibited by the strong vertical stratification in the — 
main reach of the sidearm. The vertical heat transfer is dominated by the heat 
_ loss at the water surface, and the experimental temperature profiles shown | 
_ im Fig. 4reflect this overriding influence, 
The equations developed by Jain and Yildirim can also be used to develop 
a line of reasoning that differs from theirs, as follows. Eq. 63 shows that (e" wy ; 
=, and the condition (8),-o=0 in Eq. 64 would require that (@”),.9 = 
as well. However, Eq. 62 shows that the denominator of Eq. 65 then is also 
: “zero, so (f”),.o, and therefore (du/dy),.o, may have finite, determinate values. 
_ Viewed more simply, Eq. 62 shows that (6), _. = 0 is compatible if k,, approaches — 
the constant fluid viscosity near the wall and f is no more than a third-degree 
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DISCUSSION 
— 


a 


or if”), ois otherwise sero. In of In oth In other er words, th the dis discussers’ relations 
are not incompatible with (®),-0 = =0. 
a Finally, it should be ‘emphasized that the writers’ results were not generated ; 
i for use with sidearms having lengths greater than or equal to equilibrium in ; 
_ length, but were intended also for sidearms with L < L,, in which case the 
~ equilibrium length results meery provide convenient length and heat-loss scales 
for sidearms of length L < L,, as shown by Fig. 8. In this case, the vertical 
velocities become large in the closed end reach as a result of the fixed boundary. 
The inflow current is abruptly « deflected and returns | along the bottom as the 
outflow. It is in the relatively short closed-end reach that these large a ; 
currents establish the value of T, which is then advected along the sidearm 
bottom by the outflow current. This was confirmed by the writers’ experimental 
"results w which showed that (0),.0 = = 0 is a good approximation for | L< <L, emus) = 
_ Errata.—The following « correction should be made to the discussion by Jain 
Page 523, Ref. 11: Should read ‘ “Journal al of the Energy Division,’ "instead 


RANSVERSE MIxinG Tests IN NATURAL neal 


Teor 
Closure by Spyridon Beltaos’ ak 
The writer thanks Fischer for his interest in the paper. The discussion focuses _- 
« two items: (1) Concern about the accuracy of the transverse mixing coefficients 
reported in the paper; and (2) a theoretical explanation offered by the discusser_ 
_ of why H rather than R gives somewhat more consistent results when nondimen- 
Under item (1), the discusser voiced several objections. to the approach used 
in the paper but has not explicitly stated whether he considers these objections 
sufficiently serious to Tender wed ‘Teported coefficients from Gs 


a Before dealing with items (1) and (2) in detail, the writer would like to establish sf 
_a few basic points. Firstly, it must be recognized that €,,, €,, and E, are not 


{October, 1980, by Spyridon Beltaos (Proc. Paper 15731), 
> Research Scientist, Environmental Hydr. Section, Hydr. Div., National Water — a 


Canada Centre for Inland Waters, Ontario; Research Officer, 


; B If the contrary were true, the discusser would not have attempted to theoretically atl 
| 
| 


at a point. These coefficients can only be estimated as ‘those values which — 
_ optimize agreement between measured concentrations and those predicted via 
an analytical or numerical solution of a differential equation, such as Eqs. 1 
and 2. (Note that these equations already involve several assumptions and — 
_ approximations; with these assumptions, they can be derived from a more general 
_ 3-dimensional equation by a depth-averaging operation—see also Ref. 10.) The _ 
- value of a coefficient will thus depend on the computational algorithm 
used if a numerical solution is pursued or on the simplifying assumptions made — 
in deriving an analytical solution. All coefficients so determined are ‘‘accurate’’ 
applied in conjunction with the corresponding method and “inaccurate” if 


: g thus seems to the writer that the question at hand is not so much one of 
-“‘accuracy”’ but, rather, of ‘‘acceptability.”’ For example, if a numerical algorithm | 
ie an analytical solution involve assumptions that contradict established knowl- _ 
ede or experience and if such contradictions have a significant effect on the 

solution, they should be discarded even if they occasionally succeed in producing 
4 *‘good’’ predictions by manipulating the coefficients. Those methods that are 
free of such shortcomings are acceptable while the final choice will depend 
on proven performance under varied conditions and on ease of application. 
In the subject of transverse mixing, a final choice of this type has yet to be 
- made. (One of the aims of the paper was to contribute toward making a final © 
choice and not to advocate what the final choice should be. = a a ; 


: by the discusser. On the analytical side, two basic ‘methods have been developed: 
(1) The constant velocity-depth-E,-method (see also Ref. 27) that leads to Eq. 
21; and (2) the streamtube method that has been used i in the paper. The former 
method clearly contradicts” experience | and is not popular at present. The 
streamtube method also contradicts experience (h’u,€,, varies across a stream) 
but this contradiction does not seriously affect the solution, as shown in Ref. 
17 and as additional evidence will show later in this closure. (As one illustration | 
of the latter method’ S superiority over the former, one may compare the scatter 
of the paper’s versus graphs that of versus x graphs give 
in Previous literature, such as Ref. 16, €. g.: ; other illustrations may be found 


any obvious compromises in logic, other than those involved in deriving the 

differential equation that is being solved; however, they are far from being 

solutions because of the following major shortcomings: 

_ . The region of computation must be discretized, i.e., an) equation that is 

derived for an infinitesimally small stream volume is applied to elements of 
finite and sometimes large volume. This inevitably introduces errors and leads 

to a dependency of the evaluated mixing coefficient on the particular algorithm = 
car) sbitingem Yo ila at well bas 

thy . Ifa significant advantage is to be gained from a numerical scheme’s ability 

to deal with variable A and u a> then these quantities should be known throughout — 

7 _ the region of computation. In practice, however, transverse variations of h— 

and u,are only known at a few cross sections an are considered ° aeeennanen ig 
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“(This i is caused by: limitations: e.g., to obtain a adequate depth 


information at, say, 15 points across a stream, a crew of three experienced — 

- persons using a boat for access require two to three hours; this is not counting __ 
“unproductive” work time, i.e., travel, boat loading, and launching etc.) Between 

_ “metered’’ cross sections, A and u, are generated by some kind of interpolation 

. technique; however, the detailed river geometry does not obey neat mathematical 


a and thus one can only hope that the errors introduced by such interpolations 


tions were assumed: For the reach from the injection point to measuring 7 : 
station 5 (at the end of the bend to the left), the velocity distribution — a 
measured upstream from the injection point was assumed to > apply aie 
all sections. At mile 645 the measured velocity distribution was used, aan 
# except that the 70-ft width of shallow water near the left bank was omitted 


‘same min and depths were used as at mile 645, ‘except that distances = 
- measured from the right bank were assumed to be measured from the 
left bank. At station 10 (in the middle of a bend to the right), the velocity _ 

_ distribution was assumed to be the mirror image of the velocity distribution © 
Le upstream from the injection point, except that at this point the width _ 
. of the stream is only 510 ft, so that all transverse distances were reduced ud 
by the fi factor 510/600 and all depths were increased by the inverse of = 
a same factor. For cross sections between the above sections, depths, Je 
widths, and velocities were interpolated linearly from the nearest assumed -_ 


placed the w writer. .) Where so many unverified and subjective assump 


are considered acceptable, it is only a very small additional compromise that 
one has to make in order to use the streamtube method (which on many occasions 
has been found to perform adequately). This compromise is done in exchange © 
_ of the many advantages of analytical solutions, i.e., incomparably easier applica- 

- tion; plus an ability to provide the designer with a complete appreciation of 
the effects of various pertinent factors without having to perform a separate 

_ numerical computation for each possible combination of such factors. (It is 
noted here that, usually, more than two cross sections are metered for any 

_ particular test, the number varying up to about 10. This, however, only moderates 
the basic deficiency identified in the foregoing: The region of unknown stream 
geometry and flow properties is still orders of magnitude greater than that of ii 

With this background, the discusser’s comments are now considered in detail: 


The discusser states ‘‘In earlier work (16) | the author has o- weed te 


= E,? It is s also reasonable to ‘remark that, where such gross approximations 


— 
_ quotation, taken from the discusser’s suggested example (16) for a 33,000 ft — 
long reach of the Missouri River, 
f = a point halfway from section 8 to section 9, the velocity distribution was = #§$*‘&f 


stream tube method to ascertain the value of E., but did not 1 require he assumption | 
that h’u,€,, is constant.”” (By the term “‘the author’ the discusser must be 
7 q referring to himself and not to the writer because Ref. 16 is authored by Yotsukura, 4 
Fischer, and Sayre.) The writer would like to clarify that his use of the term > 
“‘streamtube method”’ implies the mathematical transformation that results when 
is substituted by and the analytical solution that is afforded 
when the assumption h’u,€,, 
; 2. The discusser states, as a a matter of fact, that “e ~ uh ~ h?”?.” (Note 
that the “proportional to”’ sign, «, is more than the “‘approxi- 
mately equal”’ sign, ~.) The only fact in this case is that putting «,, h « 
7 is only one of a multitude of possible hypotheses concerning the variation of — 
ne across a stream; the basic cause of this state of affairs is the previously — 
‘mentioned fact that <, cannot be directly measured at a point. For two-dimensional — 
flow (h = const = H; u, = const = V; and u, = const = V,), it has been : 
established that «,, « u, h = V, H. This finding does not necessarily a 
» the discusser’s extrapolation to channels in which h, u,, and u, vary with 
_z. The only logical conclusion that results from the two-dimensional flow al flow findings 


The. depth-averaged transverse mixing coefficient an an function A 
_of z (or of qg); in the case of two-dimensional flow, this function becomes 
constant and this constant is proportional to H. 


discusser’s assumed functional relationship is ‘‘possible’’ because it satisfies 


7 _ the aforementioned requirement, but it is by no means the only possible — 
assumption. In strict logical terms, there exists an infinite number of such | 
“‘possible’’ relationships (e.g., put €,,« V, Hf(h/H) withf(h/H) being any 
“‘well behaved”’ function of h/ H). Of this multitude, a few “‘plausible’’ functions 
i have been suggested and studied in the past decade. For example, Ref. . 10 


does not vary across the stream us be 
The following quotation from Ref. 10 is particularly pertinent at this eta 


- It is recognized that there are dangers in using equations like Eq a 
and 19 (writer’s note: means equations like the foregoing) because their > 
_ derivation makes a local application of relationships which were obtained _— 
_ for cross-sectional average values. In presenting these equations, it is a 
not felt that they necessarily represent the exact variation of K, (writer’ _ 
on note: K, is equivalent to €,,). Rather, they are represented as possible 
_ Though «,, cannot be measured directly, considerable indirect evidence concern- 
_ a ing the variation of €,, is available in the literature, as described next. mn nly 7 
the soning of the equation: old , 
= aie i) Holly and Simons @ 3) tudied the comparative merits of various possible — 


|. 
| | 


al 
dependencies of in 


rather extreme cross-sectional quemetry wes deliberately sO as 
to produce strong lateral variations of h, uy, and u, 
Three possibilities for were considered, i.e.: : (i) = const; 
(2)€,, © uh = variable; and €,, = variable: s postulate. 
: “ pe The governing differential equation was solved numerically for each of | Pd 


a. these cases and the predictions were compared against the tent results. a, 
The outcome was conclusively in favor of the constant t assumption 
ae 


n 
as may be seen in the abstract of that paper: — 


- x Turbulent dispersion experiments in a triangular laboratory channel suggest 
that the transverse diffusivity does not vary with transverse position as 
J : previously suggested. The invariance of the diffusivity may reflect an 
interaction between mixing due to bed shear and due to transverse shear an 


_ the relative contributions of which cannot be determined. “ee 


A similar study with trapezoidal — ayee 


of numerical simulations using different for 
variation of e e, (writer’s note: e, is equivalent to €,,) across a a channel 
¥ showed that none of the assumptions appeared to be definitely superior. : 

For the sake of convenience and because of the availability of published 
data, it is that e, should be considered across the 


Writer’ s note: this ‘finding does not necessarily ‘moderate the : strong 


evidence presented by Holly and Simons (23) in favor of the constant — 
Sea assumption; it is likely that the superiority of this assumption was — 


"channels where €.4 should be mostly diffusive. (though effects 
"may well be present due to pooepenty currents that occur even in —_— 


(iii) The effects of meanders on transverse mixing were ina rectangular 
_ ‘meandering flume by Chang (22). Of particular interest is Chang’s Fig. 
43 where the dispersive part of E,, (defined by €,, in Eq. 11 of  . 
a paper) was plotted versus z in a nondimensional form. The data points 
Z. a considerable scatter but no consistent variation. The coefficient 

p was calculated from rary according to the well- -known hell 

hypothesis: €,, =. —(w'C Ya / (Cy in which 
= local lateral velocity | component; C’=C- 

“ss concentration. (For steady-state conditions, C and \ w may be considered 
time-averaged values; for unsteady conditions, they as 


= 
| 
s by Lau and Krish- ase 
| 
HOt aS Cvident in the Case OI Kel. Zo Decause c Channel geometry was 
¢lose to that of a two-dimensional channel.) 
_ associated helical flows have not been considered so far. The following _ 
be 
a 1 


 eross sections do not ‘adequately stream conditions. _ 
ae better reproduce typical river geometries, they used a meandering 
: _ flume with ‘‘equilibrium” bed (i. e., the bed was initially mobile and : 
final configuration was formed by the flow itself). Again, several — 
variations of E,, were examined, including the assumption E,, = const. 


lire It was found that the latter yielded the most consistent results and = 
it was concluded that, for all practical purposes, the transverse mining 7" 


/? _ coefficient can be assumed to be constant across the stream. __ woul 
hey 


(vy) Recently, Lau and Krishnappan (26) performed a comprehensive cate: 
7 son of various assumptions for E., using the results of a field test as 
7 the basis of evaluation. The governing differential equation for meandering 
flow was first transformed to x and qg coordinates and then solved by — q 
Peer a numerical scheme. The diffusion factor D, (see Eq. 10 of the paper) — 
ss ' Was assumed to take five ogee forms; of these, two are particularly 


isthes pertinent at viz 


E cons 


re Eq. 28, E -¢ iS assumed to vary in proportion to A rather than u,h 
= as hypothesized by the discusser but this variation is much closer ‘to 


 discusser’ s s hypothesis than to the constant - assumption. ‘if the 


about Eq. 27. However, Ref. 26 shows that there isno difference i in the respective 
prediction abilities of Eqs. 27 and 28.00 
_ The above evidence is indirect and thus does not conclusively ‘ ‘prove”’ that — 

4 the transverse mixing coefficient is indeed constant across a stream. However, — 
this evidence strongly suggests that even if E,,varies across a stream, its variation ia 
must be mild and nowhere near the h’/? ‘proportionality that the discusser has 
stated as a fact. It thus makes good engineering sense to use a constant E., 

g 7 since this assumption not only “‘works’’ in practice but also happens to nol 

The preceding discussion renders practically irrelevant the discusser’s 

= objections concerning the definition of and the validity of Eqs. 6, 7, and 


Ea. Academically speaking, even if €,, is assumed to vary with z (or with q), 
4 the resulting correlation term can be absorbed in the definition of e, which — 
_ will then differ from €,,. Eq. 8 is indeed incorrect, as may be seen by a check 4 
of the dimensions of the two sides. The correct form is given in the “Errata,” 


4. The quantity h* 
~ though not as h* but more likely as es. The discusser i is referred | to the text 
- following Eq. 4 where it is stated that, even though h’u, €,, does vary with 
q, this variation does not greatly influence the solution of Eq. 4. Such fortunate 
pos ssa / are not uncommon in engineering and it is a good practice a 
exploit them. ‘Since Yotsukura and Cobb’s (17) first evidence in favor of this 
assumption, additional evidence has been published (14,26). Briefly, this informa 
tion shows again that the transverse variation of h’u, €,, does not significantly 
affect the solution of the governing differential equation; more serious is the © 
longitudinal variation of h’ Me €.4 which we have no way of even guessing at — 


| 
| 
{ 


- ~ uncertainty that must be “allowed for in pertinent designs. The methods — 
in the paper provide estimates of reach-averaged values of E,; investigation 
of the longitudinal variation of E, was not attempted. (An illustration of the 
latter is however given implicitly in Figs. 6, 8, 10, and 11: e.g., join individual — 
data points by a continuous curve and consider how much the slope of this 


owe cent & Cold Lake 


beans (14. —Reach of Tests ; 3 and 4 ain 3; Sites 1 “12; oe 4: Sites A- “ad 


= FIG. 15. — of Tests 5 and 6 (Test 5: Sampling at Sites 1-7; Test 6: Slade 

Sites A-E; Primes Denote Additional Metered Sites) 

curve & at any one ‘one value of the abscissa deviates from the slope « of the straight 
_ §. There is one instance where the streamtube method may not — ~~ 
satisfactory results, i.e., where a plume is mostly within a region very near ct 
a river bank where (usually) maximum lateral gradients of h occur. Such conditions Z 
generally prevail for short distances of side sources. If accurate 


1744 
iy ARE 


- 


Me are allowed to vary as s measured. 


 * Maps of the test reaches for tests 3-4 and 5-6 are given in 2 Pik. 14 and =. 
. At the same time, it is pointed out that Ref. 2 may be obtained simply — 
_ by addressing an ordinary reprint request to: Transportation and Surface Water _ 


~ Engineering Department, Alberta Research Council, 4th Floor, Terrace Plaza : 


Office Tower, 4445 Calgary Trail South, Edmonton, Alta., Canada, T6H SR7. 
: % 7. The discusser’s explanation concerning the use of H in place of R would 


7 have been satisfactory had the experiments been performed in straight prismatic — - 


channels. In such instances, one may reasonably expect that E, reflects mostly _ 
‘(though not entirely) diffusion caused by the random motions of turbulence, 
whereby the discusser’ s arguments about length and velocity scales follow 
logically. However, the writer’s experiments were performed in natural streams 
where, as is well known and as pointed out in the paper, a large part of E, 7 
ae dispersive effects caused by helical motions. The discusser himself was 
_ one of the first investigators to identify and study this effect (8). = 
‘transverse dispersion of this kind has little to do with turbuler 
does not consider the discusser’s explanation satisfactory. 
_ 8. Finally, the writer objects to the discusser’s statement “The author’s use 
of R implies that the effect of the ice cover is to halve the length scale of _ 
‘the turbulence.”’ The suggested implication is purely the result of the discusser’s 
; own interpretation. The writer listed E,/RV, (= K,) in Table | simply because 
this enables” a direct comparison with the finding of Ref. 6 where K2, was 
_ about the same for both open- -water and ice- ~covered conditions. The writer ; 


sentences, quoted from mp. 1622: 


However, tenpineien of T Table 1 shows that values of K, are larger under 
é ice cover than with open water flow by as much as two and one- — 
‘ times. It is probably premature to attempt an explanation of this finding, 77 
7 however, it is believed that much may be learned by studying the helical | 
~ bend flow under an ice cover (see also Ref. 8). 


22. Chang, Y., ‘‘Lateral Mixing in Meandering Channels,’ ° thesis | presented to the ‘University 
a lowa, at Iowa City, Iowa, in 1971, in partial fulfillment of the requirements for 
the degree of ‘Doctor of Philosophy. 
- 2B. Holly, F. M., Jr., and Simons, D. B., “Transverse Mixing of Neutrally-Buoyant 
Tracers in Nonrectangular Channels,” of the XVIth Congress of the 


_ International Association for Hydraulic Research, 1975, Vol. 3, pp. 467-474. | wea 
24. Krishnappan, B. G., and Lau, Y. L., “‘Transverse Mixing in Meandering Channels 
a with Varying Bottom Topography,’ °” Journal of Hydraulic Research, International 
- Association for Hydraulic Research, Delfi, the Netherlands, Vol. 15, 1977, No. 4, 
Lau, Y. L., and Krishnappan, B. G., ‘‘Transverse 
4 Proceedings of Transport Processes ‘and River Modeling Workshop, he held at the N the National 7 
Water Research Institute, Burlington, Ontario, Canada, in 1978. vee 
26. Lau, Y. L., and Kishnappan, B. G., ‘Modeling Transverse Mixing in Natural Streams,” 
Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY2, Proc. Paper 16048, y 
27. Sayre, W. w., and F. “A Laboratory Investigation of f Open-Channel 
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Page 1609, Eq. 8 should read“. . 


United = 
"Geological Survey 433- United States 1es Geological Su Survey, 


Washington, D. C., 1968, 71 p 


Boome. —The following corrections s should be made to the original paper: 


in instead of 


Q 


Page i613, sth ‘ne should read “adjusted” instead of ‘ ‘adjuusted”” 
1614, Table 1, Column 8: Should read ‘*1.44 -11- 2.1° -11- 3.1° instead 


Page 1614, 14, Sth line fi from Should of reason- 


1617, Table 2, Column 6: Should read * instead of 


Page 1617, 7th line from bottom: Should read ‘ “hydraulic instead « of ‘ ‘hydralic’”’ 
Page 1620, 9th line from top: Should read sina a ‘Single . . 


‘SEDIMENTATION IN 


Discussion by S. Bruk,” ASCE, Vera Miloradov,* and V. — 


"Sedimentation i in ‘Tron Gates Reservoir has been preoccupation of writers 


‘for many years. The first forecast of sedimentation was made in 1962, at an 
early stage of the design, and since then investigations have continued to the | 
present. Several papers have been published on this, but only a few of them i 
in easily accessible international periodicals. The writers wish, therefore, to 
commend the author for bringing up this topic in the ASCE Journal. er 
_ The investigations can be divided into two phases: (1) The first, before the 
construction of the dam, which led to the formulation of the mathematical _ 
- model and a preliminary forecast of reservoir sedimentation (5); and (2) the 
“October, 1980, by Vlad Focsa (Proc. Paper 

Deputy Dir., Jaroslav Cerni Inst., Belgrade, Yugoslavia. 


“Former Head of Dept., Jaroslav Cerni Inst., Belgrade, 
*Asst., Prof., Faculty of Civ. Engrg. Subotica, Yugoslavia. 


| 
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second phase, after ‘completion of the dam, which consists of field ments 
and research work for the improvement of predictions of reservoir sedimentation. > 
The writers particularly wish to discuss this second phase of the research work, _ 
which is partly reflected in the references at the end of this discussion. == ‘ 
‘Suspended Load Measurements.—Suspended load is the main factor of mor 
_ phology of the river reach in question, as can be seen from the considerable x 
overlapping of grain sizes found in the bed material and also carried in suspension. 
Great : attention has, therefore, been | given to developing reliable methods for — 
_ In the first years of the program, standard sampling methods were used with _ 
the usual one liter sampler, as described also by the author. However, the e 
results were rather unsatisfactory; the vertical distribution of sediment concentra- | 
tions was random and irregular, and the amount of sediment in the sample _ 
_ was too small for grain size analysis, which is crucial in sedimentation calculations. ‘ 
Thus, use of the one liter samplers \ was abandoned on the Yugoslav side more _ 
discharge 
in millions of tons of sediment meters persecond 
1974 
1975 


195 


than 15 years ago. A new type of sampler was perio in the fnedae = 
Institute, which has been in use ever since (14-15). It consists of a vacuum 
tank and a suction pipe, with | a nozzle mounted on the weight of the current 
meter. The volume of samples varies between 10 and 14 L, depending on n the 
concentration of sediments (of the order of 20-300 ppm). The sampling is done 
simultaneously with velocity measurements using the current meter, so that 
information is obtained on sediment and water discharge, at the same time. 
This i is } necessary for good correlation. Iti is s worthwhile noting that the vertical 


"close to the expected theoretical cur curve. . The a amount of sediment collected at ¥ 

each point is sufficient to analyze the grain size distribution, so that the transport © 

of suspended sediments can be broken up into grain size classes. By using 

‘this method, good information about the quantity and grain size of suspended — 
BP ore input in the reservoir could be obtained. A comparative study of “7 

f 


different types of suspended sediment samplers, including the one described, 


carried out in 1979 on the Danube, with the Delft Hydraulics 


— 
Tou 
+ 
| 
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DECEMBER EMBER 1961 
Laboratory (13), The results the excellent performance o of the vacuum 
Bc: Sediment Balance 1972-1979.—The input of sediment is monitored by regular. 
_ measurements on the Danube and tributaries, and the output is measured daily 
immediately downstream of the dam. Table 13 summarizes the data obtained: 


7 _ ‘The foregoing figures are in disagreement with findings of the author, the actual 
; percentage of sediments trapped being considerably higher than that given in 
q 


the paper. Table 10 of the paper gives the yearly average discharges of suspended 
load at Orsova, downstream from the dam, for the period 1972-1975. Based — 
on this data, the total amount of sediment transported for the four years in 
question can be calculated. It amounts to 56.6 million tons. From the writers’ - 

. Table 13 it can be seen that for the same period an output of only 12.1 million 

’ tons was obtained by concentration measurements immediately downstream from 

: the dam. The difference can probably be attributed to different methods of _ 


suspended load determination, and only underlines the aad 


applying reliable methods of measurement. 


4 Density of Sediment Deposits.—Since the sediment bimesis is calculated fi, a 

_ weight, and the | morphological changes are given in volume, the density of 

deposited sediment is a very significant parameter. In order to obtain reliable 7 
_ information on the density of reservoir deposits, several hundred samples were : 
taken along the backwater reach by a new instrument specially constructed — 

7 for this purpose (12). A cylinder is pushed into the muddy deposit and the J 


empl frozen in a the cylinder. It is then brought 1 to the surface i in a a solid s state. : -> 


of each layer c can n be determined. The densities obtained vary between 0. 5 and 4-3) 
15, depending mainly on the grain size distribution. The average value of | 
ton/m’ which was used by the author, in the writers opinion, can only serve ciel 
| _as a preliminary estimate. Consolidation with time is of secondary importance, 
; as compared t with the other factors upon which the density of deposits depends. 
Improvement of Mathematical Model. —The exact prediction of how sedimenta- 
tion in the reservoir will progress in time is of significant practical importance. : 
_ Along a length of about 200 k, the flood plain is protected by levees, and 
an increase of water levels, due to the deposited sediment in the river bed, 
= reduce the factor of safety of the levees. Remedial 1 measures consist of 


= Either of ‘these » measures is very expensive thus, the prediction ‘of ie 
sedimentation is of great interest. Apart from the regular program consisting 
q monitoring the sediment input and output, research has been conducted to + 
-.. the mathematical model which is being used for predictions. Be al 


dL Calibration of the sediment transport formula used in the model by field 
2. Investigation of the mobility of sediment deposits in the reservoir. aden 
_ 3. Determination of the settling properties of the suspended sediment by z 
‘direct measurements in the field, with particluar attention given to modification 
of the g grain size distribution of suspended sediments during settling, = 
= Analysis of computational problems by applying appropriate techniques 


to the channel — selection of time apaatiened and reach lengths, oe 
7 


f 
| 

| 
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DISCUSSION 


theoretically more justified approach to sedimentation calculations, by ia 
ming 


of the modern concepts of turbulence modeling (11). ic 
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‘Tora Lew 1 IN 
Discussion William R. Brownlie,‘ A. M. 


at 
i authors have pre presented an interesting dimensional =netia of the sediment 7 
transport process in alluvial channels. Unfortunately, in working with the same 
_ data bank used by the authors, the writer has discovered a significant number 
of errors. Therefore, some of the statistical analyses presented may be inaccurate. 
_ The data bank in question is the compilation of Peterson and Howells (9). 
This first attempt at developing a computerized data base is to be commended. 
However, before it can be used with total satisfaction, all of the errors must 


A careful, item-by-item check suggests that four types of errors were made 


in the preparation ofthe datacompendium (9), 
P ies usually have in incorrectly ore 
i 2. Conversion errors—errors made in converting the data toa standard f ormat, ¥ 


typically involving conversion of transport rates to sediment format, typically 
involving conversion of transport rates to sediment concentrations. = =—/ ~ 
Bac. _Misinterpretation of data—this error usually involved whole columns of 
data, and probably occurred as a result of confusing notation in the data source. 
4. Source errors—errors originating from incorrect publication of data, 
determined by checks on internal consistency. Also encountered were omissions | 
of entries such as bed form and the gradation parameter (geometric standard 
deviation of bed particles), which could be determined, but were not given 
The authors Table 1 indicates that all four types of errors | errors are present in 
the subset of the data used here. The following is a description of some of 
the apparent errors that might be important to the present analyses. In the ~ 
, data of Sato and Kikkawa, authors set number 6, the grain size given in centimeters — 
_was read as millimeters. Therefore, the median sediment size should be 1.0 
-m, and not 0.1 m. The Straub data, set number 15, contains 3 concentration 
values which are a factor of 10 too high. For two sets of data by Vanoni 
and Hwang, numbers 17 and 19, the values given for discharge are really flow 
velocity, and the slope and depth entries are interchanged. An incorrect inter - 
‘ BD aac of the transport rate of the Williams data (set 38) as being given 
in dry unit weight per time instead of submerged weight resulted in an error 
of about 60% in the sediment concentration readings. The transport rate i 
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tons, causing a a 12% error in sediment concentration. 

_ Areview of Table | suggests that perhaps about 10% of the data used jini” 

serious errors. To reiterate, these apparent errors were not made by the authors, < 

at rather in the preparation of the data compendium of Peterson and Howells 

, The exact effect on the statistical analyses is not known. a= Ws = me 
‘= A new data compilation, with a computer format similar ‘to (9) has ia 
prepared at Caltech. Documentation for general release is now underway. ’ The | 
new compilation contains most, but not all, of the data sets in (9), as well 
This discussion is based upon 1 work supported ‘the National Science 
Foundation, under Grant CME79- 20311. 
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